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This review describes recent groundbreaking results in Si, Si/SiGe and dopant-based 
quantum dots, and it highlights the remarkable advances in Si-based quantum physics 
that have occurred in the past few years. This progress has been possible thanks to 
materials development of Si quantum devices, and the physical understanding of quan- 
tum effects in silicon. Recent critical steps include the isolation of single electrons, the 
observation of spin blockade and single-shot read-out of individual electron spins in both 
dopants and gated quantum dots in Si. Each of these results has come with physics that 
was not anticipated from previous work in other material systems. These advances un- 
derline the significant progress towards the realization of spin quantum bits in a material 
with a long spin coherence time, crucial for quantum computation and spintronics. 
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I. INTRODUCTION AND MOTIVATION 
A. Silicon Quantum Electronics 

The exponential progress of microelectronics in the last 
half century has been based on silicon technology. Af- 
ter decades of progress and the incorporation of many 
new materials, the core technological platform for classi- 
cal computation remains based on silicon. At the same 
time, it is becoming increasingly evident that silicon can 
be an excellent host material for an entirely new gen- 
eration of devices, based on the quantum properties of 
charges and spins. These range from quantum comput- 
ers to a wide spectrum of spintronics applications. Sili- 
con is an ideal environment for spins in the solid state, 
due to its weak spin-orbit coupling and the existence of 
isotopes with zero nuclear spin. The prospect of combin- 
ing quantum spin control with the exquisite fabrication 
technology already in place for classical computers has 
encouraged extensive effort in silicon-based quantum de- 
vices over the past decade. 

While there are many proposed physical realizations 



for quantum information processors (Buluta et al. 2011 



Ladd et al.\ |2010| Lloyd 1993), semiconductor-based 
quantum bits (qubits) are extremely interesting, in no 
small part because of their commonalities with classi- 



cal electronics ( Kane 1998 Loss and DiVincenzo 19981. 



Electron spins in quantum dots have received consid- 
erable attention, and significant experimental progress 



has been made since the original Loss and DiVincenzo 



( 1998) proposal. Experiments on lithographically defined 
quantum dots in GaAs/AlGaAs heterostructures have 
shown qubit initialization, single-shot single-electron spin 



single-spin (Koppens et al. 2006) and two-spin (Petta 



et al.\ |2005 



read-out ( |Elzerman et al. 2004 ) , and coherent control of 



states. One of the major issues in Al- 
GaAs/GaAs heterostructures is the inevitable presence 
of nuclear spins in the host material, leading to relatively 
short spin relaxation and coherence times. 

A way to increase the coherence time is to use materi- 
als with a large fraction of non-magnetic nuclei. Natural 
silicon consists of 95% non-magnetic nuclei (92% 28 Si and 
3% 30 Si) and can be purified to nearly 100% zero-nuclear- 
spin isotopes. Various proposals have been made for elec- 



tron spin qubits based on donors in Si (De Sousa et al. 


2004 


Hill et al.\ |2005l |Hollenberg et al. 2006| 


Vrijcn 


et al. 


2000 


1 and Si quantum dots ( 


Friesen et al. 


2003). 



The key requirement for spin quantum bits is to con- 
fine single electrons to either a quantum dot or a donor, 
thus posing a scientific challenge. In contrast with the 
technological maturity of classical field-effect transistors, 
Si quantum dot systems have lagged behind GaAs sys- 
tems, which were historically more advanced because of 
the very early work in epitaxial growth in lattice-matched 



III-V materials. Kouwenhoven et al. (1997b) studied the 



excitation spectra of a single-electron quantum dot in 
a III-V material. Even though Coulomb blockade in Si 



structures was observed very early ( Ali and Ahmed 1994 



Paul et al. 19931, it took another 5 years before regular 



Coulomb oscillations were reported (Simmel et al. 1999). 



Silicon systems needed nearly ten years to achieve single- 



electron occupation in 


quantum dots ( 


Lim et al. 


20091i 


Simmons et al. 




2007 




Zwanenburg et al. 




2009b 1 and 


dopants ( 


Fuechsle et al. 


2012 


Sellier et al. 


2006). For 



quantum dots this has laid the foundation for spin filling 



in valleys in few-electron quantum dots (Borselli et al. 



2011a Lim et al. 20111, tunnel rate measurements in 



few-electron single and double quantum dots (Thalaku- 



lam et al. 


2010 


), Pauli spin 


regime ( 


Borselli et al. 


2011b 



cillations of singlet-triplet states (Maune et al. 2012). In 



the case of dopants valley excited states ( Fuechsle et al. 



2010), gate-induced quantum-confinement transition of 



a single dopant atom (Lansbergen et al. 2008), a deter- 



ministically fabricated single-atom transistor (Fuechsle 



et al. 2012 ) and single-shot read out of an electron spin 



bound to a phosphorus donor (Morello et al. 20101 have 
been reported. The importance of deterministic doping 
has recently been highlighted in the 2011 ITRS Emerg- 



3 



ing Research Materials chapter, where a remaining key 
challenge for scaling CMOS devices towards 10 nm is the 



control of the dopant positions within the channel (ITRS 



2011 ). All these results underline the incredible potential 



of silicon for quantum information processing. 

It is tempting to project the achievements in 
integrated-circuit technology onto a supposed scalabil- 
ity of quantum bits in silicon. Even though current sili- 
con industry standards, with 22 nm features, have higher 
resolution than typical quantum devices discussed in this 
review, superb patterning alone does not guarantee any 
sort of 'quantum CMOS' (Complementary Metal-Oxide- 
Semiconductor). As one example, interface traps have a 
very different effect on classical transistors (where they 
serve as scattering centers or shift threshold voltages) 
than in quantum dots (where they also affect spin co- 
herence). Nonetheless, a fully-integrated CMOS foundry 
has been used for many steps in the fabrication of silicon 



quantum devices (Nordberg et al. 2009a). 



While silicon-based devices generate special interest for 
quantum computation, because of zero nuclear spin iso- 
topes and low spin-orbit coupling, they also face some 
special challenges and display physics that, until recently, 
has been little explored in the context of quantum com- 
putation. Examples of the challenges include the rela- 
tively large effective mass in silicon and the large differ- 
ence in lattice constant between silicon and germanium. 
An example of the unexplored physics is the presence of 
multiple conduction band valleys in silicon. 

As described in this review, there have been rapid ad- 
vances addressing the challenges and exploring the new 
physics available in silicon-based quantum devices. The 
extent to which these advances will lead to larger-scale 
quantum systems in silicon is an exciting question as of 
this writing. 



B. Outline of this review 

This review covers the field of electronic transport in 
silicon and focuses on single-electron tunneling through 
quantum dots and dopants. We restrict ourselves to 
experiments and theory involving electrons confined to 
single or double (dopant) quantum dots, describing the 
development from the observation of Coulomb blockade 
to single-electron quantum dots and single dopant atom 
transistors. Ensembles of quantum dots or dopants are 
beyond the scope of this article. Also, the review is 
strictly limited to electron transport experiments, and 
does not cover optical spectroscopy measurements. Op- 
tical spectroscopy on quantum dots and ensembles of 
dopants is a very active and emerging field, see for exam- 



ple the recent work by Greenland et al. (20101 and Steger 



et al. (20121 and references therein. 

Section [II] Quantum Confinement starts with a gen- 
eral introduction to transport through quantum-confined 



silicon nanostructures. The silicon bandstructure is de- 
scribed in Section Physics of Silicon nanostructures 
with specifics such as the valley degeneracy and split- 
ting in bulk and quantum dots, and wave function con- 
trol and engineering of dopant states. Section \TV\ Quan- 
tum dots in Si and SiGe explains the development from 
the discovery of Coulomb blockade in 1990 to single- 
electron occupancy in single and double quantum dots 
in recent years. Analogously, dopant transport in sil- 
icon has evolved from tunneling through 1980's MOS- 
FETs to current-day single-atom transistors, see Section 
\V\ Dopants in silicon. The remarkable advances of Sec- 
tions |IV] and |V] have lead to the relaxation and coherence 
measurements on single spins in Section |V7| Outlook: re- 
laxation, coherence and measurements. 



II. QUANTUM CONFINEMENT 

This section introduces quantum electronic experi- 
ments in silicon, starting with the quantum mechanical 
confinement of electrons in silicon, which can be achieved 
by a combination of electrostatic fields, interfaces be- 
tween materials, and/or placement of individual atoms. 
All of these approaches lead to single-electron tunneling 
devices consisting of a silicon potential well coupled to 
source, drain and gate electrodes. 



A. From single atoms to quantum wells 

Electrons in Si nanostructures are confined using a 
combination of material and electrostatic potentials. 
The shape and size of nanostructured materials provide 
natural confinement of electrons to 0, 1 or 2 dimensions. 
The exact confinement potential of the structure in x, 
y and z-directions sets the additional requirements in 
terms of additional electric fields. Figure [T] gives an 
overview of materials of different dimensionality and 
their integration into single-electron tunneling devices. 

Dopants 

The electrostatic potential of a single dopant atom is 
radially symmetric, resulting in the same steep potential 
well in all directions, as shown in the first row of Fig. [T] 
The Bohr radius as is the mean radius of the orbit of 
an electron around the nucleus of an atom in its ground 
state, and equals for example 2.5 nm for phosphorus 
in silicon. A dopant atom has three charge states: the 
ionized D + state, the neutral D° state (one electron 
bound to the dopant) and the negatively charged D~ 
state (two electrons bound to the dopant). Because 
the D + state corresponds to an empty dopant it does 
not appear as an electron state in the potential well. 
Measuring electron transport through a single atom has 
been a great challenge, as described in Section [Vj but 
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FIG. 1 (Color online) Combining material and electrostatic confinement to create single-electron transistors. 

First column: schematic of dopants, 0D-, ID- and 2D-structures. Second column: in the corresponding confinement potentials 
in x-, y- and z-directions electron states are occupied up to the Fermi energy Ef (dashed grey line). Occupied and unoccupied 
electron states are indicated as straight and dashed lines respectively. Third column: Schematic of the silicon nanostructure 
integrated into a transport device with source, drain and gate electrodes. Fourth column: The potential landscape of the single- 
electron transistor is made up of a potential well which is tunnel-coupled to source and drain reservoir and electrostatically 
coupled to gates which can move the ladder of electrochemical potentials, as described in Section [II. B| 



the single-dopant regime as sketched in the third column 
has been reached by several groups. Depending on the 
architecture, the source and drain reservoirs can be 



made up of highly-doped Si (Fuechsle et al. 2012 Pierre 



et al. 2010 Sellier et al. 2006), or of a two-dimensional 



electron gas (Tan et al. 2010). The same goes for the 



gates, but they can also be metallic (Tan et al. 2010). 



The resulting single-electron transistors consist of a 
steep dopant potential well connected to source and 
drain reservoirs. 



OD structures 

Like dopants, self-assembled nanocrystals provide con- 
finement to zero dimensions, but the confinement is 
better described by a hard-wall potential well in x, y 



and z-directions and is much wider (Fig. [I]). The energy 
levels of an electron in a quantum well of size L are 
quantized according to basic quantum mechanics, see 



for example Cohen- Tannoudji et al. (1992). The corre- 
sponding level spacing AE is on the order of h 2 /m e eL 2 , 
where m e ff is the electron effective mass. The separation 
between energy levels thus decreases quadratically with 
the well width: as a result, the discrete levels of e.g. 
a 30 nm size nanocrystal are expected to have energy 
spacings 2 orders of magnitude smaller than those of a 
dopant with a 3 nm Bohr radius. Making source and 
drain contacts requires very precise alignment by means 
of electron-beam lithography. The tunnel coupling 
of these devices relies on statistics; creating tunable 
tunnel coupling to self- assembled dots is very challeng- 
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ing. A highly-doped substrate can be used as a global 
backgate and metallic leads on a dielectric as a local gate. 

ID structures 

The high aspect ratio of nanowires (NWs) implies a large 
level spacing in the transverse directions, and a small 
level spacing in the longitudinal direction [L x 3> L y ^ z ), 
creating a (quasi) 1-dimcnsional channel with few 
subbands in the transverse direction (see second row 
of Fig. [TJ . Within this channel a zero-dimensional well 
can be created by local gates on the nanowire, or by 
Schottky tunnel barriers to source and drain contacts. 
In the latter case the barrier height is determined by the 
material work functions and hardly tuneable in-situ - 
the tunnel coupling will generally decrease as electrons 
leave the well and the wave function overlap with source 
and drain shrinks. Local gates, however, can tune the 
tunnel barriers since the applied gate voltage induces an 
electric field which locally pulls up the conduction band. 
Electrons tunnel from the quantum well into reservoirs 
which are part of the nanowire itself. The metallic 
leads connecting the nanowire to the macroscopic world 
must be ohmic; i.e., the contacts should have high 
transparency, to prevent the formation of multiple 
quantum dots in series (particularly if the contacts are 
very close to the quantum dot). 

2D structures 

A 2-dimensional electron gas (2DEG) can be created in 
Si MOSFETs (Metal-Oxide-Semiconductor Field Effect 
Transistors) and in Si/SiGe heterostructures. Electrons 
are unconfincd in the x-y-plane and are confined by a 
triangular potential well perpendicular to the plane, as 
sketched in in Fig. [TJ More realistic band diagrams are 
drawn in Fig. 2 in the review by |Ando et al. ( 1982 ) for 
Si MOS and Fig. 11 in the review by |Schaffler| (|1997| 



for Si/SiGe heterostructures. In a 2DEG-based quantum 
dot, the lateral confinement is a soft-wall potential de- 
fined by top gate electrodes, enabling tunnel-coupling to 
source and drain reservoirs in the 2DEG. Those reservoirs 
are connected to macroscopic wires via ohmic contacts, 
which are often highly doped regions at the edge of the 
chip. The resulting potential landscape is highly tunable 
thanks to local electrostatic gating via the top gates. 



B. Transport regimes 

Having introduced quantum-confined devices, we now 
cover the basics of quantum transport through single- 
electron transistors (SETs), which are made up of a zero- 
dimensional island, source and drain reservoirs, and gate 
electrodes. 

Electronic measurements on single electrons require a 
confining potential which is tunnel coupled to electron 
reservoirs in source and drain leads, see Fig. [2j The 




E r +AE 



FIG. 2 (Color online) Schematic diagrams of the elec- 
trochemical potential of a single-electron transistor. 

(a) There is no available level in the bias window between us 
and fiD, the electrochemical potentials of the source and the 
drain, so the electron number is fixed at N due to Coulomb 
blockade, (b) The /ijv level aligns with source and drain elec- 
trochemical potentials, and the number of electrons alternates 
between TV and N — 1, resulting in a single-electron tunneling 
current. 



SET-island is also coupled capacitively to one or more 
gate electrodes, which can be used to tune the electro- 
static potential of the well. The discrete levels are spaced 
by the addition energy E a ^(N) = Ec + AE, which con- 
sists of a purely electrostatic part, the charging energy 
Ec, plus the energy spacing between two discrete quan- 
tum levels, AE. AE is zero when two consecutive elec- 
trons are added to the same spin-degenerate level. The 
charging energy Ec = e 2 /2C, where C is the sum of all 
capacitances to the SET-islancQ 

In the limit of low temperature, if we only consider se- 
quential tunneling processes, energy conservation needs 
to be satisfied for transport to occur. The electrochem- 
ical potential ^in is the energy required for adding the 
iVth electron to the island. Electrons can only tunnel 
through the SET when /i^r falls within the bias window 
(see Fig. [2](b)), i.e. when ^5 > > [i£>. Here fis 
and hd are the electrochemical potential of the source 
and the drain respectively. Current cannot flow without 
an available level in the bias window, and the device is 
in Coulomb blockade, see Fig. [5Ja). A gate voltage can 
shift the whole ladder of electrochemical potential levels 
up or down, and thus switch the device from Coulomb 
blockade to single-electron tunneling mode. By sweep- 
ing the gate voltage and measuring the conductance, one 
obtains Coulomb peaks as shown in Fig. [3^ a). 

Usually, one measures the conductance versus source- 
drain voltage Vsd and gate voltage Vg in a bias spec- 
troscopy, as shown in Fig. [3](b) . Inside the diamond- 
shaped regions, the current is blocked and the number of 
electrons is constant. At the edges of these Coulomb dia- 
monds a level is resonant with either source or drain and 



We refer to other review articles on quantum dots and single- 



electron transistors for more background and details: Bccnakkcr 
and van Houten| (|199H; |Grabert et al. | |1993[l; |Hanson et al. 
(|2007E;|Kouwenhoven et aL|||2001||1997a ;|Reimann and Manni- 



|l |2002[ ); and |Van der Wiel et aZ.| | |2003 



G 




FIG. 3 (Color online) Zero-bias and finite-bias spec- 
troscopy, (a) Zero-bias conductance G of transport versus 
gate voltage Vg both at T > T K (solid line) and T < T K 
(dashed line). In the first regime, the full width at half 
maximum (FWHM) of the Coulomb peaks corresponds to 
the level broadening ZiF. In the Kondo regime (T <C Tic), 
Coulomb blockade is overcome by coherent second-order tun- 
neling processes (see main text), (b) Stability diagram show- 
ing Coulomb diamonds in differential conductance, dl/dVsD, 
versus eVsD and eVa at T = OK. The edges of the diamond- 
shaped regions (red) correspond to the onset of current. Di- 
agonal lines of increased conductance emanating from the 
diamonds (gray) indicate transport through excited states. 
The indicated internal energy scales Ec, AE, hT and Tk de- 
fine the boundaries between different transport regimes. Co- 
tunneling lines can appear when the applied bias exceeds AE 
(see main text). Adapted from Lansbergen 2010 



single-electron tunneling occurs. When an excited state 
enters the bias window a line of increased conductance 
can appear parallel to the diamond edges. These resonant 
tunneling features have other possible physical origins, as 



described in detail by |Escott et al.\ ([2010]) . From such a 
bias spectroscopy one can read off the excited-states and 
the charging energy directly, as indicated in Fig. [3^b). 

The simple model described above explains success- 
fully how quantization of charge and energy leads to ef- 
fects like Coulomb blockade and Coulomb oscillations. 
Nevertheless, it is too simplified in many respects. Up 
until now we only worried about the electronic proper- 
ties of the localized state but not about the physics of 
the electron transport through that state. In this sec- 



tion, based on |Lansbergen| [20101 we will describe the five 
different regimes of electron transport through a local- 
ized stated in a three-terminal-geometry. How electrons 
traverse a quantum device is strongly dependent on the 
coherence during the tunneling process and thus depends 
strongly on eVsD and fe^T. These external energy scales 



should be compared to the internal energy scales of the 
tunneling geometry that determine the transport regime, 
namely the charging energy Ec, the level spacing AE, 
the level broadening hT and the Kondo temperature Tk- 
Here, T is the total tunnel rate to the localized state 
which can be separated into the tunnel coupling to the 
source electrode Ts and to the drain electrode Frj, i.e. 
r = r,j + To . The internal energy scales are all fixed by 
the confinement potential, and the external energy scales 
reflect the external environment, namely the temperature 
T and the applied bias Vsd- 

Much literature describes the electronic transport in 
all possible proportionalities of these energy scales with 



1988) 



each other (Alhassid 2000 Beenakker 1991 Buttiker 



The internal energy scales are typically related to 
each other by Tk /ir <C AE <C Ec, and occasion- 
ally by T K < AE < hT < E c , limiting the number 
of separate transport regimes that we need to consider. 
Fig.|3Ja) is a schematic depiction of transport regimes as 
a function of eVsD and k B T . It should be noted that 
the boundaries between transport regimes are typically 
not abrupt transitions. For clarity, internal and exter- 
nal energy scales (except Tk and hT) are indicated in a 
schematic representation of our geometry, see Fig.[4]Jb). 

Here, we will not make a distinction between the ex- 
ternal energy scales fcgT and eVsD when we compare 
them to internal energy scales, as indicated by Fig.^a). 
The reason behind this equality is that both these exter- 
nal energy scales have a very similar effect on the trans- 
port characteristics. Their only relevant effect is that 
they introduce (hot) phonons to the crystal lattice, ei- 
ther directly by temperature or by inelastic tunneling 
processes induced by the non-equilibrium Fermi energies 
of the source/drain contacts. 

Next, we will describe the five separate tunneling 
regimes and their corresponding expressions for the 
source/drain current / shortly. These regimes are the so- 
called multi-electron regime, the sequential multi-level 
regime, the sequential single level regime, the coherent 
regime and the Kondo regime, see Fig. (4^a). 



1. The multi-electron regime 

Firstly there is the multi-electron regime {Ec <C 
ksT, eVsD) where Coulomb blockade does not occur, as 
mentioned in the start of this chapter. This regime is not 
relevant for this review. 



2. The sequential multi-level regime 

At AE <C k B T,eVsD -C Ec the system is in the se- 
quential multi-level regime. The transport is given by 
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FIG. 4 (Color online) The five separate transport 
regimes in a three-terminal quantum device. (a) 

Schematic depiction of the regimes in which transport through 
a localized takes place as a function of the external energy 
scales ksT and Vsd- The transitions between regimes take 
place on the order of the internal energy scales Ec, AE, hT 
and Tk- (b) Potential landscape of the three terminal ge- 
ometry, where the quantum states and the electrochemical 
potential of the leads are shown together with ksT, Vsd and 
E c , AE. 



laxation prohibits coherent transfer of electrons from the 
source to the drain electrode. 



3. The sequential single-level regime 

The next transport regime is the sequential single-level 
regime, roughly bounded by hT -C fcsT, eVsD "C AE, 
where only a single level resides inside the bias window. 
This regime is a transition between phase-coherent and 
phase-incoherent transport between source- and drain 
-electrodes and the tunneling current depends vitally 
on ksT. For Vsd = the conductance is given by 



(Beenakker 1991) 



G = 



r 1 r 1 

in out 



- F 1 : 

1 out 



(2) 



where T m is the tunnel rate into the localized state and 
r out is the tunnel rate out. Note that r; n = r,s,r out = 
T D for V SD > and r in = T D ,T out = T s for V SD < 0. 

If the localized state is strongly coupled to the contacts 
higher-order transport processes become apparent in the 
Coulomb blocked region, i.e. the so called co-tunneling 
lines indicated in Fig.[3]jb). This is the case when Ec /T 
approaches unity in the open regime. There is an elastic 



and inelastic component to the co-tunneling (Averin and 



Nazarov 1990 Nazarov and Blanter 2009). The elastic 



component leads to a constant background current in the 
Coulomb diamond. The inelastic component leads to a 
step in the current when the applied bias exceeds AE. 
The current is given by 



( |Beenakker[|1991| |Van der Vaart et aZ.[|1993 l 



/ = e- 



(Yl + T 2 

\ i ti 1 lr 



111/ out 



rt 



•r?„ 



pn 



F 1 ' 

L out 



(1) 



where the subscript denotes the direction of transport, 
into or out of the localized state, and the superscript in- 
dicates the level, where 1 refers to the ground state and n 
indicates the highest orbital within the energy window set 
by eVsD- The current thus depends on the ingoing rates 
of all levels in the bias window and the outgoing rate of 
only the ground state. Physically, electrons can enter any 
orbital state that is energetically allowed. Once a single 
electron is transferred to the localized state, Coulomb 
blockade prevents another electron from entering. For 
dopants, the bound electron will relax back to the ground 
state before it has a chance to tunnel out of the localized 
state, since the orbital relaxation times (~ ps-ns ( Lans- 



bergen et al. 2011 )) are typically much faster than the 
outgoing tunnel rates (~1 ns). For quantum dots the 
physics is similar but tunnel rates and orbital relaxation 
rates are slower, e.g. ~ l-10ns'in GaAs quantum dots 
(Fujisawa et al. 19981. The inelastic nature of the re- 
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(3) 



(4) 



for the elastic and inelastic co-tunneling respectively with 
E e + Eh = Ec, where the energies E e and Eh denote the 
distance to the Fermi energy of the filled and empty state 
and p is the density of states. The complex co-tunneling 
line shape is discussed in depth in Wcgcwi js and Nazarov[ 
I2TKJT1 



4. The coherent regime 

As soon as the external energy scales are much smaller 
then hT (T K < k B T,eV SD < hT < AE) the system is 
in the coherent regime, where the conductance is given 



by Buttiker (1988) 



G 



in out 



(r 1 

V in 



-^out ) 



(5) 



The conductance is thus given by the quantum conduc- 
tance e 2 /h multiplied by a factor that only depends on 
the symmetry between T$ and . It has been proven ex- 
plicitly that this expression, easily derived for resonances 



in ID double barrier structures (Ricco and Azbel 19841, 



also holds in three dimensions (Kalmeyer and Laughlin 



1987). 



logical relationship ( Goldhaber-Gordon et al. 1998 ) 



G(T) = (G) 



T. 



'2 



K 



T 2 + T >2 



(7) 



where T K = Tk/V^ 1 ^ — 1, Go is the zero-temperature 
Kondo conductance and s is a constant found to be equal 



to 0.22 (Goldhaber-Gordon et al. 1998) 



5. The Kondo regime 



III. PHYSICS OF SI NANOSTRUCTURES 



The final transport regime occurs when eVsD, ksT <C 
Tk- The Kondo temperature is the energy scale be- 
low which second-order charge transitions other than co- 



tunneling start to play a role in the transport ( Meir and 



Wingreen 1993). In first-order transitions, the trans- 



ferred electrons make a direct transition from their initial 
to their final state. It should be noted that the constant 
interaction model only considers first-order charge tran- 



sitions (Kouwenhoven et al. 1997a). In a second-order 



transition, the transferred electron goes from the initial 
to the final state via a virtual state of the atom or dot. A 
virtual state is an electronic state for which the number 
operator does not commute with the Hamiltonian of the 
system and therefore has a finite lifetime. The lifetime of 
the virtual state is related to the Heisenberg uncertainty 
principle, as the electron can only reside on the virtual 
state on a timescale t ~ hj (fiN — hs.d), where ^n — ^s,d 
is the energy difference between the virtual state and the 
nearest real state. The main characteristic of this trans- 
port regime is a zero-bias resonance inside the Coulomb 
diamond for 7V=odd, as we will explain next, see also 
Fig.§a) and (b). 

When N — even, the total localized spin is zero due to 
the (typical) even-odd filling of the (spin) states, result- 
ing in zero localized magnetic moment. When TV = odd, 
one electron is unpaired, giving the localized state a net 
magnetic moment. In contrast to metals doped with 
magnetic impurities, the conductance of double barrier 
structures actually increases due the Kondo effect. This 
is because the density of states in the channel at a [is, /i_d 
(associated with the newly formed Kondo singlet state) 
acts as a transport channel for electrons, as if it were 
a "regular" localized state in the channel. The Kondo 



temperature can be expressed as (Glazman and Pustil 



nik 2003) 



T K = \] E C T exp(- 



2T 



(G) 



assuming hn — hs,d <C /ijv-i ~ Hs.d- The zero-bias 
Kondo resonance is furthermore characterized by its tem- 
perature and magnetic field dependence. The conduc- 
tance of the Kondo resonance has a logarithmic temper- 
ature dependence, which is described by the phenomeno- 



Here we describe the fundamental physical properties 
of Si nanostructures. Some of these arise from the elec- 
tron confinement into a small region (tens of nanometers 
or less) and are similar to those of other semiconduc- 
tors, but other properties are present only in Si. One 
example arises because Si has multiple degenerate val- 
leys in its conduction band, described in the first section. 
The valleys play an important role in both dopant and 
quantum dot devices, although the details of the valley 
physics in those two systems are different. Moreover, in 
heterostructures, strain often plays an important role, 
and the interplay between strain, disorder, and the prop- 
erties of the valleys are important in determining the 
low-energy properties of the devices. 



A. Bulk silicon: valley degeneracy 

Because silicon is used in many technical applications, 
methods for manufacturing extremely high purity sam- 
ples are well-developed. Silicon has several stable nuclear 
isotopes, with 28 Si, which has no nuclear spin, being the 
most abundant (its abundance in natural silicon is 92%). 
This availability of a spin-zero silicon isotope is useful 
for applications in which one wishes to preserve the co- 
herence of electron spins, since the absence of hyperfine 
interaction eliminates a possible decoherence channel for 
the electron spin, see section VI. A. 4. 

The properties of electrons in silicon have been stud- 



ied in great detail for many decades (Cohen and Che- 



likowsky 1988 Yu and Cardona 20011. Here we review 



aspects of the material that will prove critical in under- 
standing the challenges that arise as one works to create 
devices with desired properties on the nanoscale. One 
such aspect is how the effects of multiple valleys present 
in the conduction band in bulk silicon appear in specific 
silicon nanodevices. The manifestations of valley physics 
in quantum dots are different from those in dopant-based 
devices, and understanding the relevant effects is critical 
for manipulating the spin degrees of freedom of the elec- 
trons in nanodevices. In the following subsections, we 
first define and discuss the conduction band valleys in 
bulk silicon, and then the behavior and consequences of 
valley physics for quantum dots and for dopant devices. 
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Si 





FIG. 5 Silicon crystal in real and reciprocal space, (a) 

3D plot of the unit cell of the bulk silicon crystal in real space, 
showing the diamond or Face-Centered Cubic lattice, which 
has cubic symmetry, (b) Silicon crystal in reciprocal space. 
Brillouin zone of the silicon crystal lattice. It is the Wigner- 
Seitz cell of the Body-Centered Cubic lattice. F is the center 
of the polyhedron. Figure from Davies (19981. 




FIG. 6 (Color online) Band structure of bulk silicon, (a) 

The conduction band has six degenerate minima or valleys at 
0.85fco. Results kindly supplied by G.P. Srivastava, University 
of Exeter. Figure from Davies (19981. (b) Zoom-in on the 



bottom of the conduction band and the top of the valence 
band (schematic, not exact). The bandgap in bulk Si is 1.12 
eV at room temperature, increasing to 1.17 eV at 4 K ( |Green 
1990[ ) . The heavy and light hole bands are degenerate for k = 
0. The split-off band is separated from the other subbands by 
the spin-orbit splitting A so of 44 meV. 



Crystalline silicon is a covalently bonded crystal with a 
diamond lattice structure, as shown in Fig. 5. The band 



structure of bulk silicon (Phillips 1962 1, shown in Fig. 6, 



has the property that the energies of electron states in 
the conduction band is not minimized when the crystal 
momentum k — 0, but rather at a nonzero value, ko, 
that is 85% of the way to the Brillouin zone boundary, 
as shown in Fig. 6(b). Bulk silicon has cubic symmetry, 
and there are six equivalent minima. Thus we say that 
bulk silicon has six degenerate valleys in its conduction 
band. 

In conventional electronic devices, the presence of mul- 
tiple valleys typically does not affect transport proper- 
ties in a profound way. However, valley physics plays a 
critical role in quantum electronics because of interfer- 
ence between different valleys that arises when the elec- 
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FIG. 7 (Color online) Valley splitting of of dopants and 
of quantum dots in silicon quantum wells, (a) For a 

quantum well, in which a thin silicon layer is sandwiched be- 
tween two layers of Si^Gei-^, with x typically ~ 0.25 — 0.3, 
the six-fold valley degeneracy of bulk silicon is broken by the 
large in-plane tensile strain in the quantum well so that two 
T-levels are about 200 meV below the four A-levels ( Schaifler 



et al. 



1992 1. The remaining two- fold degeneracy is broken by 
the confinement in the quantum well and by electric fields, 
with the resulting valley splitting typically ~ 0.1 — 1 meV. 
(b) For phosphorus dopants, strong central-cell corrections 
near the dopant break the six-fold valley degeneracy of bulk 
silicon so that the lowest-energy valley state is non-degenerate 
(except for spin degeneracy), lowered by an energy 11.7 meV. 
The degeneracies of higher-energy levels are broken by lattice 
strain and by electric helds. 



tronic transport is fundamentally quantum. For example, 
the presence of an additional valley greatly complicates 
spin manipulation because it can lift Pauli spin block- 
ade, which is fundamental for many strategies for spin 



manipulation in quantum dot nanodevices ( Hiittel et al. 


2003 Johnson et al. 


2005a 


Koppens et al.\ 


2005 Ono| 


et al. 


2002 


Rokhinson et al 


2001 


). In pure bulk sili- 



con, the valleys are degenerate (the energies of the six 
states related by the cubic symmetry are the same), but 
in nanodevices this degeneracy can be and usually is bro- 
ken by various effects that include strain, confinement, 
and electric fields. When valley degeneracy is lifted, at 
low temperatures the carriers populate only the lowest- 
energy valley state, thus eliminating some of the quantum 
effects that arise when the valleys are degenerate. 

Fig. [7] shows a summary of valley splitting in het- 
crostructures and in dopant devices. For strained silicon 
quantum wells, the large in-plane strain lifts the energies 
of the in-plane (x and y) valleys. The remaining two- 
fold degeneracy of the z-valleys is broken by electronic 
z-confinement induced by electric fields and by the quan- 
tum well itself, resulting in a valley splitting of order 
0.1 — 1 meV. The breaking of the two-fold valley degen- 
eracy is very sensitive to atomic-scale details of the in- 
terface, and is discussed in detail in Sec. |IH.B] and in the 
supplemental material. 

For an electron bound to a dopant in silicon, the valley 
degeneracy of bulk silicon is lifted because of the strong 



confinement potential from the dopant atom (Kohn and 



Luttinger 1955a). For phosphorus donors in silicon, the 
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electronic ground state is non-degenerate, with an energy 
gap of ~ 11.7 meV between the non-degenerate ground 



state and the excited states ( Andresen et al. 


2009 


Ram- 


das and Rodriguez 


1981 


). Thus, additional degeneracy 



of the electronic ground state is not a concern in dopant 
devices. However, the fact that the conduction band min- 
imum in silicon is at a large crystal momentum k$ that 
is near the zone boundary gives rise to other physical ef- 
fects that are important for quantum electronic devices. 
One such consequence arises because the wave functions 
of the electronic states in dopants oscillate in space on 
the very short length scale ~ 27r/fcn, which is roughly 
on the scale of one nanometer. These charge oscillations 
differ from the electron charge variations due to Bloch os- 
cillations because they can cause the exchange coupling 
to change sign, and thus have significant implications for 
the design of quantum electronic devices, as discussed in 
Section ITlLCl 



B. Quantum wells and dots 

In the quantum well devices we discuss here, one starts 
with a material with a two-dimensional electron gas 
(2DEG), and then lithographically patterns top gates to 
which voltages are applied that deplete the 2DEG sur- 
rounding the quantum dot. By carefully adjusting the 
gate voltages, one can achieve dots with occupancy of 



a single electron, see section IV. D Moreover, the same 
gate voltages that are used to define the dot are also 
used to perform the manipulations required for initial- 
ization, gate operations, and readout of charge and spin 
states (Maune et al. 2012), see section VI. C. 4. 



1. Valley splitting in quantum dots 

Understanding the valley degrees of freedom is impor- 
tant for ensuring that the valley splitting is in a regime 
suitable for spin-based quantum computation. Even in 
the low-density limit appropriate to single-electron quan- 
tum dots, where electron-electron interactions (Ando 



et al. 1982) are unimportant, valley splitting is complex: 
the breaking of the valley degeneracy involves physics on 
the atomic scale, orders of magnitude smaller than the 
quantum dot itself, so it depends on the detailed prop- 
erties of alloy and interface disorder. Because the loca- 
tions of the individual atoms in a given device are not 
known, statistical approaches to atomistic device mod- 
eling or averaging theories like effective mass must be 
utilized. Theory, modeling, and simulation provide in- 
sight into the physical mechanisms giving rise to valley 
splitting, so that device design and fabrication methods 
can be developed to yield dots with valley splitting com- 
patible with use in spin-based quantum information pro- 
cessing devices. 



In bulk silicon, there are six degenerate conduction 
band minima in the Brillouin zone (valleys) as depicted 
in Fig. 5. One modern strategy for fabricating Si de- 
vices for quantum electronics applications is to use a bi- 
axially strained thin film of Si grown on a pseudomor- 
phic Si x Gei_ x substrate. In such devices, the silicon 
quantum well is under large tensile strain, and the six- 



fold degeneracy is broken into a two- fold one (Schafflcr 



1997). Confinement of electrons in the z-direction in a 2- 
dimcnsional electron gas lifts the remaining two-fold val- 
ley degeneracy, resulting in four A-valleys with a heavy 
effective mass parallel to the interface at an energy sev- 



eral tens of meV above the two T-valleys ( Ando et al. 
1982), as shown in Fig. [7j The sharp and flat interface 



produces a potential step in the z-direction and can lift 
the degeneracy of the T-valleys in two levels separated 
by the valley splitting Ey ■ Built-in or externally applied 
electric fields break the symmetry of the Hamiltonian and 
can couple the various valleys and thus lift the valley de- 
generacy. Theoretical predictions for the valley splitting 
of flat interfaces are generally on the order of 0.1-0.3 



meV (Boykin et al. 



and Uemura 1977 



2004b Culcer et al. 



2010a Ohkawa 



Saraiva et al. 20111. Experimen- 



tal values in Si inversion layers mostly vary from 0.3-1.2 



meV, but some are substantially smaller 


(Koester et al. 


1997 Kohler and Roos 


1979 


Lai et al. 


2006; Nicholas| 


et al. 


1980 Pudalov et al. 


1985 


Weitz et al. 


1996 


). A 



giant valley splitting of 23 meV measured in a similar 



understood theoretically (Saraiva et al. 2011) 



structure ( Takashina et al. 2006 1 is still not completely 



The two main approaches for understanding valley 
splitting in silicon heterostructures are tight-binding cal- 



culations ( 


3oykin et al. 


2007 


2004a 2005 


Sharche 


et al. 2007 


Srinivasan et al. 1 


0081 and theories that 


use an effective mass formalism 


(Friesen et al. 




2007b 


Friesen and Coppersmith 


2010 


Saraiva et al 




2009). 



Section I in the supplemental material reviews a sim- 



ple one-dimensional tight-binding model (Boykin et al. 



2004b ) that illustrates some of the physical mechanisms 
that lead the breaking of the valley degeneracy and hence 
the emergence of valley splitting. A pictorial sketch of 
the two lowest-energy eigenstates of this one-dimensional 
model is presented in Fig. [8j The eigenfunctions have 
very similar envelopes and fast oscillations with a period 
very close to 27r/fcn, where fcn is the wavevector of the con- 
duction band valley minimum. The different alignments 
of the phases of the fast oscillations with sharp interfaces 
cause the energies of the two states to be different, thus 
giving rise to valley splitting. 

Valley splitting has a complicated dependence on envi- 
ronmental and structural conditions. Large-scale atom- 
istic tight-binding calculations can incorporate realistic 
inhomogeneity in the atomic arrangement, both in terms 
of alloy disorder and in terms of disorder in the loca- 
tions of interface steps, as discussed in section III of the 
supplemental material. Technically well controlled inter- 
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FIG. 8 (Color online) Sketch of the two lowest energy 
eigenstates in an infinite square well of the two-band 
model presented in the supplemental material. The 

envelopes of the two eigenfunctions are very similar to each 
other and to the sine behavior obtained in the absence of val- 
ley degeneracy; the effects of the valley degeneracy give rise 
to fast oscillations within this envelope. For a square well, 
one eigenfunction is symmetric and the other is antisymmet- 
ric; the symmetries are different because the fast oscillations 
have different phases, as measured from the quantum well 
boundaries. This sensitive dependence of valley splitting on 
the atomic-scale physics near the well boundary is the source 
of the sensitive dependence of the valley splitting on disorder 
at the quantum well interfaces. 



faces in Si are buffers of either SiC>2 or Si x Gei_ x , which 
are intrinsically atomistically disordered. Some of the 
effects of this disorder can be understood qualitatively 
using effective mass theory, but because of the impor- 
tance of atomic-scale physics in determining valley split- 
ting, atom-scale theory is required for quantitative un- 
derstanding. For Si x Gei_ x , there are 3 critical disorder 
effects to consider: atom-type disorder, atom-position 
disorder, alloy concentration disorder. A detailed discus- 
sion of the characterization of the effects of these different 
types is presented in section III of the supplemental ma- 
terial. 

Many features of the physics that give rise to val- 
ley splitting can be understood qualitatively and semi- 



quantitatively using effective mass theories (Kohn and 



Luttinger 


1955b 


Seitz and Turnbull 


1957 



ories are formulated carefully to incorporate the micro- 



scopic effects that give rise to valley splitting (Friesen 



et al. 



2006 



20051 JFriesen et al.\ |2007b| |Fritzsche[ [19621 |Nestoklon| 



Pantelides 19781. In the envelope func- 



tion or effective mass formalism, the theory is written 
in terms of an envelope function for the wave function, 
which is well-suited for describing variations on relatively 
long scales (such as the quantum dot confinement). The 
effects of the degenerate valleys are incorporated using a 
valley coupling parameter that is treated as a delta func- 
tion whose strength is determined by the atomic scale 



physics ( Chutia et al. 


2008 


Friesen et al. 


2007b 


Saraiva 


et al. 


2009 


). The envelope function formalism has the 



advantage that one can obtain analytic results for valley 



splitting in nontrivial geometries ( 


Culcer et al. 


2010a|b 


Friesen et al. 


2007b 


Friesen and Coppersmith 


2010). 



However, the theory must explicitly incorporate infor- 
mation from the atomic scale, either as a valley coupling 
parameter that is fit to tight-binding results, as the out- 



put of a multiscale approach ( Chutia et al. 2008 Saraiva 



et al. 2009), or by explicit atomistic calculation on large 



scales, as embodied by the NEMO tool suite (Boykin 



TTal\ |2004b| |Klimeck et~al\ [2007j [20021 |Steiger et al. 
2011 ). More details of effective mass theory treatment of 



valley splitting are in the supplemental material. 



2. Mixing of valleys and orbits 

When the valley splitting Ey is much greater than the 
orbital level spacing AE, electrons will occupy single- 
particle levels with orbital numbers 1, 2, 3, ... and val- 
ley number VI, the lowest valley state (see Fig. 9(a)). 
Conversely, if AE ^ Ey the first four electrons will 
occupy the valleys VI and V2 in the lowest orbit be- 
fore going to the next orbit with n = 2, as shown in 
Fig. 9(b). However, valleys and orbits can also hybridize 



(Friesen and Coppersmith, 2010), making it inappropri- 



ate to define distinct orbital and valley quantum num- 
bers (see Fig. 9(c)). Depending on the degree of mixing, 
the valley-orbit levels VOl, V02 etc, behave mostly like 
valleys or like orbits. Instead of referring to a pure val- 
ley splitting Ey the term valley-orbit splitting is used, 
Eyo — Eyo2 — Eyoi f° r the difference in energy be- 
tween the first two single-particle levels, Eyoi arl d Eyoi- 



This is referred to as the ground-state gap (Friesen and 



Coppersmith 2010) 



The behavior of the valley splitting in real quantum 
wells is complicated by the fact that in real devices the 
quantum well interface is not perfectly smooth and ori- 
ented perpendicular to z. The energy difference between 
the two lowest eigenstates depends on the relationship 
between the phase of the fast oscillations of the wave 
function with the heterostructure boundary, and a step 
in the interface alters this phase relationship. The low- 
est energy wave function minimizes the energy, and, as 
shown in Fig. 8, can cause the phase of the fast oscilla- 
tions to become dependent on the transverse coordinates 
x and y. This coupling between the z-behavior and the 
x-y behavior is called valley-orbit coupling. 

As discussed in subsection |HI.B.1| above, in a silicon 
quantum well under tensile strain, there are two low- 
lying conduction band valleys at wavevectors +kgz and 
—koz, whose energies are split by the effects of confine- 
ment potentials and electric fields perpendicular to z. 
In the limit of a perfectly smooth interface aligned per- 
pendicular to z, the valley splitting of a quantum well 
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(c) E v = AE: Mixing of valleys and orbits 
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FIG. 9 (Color online) Valley-orbit mixing. (a,b) If the 
valley splitting Ev and orbital level spacing AE have very 
different values, the orbital and valley quantum numbers are 
well-defined and there will be no mixing of orbital and valley- 
like behavior, (c) When Ev~AE the valleys and orbits can 
hybridize in single-particle levels separated by the valley-orbit 
splitting E V o- 




FIG. 10 Valley-orbit coupling from interface steps. 

Top: gray-scale visualization of wave function oscillations in 
the presence of a perfectly smooth interface, oriented perpen- 
dicular to z. Middle: The relationship between the phase of 
the wave function oscillations and the interface is different on 
the two sides of an interface step. When the steps are close 
together, the phase does not adjust to the individual steps, 
and the valley splitting is suppressed. Bottom: When steps 
are far enough apart, the oscillations line up with the interface 
location on both sides of the steps, which causes the phase of 
the oscillations to depend on the transverse coordinate. This 
coupling between the behavior of the wave function in the z 
direction and in the x — y plane, which arises even when the 
well is atomically thin, is known as valley-orbit coupling. 



with typical width and doping is of order 0.1 meV, a 
magnitude that can be understood using the simple one- 
dimensional model presented in section I of the supple- 
mental material. 

If the step density of the quantum well interface is 
reasonably high, then the transverse oscillations of the 
charge density cannot align with the entire interface, and 



tion straddles a step, then the valley splitting is small, 
which, together with the effects of electron-electron inter- 
actions, causes the triplet state to have lower energy than 
the singlet state. If an electron is confined to a region 
small enough that it does not extend over multiple steps, 
then the valley splitting is not affected by the steps. Over 
the past several years, measurements of valley splitting 



valley splitting is greatly suppressed ( Ando 1979 Friesen in quantum point contacts ( Goswami et al. 2007 ) and of 



et al. 2007b 2006). The physical picture that emerges 
from effective mass theory that incorporates valley-orbit 
coupling is that the envelope function for the wave func- 
tion in a silicon heterostructure is qualitatively similar to 
typical wave functions in quantum dots, but that there 
are also fast oscillations with wave vector ~ fen in the 
z-direction. The fast oscillations of the two valley states 
have different phases. In the presence of interfacial disor- 
der such as interfacial steps, the value of the valley phase 
that minimizes the energy becomes position-dependent, 
so that one fixed value of the phase cannot minimize the 
energy everywhere, and the energy difference between the 
two different valley states decreases. This suppression 



explains measurements performed in Hall bars (Khrapai 



et al. 


2003 


Koester et al. . 1997 Lai et al. 2004| Weitz 


et al. 


1996 


) that yield very small values for the valley 



splitting of only /icV, and also why singlet-triplet split- 
tings in dots with two electrons have been observed with 
both positive and negative values at non-zero magnetic 
field (Borselli et al. 2011a) — if the electron wave func- 



singlet-triplet splittings in quantum dots ( Borselli et al. 
2011a|b[ [Simmons et qX||2011| [Thalakulam et aZ.||2011| ) 



in Si/SiGe heterostructures demonstrate that these split- 
tings can be relatively large, of order 1 meV, when the 
electrons are highly confined. These splittings are large 
enough that valley excitations are frozen out at the rele- 
vant temperatures for quantum devices (~ 100 mK). 

There are two different manifestations of valley-orbit 
coupling. The first, illustrated in the bottom panel of 
Fig. 10, occurs when the phase of the valley oscillations 
depends on the transverse coordinate. The second type of 
valley-orbit coupling can be visualized by considering an 
interface with a nonuniform step density. A wave func- 
tion localized in a region with few steps has larger valley 
splitting and hence lower energy than a wave function 



localized in a region with many steps (Shi et al. 2011). 
Therefore, the presence of the valley degree of freedom 
leads to translation of the wave function in the x-y plane. 
Valley-orbit coupling is important when the scale of the 
variations of the orbital and valley contributions to the 
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FIG. 11 A silicon-based nuclear spin quantum com- 
puter (a) Schematic of Kane's proposal for a scalable quan- 
tum computer in silicon using a linear array of 31 P donors in 
a silicon host. J-gates and A-gates control respectively the 
exchange interaction J and the wave function, as shown in 
(b). Reproduced from Kane (19981. 



energy are similar, a situation that occurs frequently in 
few-electron quantum dot devices. 

Because valley-orbit coupling and valley splitting de- 
pend on interface details, the observation of valley split- 



tings that vary substantially between devices (Borselli 



et al. 2011b I is not unexpected. Understanding and con- 
trolling this variability is important for being able to scale 
up the technology and for the development of devices that 



exploit the valley degree of freedom ( Culcer et al. 2009a 
|2012| |Li et al] |2010| |Shi et al] |2012|). Therefore, im- 



proved understanding of the physical mechanisms that 
affect valley splitting in real devices remains an impor- 
tant topic of active research. The valley-orbit coupling 
also contains phase information, which can be used for 



quantum computation (Wu and Culcer 2012) 



C. Dopants in Si 

1. Wave function engineering of single dopant electron states 

The central theme of quantum electronics applica- 
tions using single dopants is the ability to modify the 
dopant electron wave function using external electric 
fields and/or to manipulate the spin degrees of freedom 
using magnetic fields. In many proposals for dopant 
based qubits using either electron or nuclear spins as the 
qubit states, dopant electron wave function engineering 
is critical to effect single and two qubit gates. Since most 
work has been done on n-type dopants, this section will 
focus on donors. The original idea comes from the Kane 
proposal for a nuclear-spin based quantum computer in 



silicon ( Kane , 2000 ) where the single qubit operations are 
implemented by tuning the contact hyperfine interaction 
to bring the donor electron into resonance with a trans- 



verse oscillating magnetic driving field (see Fig. 11). To 
see this we write the effective spin qubit Hamiltonian of 
a single donor nucleus-electron system in the presence of 
a gate potential with strength V at the donor position as 
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FIG. 12 (Color online) Relative Stark shift of the con- 
tact hyperfine interaction for different donor depths 
(z) calculated for a uniform field in the z direction. 

(a) Using the tight-binding approach dMartins et aZ. [ |2004 1 , 



Using the tight-binding approach (Martins et al. 
and (b) Direct diagonalization in momentum space ([ 



Wellardl 



sonable 
ionization at 



and Hollenberg 20051. Agreement in overall trends is rea- 
and for the z — 10.86 nm case both methods predict 



-6 MV/m. 



Goan 


2005 


Kane 


1998) 



Hiq = [i B B z al - g n ^nB z a z n + A(V A )a n .a e , (8) 

where fj,B is the Bohr magneton, g n the Lande factor for 
P, and fi n is the nuclear magneton. The contact hyper- 
fine interaction strength A can be tuned by an applied 
electric field arising from a bias Va on an A-gate as: 



(9) 



where [1q is the permeability of silicon and ^(O, Va) is 
the donor electron wave- function evaluated at the nucleus 
under the A-gate bias V A - 

The change in the strength of the contact hyperfine 
coupling due to the application of a gate bias has been 
studied by several authors since Kane's proposal. To 
determine the change in the contact hyperfine coupling 
strength it is necessary to calculate the shift in the donor 
electron wave function at the position of the donor nu- 
cleus. Depending on the applied bias polarity, an A-gate 
control electrode will either draw the wave function to- 
ward, or away, from the gate. In either scenario the wave 
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function at the donor nuclear position is perturbed to 
some extent. The resulting tuning of A depends crit- 
ically on device parameters such as the depth of the 
donor from the interface, and the gate/interface geom- 
etry. The level of sophistication of the treatment of the 
donor electron wave function in these devices has steadily 
improved since the original calculations following |Kane| 
(1998). The earliest approaches used fairly simple hy- 



of the long range part of the donor potential ( Gerstmann 



2011 Overhof and Gerstmann 2004) 



In more recent years, the effect of depth and proximity 



to the interface on donor orbital states ( Calderon et al. 



2008 2006b Hao et al. 2009 Rahman et al. 2009a) 



drogenic wave functions scaled by the dielectric constant 
of silicon 



has received more attention as key experimental mea- 
surements became available. A turning point was the 
measurement of donor orbital states through transport in 
FinFET devices. The observed donor energy levels were 



Larionov et al. (2000) treated the bias poten- very different from the bulk spectrum (see section V.C) 



tial analytically, and the shift in the hyperfine interaction 
constant as a function of applied bias voltage was calcu- 
lated using perturbation theory. Wellard et al. (2002), 



again using scaled hydrogenic orbitals treated the prob- 
lem using a more realistic gate potential (modeled using a 
commercial semi-conductor software package, with built 
in Poisson solver). The donor electron wave function was 
expanded in a basis of hydrogenic orbitals in which the 
Hamiltonian was diagonalized numerically. Ke ttle et al.\ 
( |2003[ ) extended these calculations using a basis of non- 
isotropic scaled hydrogenic orbital states. |Smit et ~aL\ 



Extensive tight-binding calculations were used to explore 
the space of electric field and donor depth on the quan- 
tum confinement conditions of the donor-associated elec- 
tron, identifying Coulombic, interfacial, and hybridized 
confinement regimes. These calculations provided an 
excellent description of the low lying donor states ob- 



served and determination of the donor species (Lansber- 



(2003 2004) used group theory over the valley manifold 



and perturbation theory to describe the Stark shift of the 



gen et al. 2008). It would appear that the theoretical 
description of electric field "wave function engineering" 
of the donor electron across device dimensions is now 
well understood. The context of the Kane donor qubit 
has spurred further refinements of the theoretical descrip- 
tion of donor states, including the site-specific contact 



donor electron while Martins et al. (2005 2004) applied 



tight-binding theory to obtain the first description of the 
Stark shift of orbital states and the hyperfine interaction 
incorporating Bloch structure. Meanwhile, the effective 
mass treatment was further developed in a combined vari- 
ational approach ( Friesen 2005 1 and ( |Calder6n et al. 
2009), and in (Debernardi et al. 2006) using a Gaussian 



and non-isotropic hyperfine interaction terms (Ivey and 
Mieher[ |i"9~75a|bD for wave function mapping under elec- 
tric fields ( Park et al. 2009 ) , interaction with magnetic 



fields and gate control of the g- factor (Rahman et al. 



expansion of the effective-mass theory (EMT, see section 
II of the supplemental material) envelope functions. This 
was followed by the application of direct diagonalization 



2009b Thildcrkvist et al. 1994), dynamics of molecular 



donor-based systems (Hollenbcrg et al. 2004 Hu et al. 



2005| |Rahman et al.\ |2011b| [Wellard et aL[|2006|), c ross- 
talk in hyperfine control (Kandasamy et al. 2006), co- 



herent single electron transport through chains of ion- 



ized donor chains (Rahman et al 



in momentum space (Wellard and Hollenberg 2005) al- readout mechanisms (Fang et al. 



lowing the potential due to the A-gate to be included 
at the Hamiltonian level and gave a similar picture of 
the Stark shift of the hyperfine interaction as a function 
of external field strength and donor depth as the ear- 



2009b), spin-to-charge 



2002 Hollenberg et al. 



2004), and the calculation of donor levels in the pres- 



ence of STM-fabricated nanostructures providing modi- 
fications to the overall potential in a single-atom transis- 



tor, as shown in section V.B.3 (Fuechsle et al. 2012). 



lier tight-binding treatment of Martins et al. ( 2004 ) (see 
Fig. 12). Although not optimized computationally, the 
momentum space diagonalization approach has served as 
a consistency check against larger scale real-space tight- 
binding calculations of the Stark shift of the donor hy- 



perfine interaction at low fields (Rahman et al. 2007) in 



the overall benchmarking against experiment (Bradbury 



et al. 2006 1 which shows the theoretical description has 



converged to a reasonable level in terms of internal con- 
sistency and comparison with experiment (see Fig. 13). It 
should be noted that in such descriptions encompassing 
the overall donor electron wave function it is the relative 
change in the contact hyperfine interaction as a function 
of electric field that is computed since these approaches 
do not describe well the details of the electron state at 
the nucleus. Absolute calculations of the contact hyper- 
fine interaction are the domain of ab-initio theories where 
they have had remarkable success despite the truncation 



2. Two-donor systems and exchange coupling 

In the quantum computing context, the two main ap- 
proaches to directly couple the spins of donor electrons 
are through the Coulomb-based exchange interaction be- 
tween proximate donor electrons, or the magnetic dipole 
interaction. The Kane model uses gate control of the 
exchange interaction as per the two-qubit effective spin 
Hamiltonian: 

H 2 q = HBB z <r z ei - g^nB^a^ + Ai(V"Ai)<7 ril .<7 ei (10) 
+ /i B B 2 crf 2 - g n ^ n B z a^ 2 + A 2 (V A 2)vn 2 ■<Je 2 

+ J(Vj)ff ei .<7 e2 . 

In this equation we apply equation [8] on two dopants 
and add the exchange-coupling J between the dopants. 
There have been a number of papers investigating the 
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FIG. 13 Low-field Stark shift of the hyperfine inter- 
action for momentum space diagonalization (BMB) 
and tight-binding (TB) methods, (a) Electric field re- 
sponse of hyperfine coupling at various donor depths (BMB 
and TB). (b) Quadratic (lefthand axis) and linear (right-hand 
axis) Stark coefficients as a function of donor depth (TB) . (c) 
Shift of the ground state electron distribution (dipole mo- 
ment) as a function of the electric field (TB). (d) The electric 
field gradient of t he dipole moments as a function of donor 
depth (TB). From|Rahman et al.\$2007\. 
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FIG. 14 (Color online) J-oscillations in the exchange 

coupling. Calculated exchange coupling between two phos- 
phorus donors in Si (solid lines) and Ge (dashed lines) along 
high-symmetry directions for the diamond structure. Values 
appropriate for impurities at substitutional sites are given by 
the circles (Si) and diamonds (Ge). Off-lattice displacements 
by 10% of the nearest-neighbor distance lead to the perturbed 
values indicated by the sq uares ( Si) and crosses (Ge). Repro- 
duced fromlKoiller et al. (2002a I. 



construction and fidelity of two-qubit gates (e.g. such as 
the controlled-NOT) from this Hamiltonian (|Fang et al. 



ridge et al. 2006 Tsai et al. 2009 Tsai and Goan, 2008). 



2005] |Fowler ~eTal\ [20031 |Hill and Goan[ |2003t |2004| |Ker-| 



From a microscopic physics viewpoint, in general the ex- 
change energy J is stronger than the dipole interaction 



for smaller separations, however it behaves as (Herring 



and Flicker 1964) 



J(R) ~ (i?/a*) 5/2 exp(-2i?/a*), 



(11) 



where R is the donor separation and a* is the effective 
Bohr radius of the electron wave function. The exchange 
coupling dominates over dipole coupling for donors that 
are separated by less than approximately 20-30 nm. 

The valley degeneracy of the silicon conduction band 
gives rise to a far more complicated dependance of J on 
the donor separation (so-called "exchange oscillations") 



as noted in the early work of Cullis and Marko (Cullis 



and Marko 1970), and is particularly relevant in the 



Kane quantum computer context (Koiller and Hu 2005 



Koiller et al] |2002a| |2003[ ) (see Fig. 14). The effect 
persisted in effective mass treatments in which the ex- 
change integrals over Bloch states were carried out nu- 



merically (Koiller et al. 2004 Wellard et al. 2003). For 



some time these "exchange oscillations" were seen as a 
fundamental limitation of donor based quantum comput- 
ing as it was thought that to achieve a given exchange 



coupling the donors would have to be placed in the lat- 



tice with lattice site precision (Koiller et al. 2002a), al 



though Koiller et al. (2002b|) found that strain could be 



used to lift the valley degeneracy and alleviate the prob- 
lem to some extent. In these treatments the exchange 
coupling is calculated in the Heitlcr-London approxima- 



tion ( |Calder6n et al. 2006a Koiller et al. 2004) using 



effective mass wave functions containing a single Bloch 
component from each valley minimum, hence it is per- 
haps not surprising that the overlap integral results in 
an oscillatory behavior in the donor separation at the 
level of the lattice constant. When the exchange inte- 
gral is computed using a more accurate wave function 
including many such Bloch states to reproduce the ob- 
served donor levels and valley splittings, the interference 
effect is somewhat smeared out ( |Wellard and Hollenberg| 
2005 ) over the background Herring-Flicker dependence in 



equation 11 (see Fig. 15). Nonetheless, the issue remains 
that in fabricating donor devices there will be some level 
of imprecision in the donor atom placement and hence 
a variation in the (un-gated) value of J between donor 
pairs, however, using STM fabrication these errors might 
be constrained to the single lattice site level. 

In any case, all components of a quantum computer 
will need some form of characterization. For all donor 
qubit logic gates (single and two qubit), considerations 
of background noise sources and decoherence also need to 



be taken into account, e.g. see Fowler et al. (2003); Hill 
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FIG. 15 (Color online) Smoothing out the exchange 
oscillations - the exchange coupling J as a function 
of donor separation along [110]. Top curve: Calcula- 
tion using the effective-mass wave function. Middle curve: 
Calculation of J based on wave functions obtained using di- 
rect momentum diagonalization over a large basis of Bloch 
states (BMB) with no core-correction of the impurity poten- 
tial (77 = 0). Bottom curve: BMB calculation of J with a core- 
correction (r) — 5.8) that reproduces the donor ground-state 
and valley-splitting. Note that the points refer to substitu- 
tional sites in the silicon matrix. Although the donor separa- 
tions are relatively small in this case, the spatial variation of 
the exchange interaction appears to be significantly damped 
compared to the effective mass treatment. All J values are 
calculated in the Heitler-London approximation. Reproduced 
from |Wellard and Hollenberg| (|2005|) . 



and Goan 


(2003 


); 


Saikin and Fedichkin 


( 


2003 


); 


Wellard 


and Hollenberg 


(2001 2002) (the decoherence of donor 



electron spins is covered in Section VI). Robust control 
techniques have been developed specifically for the even- 
tuality of some level of variation in the exchange coupling 



(Hill 2007| , which in conjunction with gate characteriza- 



tion protocols ( Cole et al. 2006 Devitt et al. 2006 ) have 



the potential to produce high fidelity two qubit gates in 



the Kane scheme (Testolin et al. 2005). Tsai et al. (2009) 



have applied control techniques to optimize the CNOT 
gate in the Kane scheme. A more serious impediment to 
employing the exchange interaction for quantum gates 




FIG. 16 Gate control of the two-donor system. Aver- 
aged charge distribution along the interdonor axis for various 
strengths of the J-gate potential (fi) for the singlet (a) and 
triplet (b) states (fixed donor separation at 10 as). Repro- 
duced from|Fang etal\ p002|). 



out ( 


Calderon et al. 


2007 


Fang et al.\ 
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Kettle et al. 


2006 




2004 


Wellard and Hollenberg 


2004) 


given further 



insight into the controllability of the exchange interac- 
tion. However, the gate modification of the overlap be- 
tween electron states is a difficult calculation and most 
likely a full configuration interaction framework incorpo- 
rating valley physics and Bloch structure is required to 
obtain quantitative results to compare with experiments 
once measurements are made. A related problem is the 
calculation of the two-electron donor state (D~), notori- 
ously difficult in the case of a hydrogen ion in vacuum, 
but even more so when the non-trivial valley physics is 
added in to complicate such simple points of reference 
as Hund's Rule. In the context of single donor quantum 



computing Fang et al. (2002) calculated the effect of elec- 



tric fields on the D state, which was a key component of 
the spin-to-charge conversion read-out scheme of Kane. 



In Hollenberg et al. ( 2004 ) time-dependent calculations of 
the D u D u ->• D+ D transition were undertaken in a pro- 
posal for resonant based spin-to-charge conversion. More 
recent calculations have focussed on the complication of 
valley physics in the D~ bound states particularly un- 



der electric fields (Calderon et al. 2010a Rahman et al. 



2011a), with some notable success in comparisons with 



is the effect of charge noise ( |Hu and Das Sarma| |2006| iLansbergen et all 120081). 



recent experimental measurements (Fuechsle et al. 2012 



Vorojtsov et al. 2004 1. Because the exchange interac- 
tion is ultimately derived from an overlap of electronic 
wave functions, variations in the background potential 
such as from charge noise in the device can affect the 
exchange coupling and may require further development 



of the materials design (Kane 2005), and/or quantum 
control techniques. 

The control of the exchange interaction J has also 
received considerable attention since the original Kane 
paper. Early calculations of the dependence of J on 
an external J-gate bias were carried out by |Fang et al.\ 
(2002) using a Gaussian expansion (see Fig. 16). Sub- 



sequent calculations of the J-gate control in various ap- 
proaches describing the two-electron physics were carried 



3. Planar donor structures: delta-doped layers and nanowires 

The atom-by-atom fabrication of monolayer donor 
structures using STM techniques represents the state- 
of-the-art in precision silicon devices (see section V.B.3). 
From a theoretical point of view these structures present 
new challenges in order to describe not just their inher- 
ent physics (band structure, Fermi level, electronic ex- 
tent, valley splitting, effect of disorder etc), but their use 
as in-plane gates in quantum electronic devices, includ- 
ing quantum computing. In understanding the physics 
of these highly doped monolayer systems ab-initio tech- 
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niques have been used to good effect. Paradoxically, ab- 
initio techniques whilst being severely limited to rela- 
tively small numbers of atoms can handle planar sys- 
tems with a high degree of symmetry, exploiting peri- 
odic boundary conditions of the supercell in the plane 
of the structure with sufficient silicon "cladding" verti- 
cally for convergence. The earliest calculations in this 



context were by Qian et al. (2005) for the infinite 2D 



planar ("delta-doped") ordered layer using a Wannier- 
based Density Functional Theory (DFT) approach (see 
Fig. 17(a)). |Carter et al. (2009) carried out an exten- 
sive DFT calculation of the same Si:P structures us- 
ing a single zeta polarized basis providing a comprehen- 
sive picture of the band structure, effective potential, 
Fermi energy and electronic width as a function of pla- 
nar doping density, finding converged results for cladding 
above 80 layers (see Fig. 17(b)). More recently the ef- 
fect of disorder on the physics of the delta-doped layer 



et al. 2011 



has been investigated both in a DFT approach (Carter 
2009[ ), and in a self-consistent tight-binding 



approach which can handle much larger supercell sizes 
and hence more accurately represent instances of disor- 
der (Lee et al. 2011). These calculations indicate that 



the valley spitting of the sub-Fermi bands is quite sensi- 
tive to the degree of disorder and will play an important 
role in eventual device applications. 

The question of convergence between methodologies 
still remains on important quantities such as valley split- 
ting. Drumm et al. (2012a I have applied distinct DFT 



approaches based on localized and de-localized basis sets 
to calculate the properties of delta-doped layers. They 
obtain convergence in the description of the valley split- 
ting and Fermi level only when the localized basis set is 
extended to the double zeta polarized level. The DFT 
calculations of the band structure have informed a self- 
consistent effective mass description of Si:P monolayer 



structures (Drumm et al. 2012b), which has been effec- 



tive in describing states observed in a STM fabricated 



7-donor planar quantum dot (Fuechsle et al. 2010). The 
self-consistent tight-binding approach has also been em- 
ployed beyond the delta-doped layer to describe recent 
STM fabricated devices. In |Weber et al] ( |2012a[ ) the 
electronic structure of Si:P monolayer wires only four 
atoms wide was calculated and gave results in terms 
of the number of conduction modes in good agreement 
with experiment. The most ambitious calculation to date 



was a simulation of the single- atom transistor (Fuechsle 



et al. 2012 ) where the same self-consistent tight-binding 
approach was used to determine the effective potential 
due to planar gates at the channel-donor site and subse- 
quently coupled with a tight-binding description of the 
donor electronic levels. The agreement of the calculated 
D° and D _ charge transitions with the measurements was 
indeed remarkable given the complexity of the device and 
is a strong indication that the theoretical description of 
donor based quantum electronic devices is well in hand. 
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FIG. 17 Band structure of the 1/4 monolayer phos 
phorus (5-doped layer Top: (a) The calculation by Qian 



et al. (2005 I: the solid lines show the band structure without 
exchange-correlation and short-range effects, while the dotted 
lines show the band structure obtained in the full model, (b) 
T he DFT calcu l ation in a supercell with 200 cladding layers 
by Carter et al. (20091. The plane projected bulk band struc- 
ture of Si is represented by the gray continuum. The Fermi 
level is indicated by a horizontal dashed line. Reproduced 



from Qian et al. (20051 and Carter et al. (2009) 



IV. QUANTUM DOTS IN SI AND SIGE 

Quantum dots showing Coulomb blockade and display- 
ing single-electron physics can be created in Si and SiGe 
in many different ways. In this section we first briefly 
review the early work aimed at the demonstration of 
Coulomb charging effects in Si and SiGe. An empha- 
sis in this work was the quest to see Coulomb effects 
at as high a temperature as possible. We then discuss 
modern approaches to quantum dot fabrication. The ap- 
plication of charge sensing methods is shown to enable 
a wide range of experiments, including calibration of the 
absolute electron number, spin-state spectroscopy, and 
the measurement of spin filling as a function of electron 
number. We close this section with a discussion of both 
transport and charge sensing measurements in silicon- 
based double quantum dots. 



A. Early work: Coulomb blockade in silicon 

In this section we discuss early experiments studying 
Coulomb blockade in Si devices. Additional background 



and details can be found in 


Ahmed 


( 


1997); 


Likharev 


(1999); 


Meirav and Foxmanj ( 


199G 


); Ono et al. 


(2005); 


|Takahashi et al. 


(2002 


1; and Tilke et al. 


(2001 
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FIG. 18 (Color online) Calculated electronic spectrum 
of a single-atom transistor. Top left: Calculated energies 
of the D° and D~~ ground states (GS) as a function of the 
applied gate voltage Vg- The difference in the energy of these 
two ground states gives a charging energy of Ec ~ 46.5 meV, 
which is in excellent agreement with the measurement in this 
device. Potential profiles between source and drain electrodes 
calculated for Vg = 0.45 V (top middle) and 0.72 V (bottom 
left). The calculated orbital probability density of the ground 
state for the D° potential (top right) is more localized around 
the donor than for the D~ potential (bottom right), which is 
screened by the bound electron. Reproduced from |Fuechsle] 
et fflZ.M2012l. 



Experiments exploring intentional Coulomb blockade 
and transport through Si/Si02 and Si/SiGe quantum 
dots dates to the early 1990s, shortly after the discov- 



ery of Coulomb blockade (Field et al. 



1990 



Fulton and 



Dolan[ |1987| |Meirav et at] |1990| |Scott-Thomas et al. 
1989). The primary requirements for the observation 
of Coulomb blockade are to isolate a small island while 
maintaining a weak but nonzero tunnel coupling to the 
leads. The addition of one or more gates to control the 
charge on the dot is essential for more complicated ex- 
periments. 

Coulomb blockade was achieved very early in struc- 
tures formed by etching delta-doped SiGe or doped 



silicon-on-insulator (SOI) structures ( 


Ali and Ahmed 


1994 


Paul et al. 


1993). Ali and Ahmed ( 


1994 


) made use 



of two separate lithography and etching steps to mod- 
ulate the thickness of a patterned silicon-on-insulator 
layer, resulting in a weakly coupled island between two 
leads. Coulomb blockade, which was observed in mea- 
surements of current versus source-drain voltage that 
showed a Coulomb gap, persisted up to T = 3.8 K. The 
Coulomb gap could be modulated by an integrated side 
gate. In this type of highly-doped SOI structure, current 
in the doped leads was three-dimensional, as the mean 
free path was smaller than the lead thickness. 

Silicon nanowires formed in SOI can be transformed 



into a quantum dot by pattern-dependent oxidation (PA- 
DOX), a process that makes use of the dependence of 
oxidation in silicon on the exposed surface area and 
strain (Takahashi et al. 1994 1995). One of the fea- 



tures of this process is that very small quantum dots 
can be formed, enabling measurement of Coulomb os- 
cillations at high temperatures, with a demonstration of 
some modulation persisting to room temperature as early 



as 1994 (Takahashi et al. 1994 1. Fujiwara and co-workers 



studied the few-electron regime in similar devices using 



photoexcitation techniques (Fujiwara et al. 1997). Elec 



tron beam lithography can be used to help control the 
shape of small silicon dots that show Coulomb effects at 



temperature above 100 K (Leobandung et al. 1995). Very 



narrow triangular cross-section wires also can be formed 
by anisotropic etching on SIMOX, resulting in Coulomb 
effects at room temperature from disorder-induced dots 



along the length of each wire ( Ishikuro et al. 1996 ) . 

Coulomb blockade can in fact be observed in devices 
that are similar to production FETs, provided a small 
island of electrons can be isolated in the device. Isola- 
tion of such an island of electrons can be accomplished 
by the use of a gate that does not overlap the source 
and drain, leading to Coulomb blockade in CMOS de- 
vices ( Boeuf et al. 2003 ) . This approach has culminated 



very recently in a single-electron transistor operating at 



room temperature (Shin et al. 2010 2011a) 



In 1994 Matsuoka and co-workers proposed using "two- 



story gates" to create single-electron devices (Matsuoka 



et al. 1994 ) . A single gate was used to form an inversion 



layer for transport, and an upper gate was reverse-biased 



to generate barriers and define a quantum dot ( Matsuoka 



and Kimura 19951. While this structure has only a single 
gate to control the tunnel barriers and differs in signifi- 
cant ways from later work, it anticipates the use of two 
layers of gates that would be used more than a decade 
later for experiments on spin blockade, spin measure- 



ment, and spin manipulation (see Sections IV. F. 2 and 

vib. 



B. Single quantum dots 

This section assesses the experimental analogues of the 
quantum dot concepts different in silicon nanostructures 
as explained in section [II. A| 



1. Self-assembled nanocrystals 

The material dimensions of nanocrystals can easily be 
made as small as 10 nm, resulting in large and thus 
easily observable level splittings, even at room temper- 
ature (Otobe et al. 1998). On the other hand, those 
dimensions make electron transport measurements cum- 
bersome because the crystals are not easily connected 
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to source and drain reservoirs. Self-assembled silicon 
nanocrystals with diameters varying from 3-12 nm have 
been grown by chemical vapor deposition techniques 



(Baron et al. 2000 Steimle et al. 20071. Coulomb os- 



cillations have been observed by electrostatic trapping 



between Al source and drain electrodes (Dutta et al. 



2000). Zaknoon et al. (2008) showed charging energies of 



~ 50 meV using scanning tunneling spectroscopy. Twelve 
resonances in the conductance versus bias voltage were 
attributed to the twelve-fold conduction band degener- 
acy owing to spin and the six-fold valley degeneracy as 
described in section UlLAI 
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FIG. 19 (Color online) Self-assembled nanocrystals. (a), 
STM image of a Ge/Si(001) cluster with a height of 2.8 nm. 



Scan area is 40 x 40 nm, from (Mo et al. 1990). (b), Band di- 



agram for a Si/Ge/Si heterostructurc, showing the accumula- 
tion of holes owing to the valence band offset between Ge and 
Si. (c), Schematic of a quantum-dot device obtained by con- 
tacting a single SiGe nanocrystal to aluminum source/drain 
electrodes. The heavily doped substrate is used as a back- 
gate for the measurements in (d) where 7sd is plotted a s a 
function of Vb and Vsd- (c,d) from |Katsaros et aZ. | ( |20To| ) . 



Small Ge islands can be grown on Si(001) via Stranski- 
Krastanov growth resulting in huts, pyramids and domes 
with heights of 5-70 nm and lateral dimensions varying 



Due to the limited tunability reaching the few-charge 
regime in self-assembled nanocrystals will be a great 
challenge. 



2. Bottom-up grown nanowires 

Bottom-up grown nanowires are generally synthesized 



by means of a vapor-liquid-solid process (Wagner and 



Ellis 1964), allowing for growth of single-crystal Si and 



Ge nanowires (Morales and Lieber 1998) with diameters 



varying from 3-100 nm and lengths up to tens of microns, 
see Fig. [20[a,b). Both n-type and p-type dopants have 
been incorporated, and their location depends on the di- 
ameter ( Xie et al. 2009 ) . The doping can be varied dur- 



ing growth: such modulation doping has been used to in- 
tersect heavily-doped n-Si regions with two short lightly- 
doped regions, resulting in single-electron tunneling at 
1.5 K (Yang et al. 2005). Within one nanowire het- 
erostructures of different materials can be created both 
radially and axially, such as core/shell Ge/Si nanowires 
( Lauhon et al. 2002 ) . In the latter case the valence band 



offset will induce hole population in the Ge core, see 



Fig.[19Jb). 

When metallic contacts are made to nanowires the 
Schottky tunnel barriers can define the quantum dot 
length as shown in core/shell Ge/Si nanowires (Lu 

see 



et al. 



Fig. 20 



2005) 



20051 and Si nanowires (Zhong et al. 
c). The Si nanowire quantum dot length can 
be shortened by silicidation transforming the device into 
e.g. a NiSi-Si-NiSi nanowire as shown in Fig. [20^d) 



et al. 2009a I 



(Mongillo et al. 2011 Weber et al. 2006 Zwanenburg 



After the demonstration of Coulomb blockade oscilla- 
tions in Ge/Si nanowires by the Lieber group from Har- 
vard (Lu et al. 2005 ), they joined forces with the Marcus 



group and created double quantum dots with tuneable 
tunnel barriers, see Section IV. F Here the source and 
drain contacts were ohmic, while the tunnel barrier were 



defined by local top gates ( Hu et al. 2007 1 . Roddaro 



from 20-80 nm (Eaglesham and Cerullo 


1990 Kat- 


saros et al.\ 2008| Medeiros-Ribeiro et al. 


1998[ Mo 


et al. 1990 Ross et al. 1999| Stangl et al. 


2004), see 


Fig. 19 a) . The group of De Franceschi in Grenoble 



et al. ( 2008 ) used the same configuration to create single 
quantum dots and probe the hole spin states, see sec- 
tion IV. E Ge/Si nanowires were found to have a strong 
spin-orbit interaction, which can be tuned by means of an 
made Al contacts to Ge domes with an additional 2 electric field ( Hao et al. 20101. Recent spin lifetime mea- 



nm Si capping layer (Katsaros et al. 
Fig. [T9] In this configuration the SiGe nanocrystal acts 
as a confining potential for holes due to the valence band 
offset between Ge and Si at the heterostructure interface 
dSchafflerl 119971 [Van de Walk and Martini 119861). Free 



2011 2010), see surements ( Hu et al. 2011 ) indicate spin-orbit interaction 



as the dominant mechanism for spin relaxation. Accord- 



ing to the work by Klocffcl et al. (2011), the unusually 



strong spin-orbit coupling makes them particularly at- 
tractive candidates for quantum information processing 



holes will accumulate in the Ge when the Fermi level 
lies below the valence band edge of the Ge center, see 
Fig. [l9{h). Electron transport measurements at 15 mK 
show Coulomb diamonds with charging energies varying 
from few to 20 meV as 8 holes leave the quantum dot. 



via electric-dipole induce spin resonance (Golovach et al. 



2006 Nadj-Perge et al. 2010 Nowack et al. 2007), and 



for research on Majorana fermions (Majorana 1937). 

Very recently, Ge/Si nanowires with a triangular cross 
section and a height of just three unit cells were realized 
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FIG. 20 (Color online) Bottom-up grown nanowires. (a), 
TEM image of a Si nanowire; crystalline material (the Si core) 
appears darker than amorphous material (SiOa; sheath) in this 
imaging mode. Scale bar, 10 nm. (b) High-resolution TEM 
image of the crystalline Si core and amorphous SiO x sheath. 
The (111) planes (black arrows) are oriented perpendicular 
to the growth direction (white arrow), (a) and (b) adapted 
from Morales and Lieber ( |1998 1, (c) Stability diagram of 
a p-Si nanowire quantum dot, from Zhong et al. (20051 (d) 
SEM image of a nanowire quantum dot with NiSi Schottky 
contacts, taken from Zwanenburg et al. (2009a I. 



by molecular beam epitaxy (Zhang et al. 2012). These 



wires are directly grown on planar Si without the use of 
any catalyst, and preliminary low-temperature measure- 
ments show Coulomb blockade. 



3. Electrostatically Gated Si/SiGe quantum dots 

A powerful way to achieve tunability of tunnel cou- 
plings in quantum dots is to provide confinement in one 
or more directions through the use of electrostatic gates. 
Using Si/SiGe heterostructures or MOS structures, it is 
possible to form high-quality two-dimensional electron 
systems that can be partitioned into tunable quantum 
dots using depletion or accumulation gates, a procedure 
described in detail in this section and the next. In gen- 
eral, at least direction of confinement must be provided 
by a non-electrostatic method; usually a materials inter- 
face is used, the two most common being the interface 
between single-crystal silicon and its amorphous oxide 
(in MOS structures, see next section), and the epitaxial 
interface between single- crystal Si and Sii^Ge^. When 
the precise composition x is unimportant and no con- 
fusion will arise, we refer to these heterostructures as 
Si/SiGe. Both MOS devices and Si/SiGe devices have 
been reviewed extensively: see, for example, (Sze and 
~ Wolfl[T990| for the former, and ( |Mooney||1996 



Ng 1981| 



Schaffier 



A convenient, if incomplete, figure of merit for two- 
dimensional electron systems is the mobility fi. For Si 
MOS, mobilities in the range 5, 000 — 15, 000 cm 2 /Vs are 



quite good (see e.g. Eng et al. (2005 2007)), and mobil- 
ities in excess of 40,000 have been reported (Kravchenko 



and Sarachik 2004). The low-temperature mobility in 



Si/SiGe two-dimensional electron gases is not limited by 
defects at the interface and has been improving rapidly 
in recent years. In 1995, Ismail and coworkers reported a 
low-temperature mobility of 520, 000 cm 2 /Vs in a mod- 
ulation doped Si/SiGe heterostructure. Even higher mo- 
bility 800, 000 cm 2 /Vs was reported by a group from Hi- 
tachi in 1998 ( |Sugii et aZ.[|1998| ). Very recently, Si/SiGe 
two-dimensional electron systems have been formed using 
undoped structures with a positively-biased accumula- 
tion gate. In this approach, an intervening oxide such as 



AI2O3 (Lai et al. 2005 ) is used to separate the accumula- 



tion gate from the semiconductor surface to avoid inject- 



ing current into the heterostructure ( Lu et al. 2007). The 



positively biased accumulation gate removes the need for 
any doping in the structure, removing a source of back- 
ground impurities and eliminating the modulation doping 
layer altogether, both of which cause scattering. Result- 
ing mobilities as high as \i = 1.6 x 10 6 cm 2 /Vs have 
been reported (Lu et al. 2009). Further, the removal 



of intentional doping appears to significantly reduce low- 
frequency charge noise in the devices. 

Because both Si and Ge have isotopes with zero nuclear 
spin, the proposal by Loss and DiVincenzo to use quan- 



tum dots as hosts for semiconductor spin qubits (Loss 



and DiVincenzo 19981 led to great interest in the de- 



velopment of high-quality quantum dots in Si/SiGe het- 
erostructures (Friesen et al. 2003 Vrijen et al.\ 2000). 
The challenge in the early work in this field was to find 
ways to fabricate such dots with low-leakage gates, suffi- 
cient tunability, and in such a way as to yield stable, low- 
noise devices. As we discuss later in this review, modern 
Si/SiGe quantum dots have achieved performance that 
rivals that of any materials system available. In this sec- 
tion we discuss the materials and device research that 
enabled this advance. 

Here we discuss a few critical materials issues relevant 
to Si/SiGe heterostructures. Interest in Si/SiGe arises 
because of the inevitability of defects at the interface 
between crystalline Si and its amorphous oxide. Het- 
erostructures formed from Si and Sii^Gc^ offer a nat- 
ural alternative with, in principle, no interfacial traps 
(although other types of disorder, such as atomic steps 
and strain variation are certainly present). 

Although both Ge and Si have the diamond structure, 
Ge sits one row beneath Si in Group IV of the periodic 
table, so that the lattice constant of Sii-zGez increases 
as x increases, achieving a mismatch between pure Si and 



Ge of approximately 4.17% (Schaffier 19971. Because of 



1997) for the latter. 



this mismatch, pure Ge will grow epitaxially only three 
monolayers on Si (REF). Beyond this critical thickness, 
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self-assembled quantum dots or "huts" form (Mo et al. 



1990), as discussed in Sec. IV.B.l preventing the growth 



of uniform quantum wells. 

Because the lattice constant of Sii^Ge^ depends on 
x, a full description of a heterostructure of these two ma- 
terials must include the strain of the various layers. For 
the structures considered here, the layers of interest typ- 
ically include a Si quantum well with Sii-^Ge^ barriers 
on either side, as shown in Fig. [21] typically, x ~ 0.3. 
If the quantum well is below the critical thickness for 
dislocation formation, the in-plane lattice constant will 
remain unchanged passing vertically from the Sii-^Ge^ 
through the Si quantum well and into the upper barrier. 
The band offsets at the Si/Sii-xGejc interfaces depend on 
this in-plane lattice constat. For an unstrained, relaxed 
Sio.7Geo.3 barrier layer, the minimum in the conduction 
band is approximately 160 meV lower inside a Si quan- 



tum well compared with the barriers (Schaffler 1997) 




FIG. 21 Layer design and corresponding band dia- 
gram of a Si/SiGe modulation doped heterostructure 
used to form top-gated quantum dots. Reproduced from 
Berer et al.\ \2QQ7\ . 



Because it is very challenging to grow bulk, re- 
laxed Sii-zGe^ with even moderately large x, relaxed 
Sii-^Gea; substrates conventionally are formed by slowly 
increasing the Ge concentration x from zero to the de- 
sired final value over a thickness of several microns. This 
procedure induces the formation of misfit dislocations, 
increasing the overall lattice constant, and can yield low- 
defect structures (Mooney 1996). The relaxation pro- 



cess itself does result in small inhomogeneities, which can 



be observed with nano-beam x-ray measurements ( Evans 



et al. 2012) 



Quantum dots in Si/SiGe demonstrating Coulomb 
blockade were first formed using a combination of etching 



and electrostatic gating. |Notargiacomo et al. (2003) ob 



etched into a Si/SiGe heterostructure. This early de- 
vice had a single overall top gate used to control the 
number of electrons in the quantum dot. |Klein et al.\ 
(2004) formed a quantum dot with three separate elec- 
trostatic gates. These gates were formed of the same 
two-dimensional electron gas as the quantum dot, source 



and drain leads ( Eriksson et al. 2004 ) . To avoid current 



flowing from the gates to the dot, deep trenches were 
etched between the gates and the dots; the intervening 
gaps make it difficult to apply local fields and separately 
gate the quantum dot and the tunnel barriers. This draw- 
back was partially ameliorated by the demonstration that 
gates could be formed by metal deposited into etched re- 



gions surrounding the dot ( Sakr et al. 2005 ) , and by the 



use of extremely small top gates used to break an etched 



wire into a gated quantum dot ( Slinker et al. 2005 ) . The 
drawback of etching, however, is the potentially large de- 
gree of side- wall depletion ( Klein et al. 20061) . 
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FIG. 22 (Color online) (a) Scanning electron micrograph 
of the Schottky gates used to form a gated quantum dot in 
Si/SiGe. (b) Coulomb diamonds: conductance of the dot as a 
function of the voltage Vb applied to gates Gl and G2 and of 
the drain-source voltage Vds- Figure from Berer et al. (20061. 



Berer et al. (2006) demonstrated a fully top-gate de- 
fined quantum dot formed in a modulation-doped 2DEG, 
as shown in Fig. 22 They showed that Pd Schot- 
tky gates, when fabricated on heterostructures like that 
shown in Fig. |21| in which care was taken to reduce 
the dopant density near the surface, enabled low-leakage 



gates (Berer et al. 2006). There had been great concern 
about leakage between the top gates and the electron- 
gas, but the Pd Schottky gate approach has proven to 



be very robust (Klein et al. 2007 Payette et al, 2012 



Wild et al. 2010). The Schottky gate approach has 



also been used successfully to gate heterostructures with 



enhanced concentration of 28 Si and 70 Ge (Sailer et al. 



2009 ) . A second approach to eliminating leakage is to use 
a dielectric material beneath the gates, creating metal- 
oxide-semiconductor split gates to define the quantum 
dot (IShin et aZ.l |2011b|). 



served Coulomb blockade oscillations in a gated nanowire 



The primary advantage of top-gated quantum dots, in 
which the lateral confinement is entirely provided by ad- 
justable gate voltages, is their extreme tunability. At 
zero gate voltage in most cases current can flow directly 
under a gate, enabling a smooth transition from a com- 
pletely open two-dimensional electron gas to a fully con- 
fined quantum dot. This tunability led both to the obser- 
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vation of the Kondo effect in a Si/SiGe top-gated quan- 



tum dot (Klein et al. 2007) and the demonstration of 



single-electron occupation, as shown in Fig. [28] below. 



4. Quantum dots in planar MOS structures 

The silicon MOSFET is arguably the world's most im- 
portant electronic device, being the basic component of 
all modern microprocessor chips. Its success has been 
built on the ability to grow a high-quality SiC>2 layer on 
the Si(001) surface by thermal oxidation, forming a high 
band-gap insulator that isolates the gate from the sil- 
icon channel. In current processor chips a SiC>2 layer 
of ~1 nm is sufficient to maintain gate voltages that 
are a significant fraction of a volt with negligible leak- 
age. The Si/SiC>2 interface, which confines the electron 
layer in a MOSFET, can also have relatively low disor- 
der, with reported electron mobilit i es as high as 40,000 
cm 2 /Vs (Kravchenko and Sarachik 2004), although the 



imperfect lattice match between the Si and SiC>2 creates 
defects at the interface, thus constraining the electron 
mobilities below those attainable at Si/SiGe interfaces. 
Despite this, it is possible to form quantum dots in MOS 
structures that can be controlled down to the single elec- 
tron level with high tunability. 

In this section we focus on quantum dots formed at 
the Si/Si02 interface via the use of multiple surface gates 
that provide electrostatic confinement in all three dimen- 
sions. In general an upper gate is used to induce an 
electron layer at the interface (as in a 'traditional' MOS- 
FET), while two or more lower gates provide tunable 
tunnel barriers between the electron reservoirs and the 
dot. As already described in Section IV. A, one of the 



earliest such structures (Matsuoka et al. 1994) exhibited 



Coulomb blockade oscillations, although these prelimi- 
nary results were rather irregular. 

One of the first well controlled MOS quantum dots was 
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FIG. 23 Si MOS quantum dot with large-area top 
gate, (a) Cross-sectional schematic, showing two oxide and 
two gate layers, formed on a silicon substrate. The lower 
Si02 layer is thermally grown, while the upper oxide layer is 
formed using plasma deposition. The large-area upper gate 
induces a 2DEG at the Si/Si02 interface, while the lower 
gates locally deplete the 2DEG to form a quantum dot. (b) 
Top-view schematic, showing lower depletion gates (black) 
and induced electron layer (grey), (c) Normalized spacings 
8 between Coulomb peaks in dot conductance as a function 
of upper gate voltage. Inset: Raw Coulomb oscillations in 
dot conductance as a function of upper gate voltage. Data 
reproduced from Simmel et al. (1999 1. 



Si) is used for the lower gates, and a large area upper 



metal gate is used to induce the 2DEG layer (Nordberg 



et al. 


2009a. 


Tracy et al. 


2010) 



is appealing from the perspective of future manufactur- 
ing, since it opens the way towards the use of CMOS 
process technologies. Similar MOS quantum dots also 
have been used to confine single electrons, enabling di- 



rect measurement of electron spin relaxation times ( Xiao 



et al. 2010a) 



By reducing the upper MOSFET gate to nano-scale di- 
mensions, a group at the University of New South Wales 
developed a highly compact multi-gate MOS architecture 



demonstrated by Simmel and co-workers (Simmel et al. ble (Lai et al. 2011 Lim et al. 



(Angus et al. 2007 ) that has since been used to construct 
a wide range of single (|Lim et al. |2011| [2009b]) and dou 



1999), see Fig. 23 In this structure a continuous upper 
gate was used to induce a 2DEG over a large area, while 
four lower gates were used to confine the dot and form 
tunnel barriers. The resulting lower gate structure mim- 
ics those used to confine GaAs/AlGaAs quantum dots, 
although in the latter case the 2DEG is created by modu- 
lation doping. The resulting Coulomb oscillations in this 



2009a) quantum dot 



MOS device were quite regular (Fig. 23;) and provided 
promise for future MOS quantum dot studies. The lower 



gates of the device in Fig. 23 were made using refrac- 
tory metal, since a high-temperature process was used to 



deposit the upper oxide isolation layer (Fig. 23 1) . 

This type of architecture, employing a large-area up- 
per gate, has since been used by a number of groups to 
construct MOS quantum dots. A group at Sandia Na- 
tional Laboratory has demonstrated a range of quantum 
dot devices in which etched polycrystalline silicon (poly- 



structures. This architecture uses aluminum (Al) upper 
and lower gates, with a thin (3-5 nm) AI2O3 insulat- 
ing layer between the gates, formed by thermally oxidiz- 
ing the lower gates at the relatively low temperature of 
150 C . Despite being very thin, the AI2O3 insulator can 
maintain inter-gate voltage differentials of up to 4 volts, 
allowing for high gate tunability and the formation of 
very small (sub-50 nm) multi-dot structures. Figure 24 
shows a quantum dot device based on this technology, in 
which a third layer of gate metal is used. This allows one 
upper gate to be used as a 'plunger', to control the dot's 
electron occupancy, while separate upper gates are used 
to induce the source and drain electron reservoirs - see 



Fig. 24 'b) . In this way the dot occupancy can be reduced 
to the single electron level, as confirmed by the bias spec- 
troscopy measurements in Fig. |24[ c), while maintaining 
a high density of states in the reservoirs. Such inde- 
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FIG. 24 (Color online) Si MOS quantum dot with com- 
pact multi-layer gate stack, (a) Scanning electron micro- 
scope image of device, (b) Cross-sectional schematic, showing 
three oxide layers and three Al gate layers, formed on a sili- 
con substrate. The Si02 layer is thermally grown in a high- 
temperature process, while the thin AI2O3 layers between the 
gates are formed by low-temperature oxidation of the alu- 
minum, (c) Stability map obtained by plotting differential 
conductance through the device as a function of source-drain 
bias Vsd and plunger (P) gate voltage Vp. The first diamond 
opens up completely, indicating that the dot has been fully 
depleted of electrons, (d) Coulomb oscillations as a function 
of plunger gate voltage Vp for the first 23 electrons in the dot. 
Data reproduced from Lim et al. (2011 1. 



pendent tuning of the dot occupancy and the reservoir 
electron density is not possible when a large-area upper 
gate is employed. 

The metal-oxide-semiconductor techniques just dis- 
cussed can be applied to Si/SiGe heterostructures, yield- 



ing extremely stable and tunable quantum dots ( Borselli 



et al. 


2011a 


Hayes et al. 


2009) 



The device design, as 
shown in Fig. |25| uses a Si quantum well surrounded by 
epitaxial SiGe barriers to provide a clean environment for 
the electrons in the device. Those electrons are induced 
by an accumulation gate at the top of the structure. De- 
pletion gates in between the accumulation gate and the 
heterostructure surface are used to control size and shape 
of the dot. 

MOS-based quantum dots, using architectures like 
those in Figs. [23] and |24j have since been used in a 
range of advanced measurements, including single-spin 
measurement, and spin- and valley-state spectroscopy, as 
will be discussed in Sections IIV.EI and IIV.F1 
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FIG. 25 (Color online) Gated quantum dot formed from 
a Si/SiGe heterostructure with a global accumulation gate, 
(a) Cross-sectional view of the heterostructure and the two 
layers of gates, (b) Top-view SEM image of the gates with 
a numerical simulation of the electron density superimposed. 
Figure from iMaune et al.\ (12012 1 . 



5. Quantum dots in etched silicon nanowires 



As discussed in Section |IV.A[ some of the earli 



est silicon-based single-electron devices (e.g., Takahashi 



et al. ( 1994 1995 ) were based upon narrow nanowires, 
patterned using traditional top-down lithographic tech- 
niques, and etched from thin (typically < 50 nm) silicon 
layers that form the upper layer of silicon-on-insulator 
(SOI) wafers. These early devices used the pattern- 
dependent oxidation (PADOX) technique to create ad- 
ditional confinement along the length of the nanowire, 
but in subsequent structures researchers have incorpo- 
rated 'wrap-around' gates, positioned along the wire to 
provide additional confinement. 

One of the first examples of this type of gated silicon 
nanowire was demonstrated by a group at NTT in Japan 



(Fujiwara et al. 2006 ) - see Fig. 26 Here, confinement in 
the y and z directions was provided by the narrow wire, 
of width 20 nm and thickness 20 nm. Confinement along 
the wire was created by wrap-around lower gates, which 
in this case were made from poly-Si. Finally, a large-area 
poly-Si upper gate, isolated from the gates below using 
Si02, was patterned above the entire structure to induce 
carriers in the nominally un-doped nanowire. The re- 
sulting structure is entirely CMOS compatible, making 
it convenient for production using well established man- 
ufacturing processes, and also utilizing the high-quality 
thermally grown Si02 insulator, which is known for hav- 
ing very low charge noise. In subsequent measurements 
on these devices it was found that they exhibited ex- 
tremely high charge stability, with a drift of less than 



O.Ole over several days (Zimmerman et al. 2007) 



As seen in Fig. 26 'd) , a quantum dot could be formed 
by using the outer gates LGS and LCD to create tunnel 
barriers, with the central gate LGC acting as a 'plunger' 
to control the dot occupancy. The Coulomb oscillations 



(Fig. 261) were highly periodic over a large gate volt- 
0.5V < Vlgc < 1-0V), with a deviation of 
although the dot occupancy N e in 



age range ( 
less than 1 percent 
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FIG. 26 (Color online) Multi-gated quantum dot in 
etched silicon nanowire. (a) Schematic top- view and cross- 
sectional view of the device. Three lower 'wrap-around' gates 
(LGS, LGC, LGD) are used to form tunnel barriers in an 
etched silicon nanowire. (b) Top-view scanning electron mi- 
croscope image of the device before the upper gate is de- 
posited, (c) Equivalent circuit of the device, (d) Coulomb 
blockade oscillations in device conductance as a function of 
central gate voltage Vlgc when the two outer gates (LGS, 
LGD) are biased to set each tunnel barrier to G = 1 fiS. 
Inset: Coulomb oscillations for a range of values of barrier 
conductance from 20 nS to 8 fiS. Data reproduced from |Fujr] 
wara et al. ( 2006 1. 



this case was relatively large, with N e ~ 200 electrons 
at Vlgc — 0V. The peak conductance could also be 
tuned over more than three orders of magnitude by vary- 
ing the barrier gate voltages. For central gate voltages 
Vlgc < — 1-0V, an additional tunnel barrier was formed, 
breaking the quantum dot into two dots in series. Using 
similar device structures this group could therefore oper- 
ate double quantum dots, demonstrating effects such as 
Pauli spin blockade ( Liu et al. 2008b ) - discussed further 
in Section VI.C.4. 

It is also possible to form a quantum dot in a silicon 
nanowire using just a single gate, by making use of tech- 
nology that has been developed for the manufacture of 
FinFET-type MOSFETs. FinFETs are considered likely 
replacements for planar CMOS technology, due to their 
ability to operate as FETs with good ON/OFF ratios 
at much shorter channel lengths. Figure |27^a,b) shows 
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FIG. 27 (Color online) Single-gated quantum dot in 
etched silicon nanowire. (a) SEM images and (b) cross- 
sectional schematics taken perpendicular to the nanowire 
(upper) and along the nanowire (lower), (c) Stability map 
(Coulomb diamonds) obtained by plotting differential current 
through the device as a function of source-drain bias Vd and 
wrap-gate volta ge V g . Data re produced from |Sellier et al.\ 



([2006} and from |Hofheinz et al.\ (|2006bJ 



a FinFET structure, which is based upon a nanowire 
(the 'fin') that is etched from a SOI wafer, as previ- 
ously described. A single wrap-around poly-Si gate is 
encapsulated on either side by an insulating 'spacer', 
made from either Si02 or SiaN4. The gate and spacer 
act as a mask for subsequent ion implantation of the 
n+ source and drain regions, which is a standard 'self- 
aligned' gate process used in CMOS production. By ap- 
plying a positive voltage to the poly-Si gate electrons 
can be induced below, to form a quantum dot, isolated 
from the source and drain due to the natural barrier cre- 
ated by the spacer regions - see Fig. 27 Id) . Such quan- 
tum dots can be extremely stable in the many-electron 
regime, as shown in Fig. [27^c), which demonstrates bias 
spectroscopy ('Coulomb diamonds') taken over a wide 
range of electron occupancy, with high stability and al- 



2006b). 



most constant charging energy (Hofheinz et al 
Similar FinFET structures have also been used for single 
dopant tunneling studies - see Section [V.B.2| 



C. Charge sensing techniques 

The non-invasive sensing of charge displacements in 
quantum nanostructures was first demonstrated in a 



GaAs/AlGaAs heterostructure device (Field et al. . 1993 1, 
when a quantum point contact (QPC) was used to de- 
tect the change in occupancy of a quantum dot. Here, 
the QPC is biased close to pinch-off, where its transcon- 
ductance dl/dVc can be very large. Any small charge 
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displacement in the vicinity of the QPC channel can 
then lead to a significant change in QPC current, via 
its capacitive coupling. This technique has since been 
applied widely, enabling the direct probing of single elec- 
tron charges and the indirect probing of single spins in 
nanostructures based on a variety of materials systems, 
including silicon. 



sure the spin state of a single electron in a GaAs/AlGaAs 



Sakr et al. ( 2005 ) fabricated a QPC adjacent to a quan- 
tum dot in a Si/SiGe heterostructure using a combination 
of isolation etching and metal gates aligned to the etched 
trenches. While this structure enabled sensing of the 
dot's electron occupancy in the many-electron regime, 
it did not have sufficient sensitivity to probe down to 



the last electron. Simmons et al. (20071 used Pd metal 



surface depletion gates on a Si/SiGe heterostructure to 
define a similar geometry - see Fig. |28^ . By monitoring 
the differential conductance of the QPC sensor they were 
able to accurately probe the depopulation of electrons in 
the adjacent quantum dot, even when the transport cur- 
rent /not through the dot had fallen below the noise level 
(Fig. [28})) . In this way they were able to track the occu- 
pancy of the dot down to the final electron, as shown in 
Fig.ptc). 
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FIG. 28 (Color online) Non-invasive charge sensing of a 
Si/SiGe quantum dot using a quantum point contact 
(QPC) sensor, (a) SEM device image, (b) (Top) Derivative 
of the QPC current (IIqpc /dVa as a function of gate voltage 
Vg- The peaks correspond to changes in the number of elec- 
trons in the dot. (Bottom) Current loot through the quantum 
dot as a function of Vg- (c) QPC sensor output in the few- 
electron limit. No further transitions occur for Vg < 1-68 
V, indicating an empty quantum dot. From Simmons et al. 
(2W7\. 



More recently, QPC sensors have been used with great 
versatility in both Si/SiGe and Si MOS quantum dots 



systems for measurement of both charge ( 


Nordberg et al. 


2009b 


) and spin ( 


Hayes et al. 


2009 


Simmons et al. 


2011 


Xiao et al. 


2010a 


) states. A technique developed to mea- 



quantum dot (Elzerman et al. 2004) has been success- 



fully applied to dots in silicon. This involves loading an 
electron (of indeterminate spin) into an empty quantum 
dot and positioning the Fermi level so that only a spin- 
up electron is able to tunnel out, with the charge dis- 
placement monitored by a QPC sensor. The technique 
has been used to measure the spin lifetime of single elec- 



trons loaded into Si/SiGe (Hayes et al. 2009 Simmons 



et al. 2011 land Si MOS (Xiao et al. 2010a) quantum 



dots. These experiments are discussed in more detail in 
Section ED 

Single electron transistors (SETs) can also been used as 
highly sensitive electrometers in nanostructure devices. 
The most sensitive such electrometers employ Al metal 
islands, with AI2O3 tunnel barriers, which can be in- 
tegrated with both MOS ( |Andresen et al.\ |2007[ ) and 
Si/SiGe-based quantum dots ( Yuan et al.\ 2011). Inte- 
grating such SETs into a radio-frequency (rf ) tank circuit 



forms an rf-SET (Schoelkopf et al. 1998), which can op- 



erate at frequencies above 100 M Hz with charge sensitiv- 
ities approaching ^10~ 6 e/\/Hz. Andresen et al. (2007) 
fabricated such an AI-AI2O3 rf-SET on the surface of a 
phosphorus-doped silicon (Si:P) device to study the gate- 
controlled transfer of an electron between two implanted 
phosphorus donors, with a measurement bandwidth ex- 
ceeding 1 MHz. They were able to study the charge relax- 
ation rate as a function of gate- induced detuning between 
the two donor levels, measuring an oscillating relaxation 
rate consistent with acoustic phonon emission in silicon. 

While A1-A1 2 3 rf-SETs are well established as fast 
charge sensors, it is advantageous to integrate the SET 
sensor into the silicon device itself, as has been done with 
silicon-based QPC sensors, since this can improve the ca- 
pacitive coupling to the system being measured and can 
also simplify fabrication. Furthermore, the larger charg- 
ing energies that can be obtained with silicon quantum 
dots, compared with Al metal islands, provides the po- 
tential for increased sensitivity and higher operating tem- 
perature. Figure[29]shows an example of a silicon SET in- 



tegrated adjacent to a Si MOS quantum dot ( |Yang et al. 



2011). In this experiment, Yang and co-workers also em- 
ployed a dynamic feedback technique to keep the SET 
sensor at a point of constant sensitivity, allowing for more 
robust measurements that can tolerate random charge 
displacement events. Podd and co-workers in Cambridge 
also demonstrated a capacitively coupled pair of Si MOS 
quantum dots, in which one of the dots could be used to 



sense the potential of the other (Podd et al. 2010) 



Angus et al. (2007) configured a silicon-based rf-SET 
by using a double-gate structure to induce a Si-MOS 
quantum dot and connecting this within a radiofrequency 
tank-circuit. They demonstrated a charge sensitivity of 
better than 10 -5 e/vHz at a bandwidth up to 2 MHz, 
which compares well with metallic rf-SETs. In their de- 
vice the bandwidth was limited by a high gate resistance, 
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FIG. 29 Non-invasive charge sensing of a Si MOS 
quantum dot using a single electron transistor (SET) 
sensor, (a) SEM device image, showing a Si MOS SET sensor 
(upper device) that is capacitively coupled to a Si MOS quan- 
tum dot (lower device), (b) Transport current Id through 
the quantum dot shows Coulomb peaks as a function of dot 
plunger gate voltage Vpd- The changing potential on the 
dot is detected by monitoring the uncompensated current 7s 
through the SET sensor, which shows charge transfer events 
superimposed on a rising background, due to the coupling of 
the SET to Vpd- This background can be largely removed 
by adding a linear correction (fixed compensation) to the 
SET gate voltage Vps, and then further enhanced by plot 
ting the derivative dls/dVpD- Data reproduced from Yang 
et aZ.I (12011 1. 



but there is no reason why such a structure could not be 
designed to operate at bandwidths above 100 MHz. One 
advantage of a Si-MOS SET compared with its AI-AI2O3 
counterpart is that the tunnel barriers of the Si-MOS de- 
vice are gate controlled, meaning that the resonant fre- 
quency of the tank circuit can be easily tuned to optimize 
its operation. 

For studies of spin dynamics, which can be orders of 
magnitude slower than charge dynamics in silicon, the 
need for high-frequency sensing becomes less critical and 
standard low-frequency (sub-MHz) SET operation can 



be used (Hofheinz et al. 2006a I. Most notably, Morcllo 



et al. (20101 used a Si-MOS SET, similar to the struc- 



ture used by Angus et al. (2007), to detect charge mo- 



tion between the SET island and implanted phosphorus 
dopants, thus enabling single-shot spin readout of an elec- 
tron bound to a phosphorus donor. This experiment is 
discussed further in Section VI. C. 3. 



D. Few-electron quantum dots 

For many years it was difficult to achieve single- 
electron occupation in gated quantum dots, in spite of the 
tunability of such dots. The fundamental problem was 



the difficulty maintaining reasonably fast tunnel rates be- 
tween a quantum and nearby charge reservoirs. A com- 



mon gate design (see, e.g., Waugh et al. (1995)), is shown 



schematically in Fig. |30[a) . As the quantum dot is made 
smaller, by making the gate voltages more negative, the 
tunnel barriers to one or both reservoirs must become 
wider. 
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FIG. 30 Gate design enabling few-electron occupa- 
tion. The gate design in (a) is a natural way to form a 
quantum dot tunnel-coupled to two reservoirs, as shown by 
the arrows. As the dot becomes smaller, however, it is very 
difficult to maintain a high tunnel rate to both reservoirs. 
The gate design in panel (b), based on Fig. 1 of |Ciorga et al. 



(2000), enables a small dot to be coupled to both reservoirs. 



Fig. 30 'b) shows an alternative approach for the for- 
mation of few-electron quantum dots in GaAs, developed 
by the group in Ottawa ( Ciorga et al. 2000 ) . The advan- 



tage of this gate design is that it enables strong tunnel 
coupling to both reservoirs even when the quantum dot 
is small. This gate design is equally useful for gated dots 
in Si, and it was first implemented in a Si/SiGe het- 
erostructure in Sakr et al. ( 2005 1 , enabling observation 



of both Coulomb blockade and charge sensing, but not 
single-electron occupation. 

The challenge to achieving single electron occupation 
in both single and double one-electron dots in Si/SiGe 
has been to bring under control instability in the back- 
ground offset charge of the quantum dots. In 2007 Sim- 
mons et al. (2007) demonstrated single-electron occupa- 



tion in a top-gated, Si/SiGe quantum dot. In that work, 
care was taken to ensure that the doping of phosphorous 
in the modulation doping layer was not larger than neces- 
sary; limiting the doping in this layer appears to improve 
the stability of devices. The primary evidence for single- 
electron occupation was the absence of additional charge 
transitions, as shown in Fig.[28j for a change in gate volt- 
age more than 3.5 times as large as that required to add 
the last observed electron. 

Metal-oxide-semiconductor quantum dots can also ap- 



proximate the few-electron regime (Prati et al. 2011 ). In 



the approach of Xiao et al. (2010b), the depletion gates 



underneath a global accumulation gate form the quantum 
dot. Using an approach analogous to this type of MOS 
Si structure, Borselli and collaborators have shown that 
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single-electron occupation can be achieved in very stable 
Si/SiGe quantum dots when the doping is removed from 

see section 



the structure (Borselli et al. 2011a 



IV.B.4 



A novel approach to achieving single-electron occupa- 
tion was demonstrated by Borselli and colleagues at HRL 
Laboratories ( Borselli et al. 201 lb I . As shown in Fig. 31 



the device structure uses two quantum wells, the lower 
of which is doped. An air bridge is used to apply a posi- 
tive voltage to an isolated, circular surface gate, pulling 
electrons into the upper quantum well. Nearby surface 
gates are negatively biased, enabling the formation of a 
charge-sensing channel in the lower electron layer. Such a 
device forms an extremely symmetric quantum dot that 
is easily tuned to the one-electron charge state. 




FIG. 31 (Color online) Schematic diagram of a few-electron 
quantum dot formed from a Si/SiGe heterostructure with a 
double quantum well and an accumulation gate contacted by 
an air bridge. Inset: SEM microgr aph of the gate region of a 
corresponding device. Figure from Borselli et al. (2011b I. 



Few-carrier occupation can be accomplished even 
in the absence of charge sensing, as demonstrated 
in nanowire-based hole quantum dots for which the 
Coulomb diamonds open to very large gate voltages at 



sufficiently positive gate voltage (Zhong et al. 2005) 



Zwanenburg and collaborators have reached the one-hole 
state in a very small Si quantum dot in a nanowire, en- 



abling them to perform spin spectroscopy (Zwanenburg 



et al, 2009b). The device made use of NiSi contacts, 



in which a Schottky barrier defines the quantum dot, as 
shown in Fig. |32| The few-electron regime was also ob- 
served without charge sensing in planar MOS Si quantum 
dots, thanks to the high degree of tunability of these de- 
vices (|Lim et al] 120111 |2009b|), and in MOSFETs built 



within a pre-industrial Fully Depleted Silicon On Insula- 



tor technology (Prati et al. 2012a). 



E. Spins in single quantum dots 

In the previous sections we have established the evo- 
lution in recent years from the observation of simple lo- 
calization and coulomb blockade to few-electron quantum 
dots in silicon. With the understanding and control of the 
charge side of electrons one can also probe their spins. In 




FIG. 32 (Color online) Transport data showing the last hole 
in a Si nanowire based quantum dot. Inset: SEM image of 
the device showing the NiSi contacts and a Cr/Au side gate, 
device. Figure from Zwanenb urg et al\ ( |2009b| . 



this section we first discuss experiments on ground-state 
and excited-state magnetospectroscopy in silicon quan- 
tum dots. The existence of valleys in silicon make the 
spin filling non-trivial: the configuration and mixing of 
valleys and orbits determines how electrons will consec- 
utively occupy the available spin-up or down states. 



1. Spin-state spectroscopy 

The most straightforward methods of measuring elec- 
tron spin states in quantum dots are ground-state and 



excited-state magnetospectroscopy (Hanson et al. 2007). 



Excited-state magnetospectroscopy allows observation of 



spin excited states at a fixed magnetic field (Cobden 



et al. 1998), as long as the Zeeman energy can be re- 



solved. Four experimental demonstrations in silicon sys- 



tems are bottom-up Si and SiGe nanowires (Hu et al. 



2011 Roddaro et al. 2008 Zwanenburg et al. 2009b) 



and SiGe nanocrystals (Katsaros et al. 2010 1, see Fig. 33 



When the spin-excited state is measured at different mag- 
netic fields, one can extract the g-factor by plotting the 



Zeeman energy versus magnetic field, see Fig. 33 'b). The 



first two holes in a Si nanowire quantum dot were found 
to have a g-factor of 2.3 ± 0.2 in perpendicular magnetic 
field. In SiGe nanocrystals and nanowires the g-factor 
is anisotropic: the results in Fig. 33 'c) show g-factors of 
mi = 1.21 and g± = 2.71 for respectively parallel and 
perpendicular field. 

In case of ground-state magnetospectroscopy, the spin 
filling is investigated by measuring the magnetic field 
dependence of the electrochemical potential fiN, which 
is by definition the energy required for adding the N th 
electron to the dot. The slope of /ijv(-B) is given by 
= —gnB&-S tat {N), where g is the g-factor, the 
Bohr magneton [is = 58 /^eV/T and AS to t(N) is the 
change in total spin of the dot when the N th electron is 
added (Hada and Eto 2003). The electrochemical po- 



tential has a slope of +c//is/2 when a spin-up electron is 
added, whereas addition of a spin-down electron results 
in a slope of — <7/is/2. The rate at which ^jv changes with 
magnetic field thus reveals the sign of the added spin. 
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FIG. 33 (Color online) Excited-state magnetospec- 
troscopy in Si quantum dots, (a) Zeeman splitting at 
the 0-1 and 1-2 transition in a few- hole Si nanowire quan- 
tum dot and (b) the corresponding magnetic field dependence 
of the Zeeman energy, data from Zwanenburg et al\ ( 2009b) 
(c), Anisotropic g-factors in SiGe nanocrystals, and (d) the 
corresponding excited-state magnetospectroscopy, data from 
|Katsaros et al.\ ( |2010| >. 
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FIG. 34 (Color online) Ground-state magnetospec- 
troscopy. Three examples of even-odd hole spin filling in (a) 



2005) 



a many-hole Si nanowire quantum dot (Zhong et al. 
(b), a few- hole Si nanowire quantum dot ( Zwanenburg et aT] 
2009b) and (c), a ma ny-hole Ge/Si nanowire quantum dot 
(IRoddaro et aZ.1 12008 1. 



2. Spin filling in valleys and orbits 



Rokhinson et al. (2000) were the first to observe the 



theoretically expected slopes in multiples of <?/zb/2 in an 
n-type Si quantum dot. They show the peak shift with 
magnetic field of 29 electrons entering the dot, and more 
detailed measurements on two sets of Coulomb peaks 
with slopes of ±l/2g/is and ±3/2g/is. The charge tran- 
sitions display an unexpected large number of kinks at 
which the slope changes sign, and thus the spin state as 
well. They conclude that the spin filling is inconsistent 
with a simple picture of non-interacting electrons in four 
single-particle levels. Later reports are more straightfor- 
ward to interpret and will be discussed below. 

The spin filling of holes has been investigated in 
nanowire quantum dots. 



The even-odd spin filling as observed in p-type silicon 



quantum dots (see Section IV.E.l) is not very different 



from similar devices in other material systems. However, 
the valleys in the silicon conduction band make the spin 
filling of electrons non-trivial. Valley physics in silicon 
has been studied extensively both theoretically (|Culcer 
et al.[ |2010a|b| |Friesen and Coppersmith) |2010| |Saraiva 



et al. 



2009 



2010: 



2011) and experimentally (Fuechsle et al. 



Goswami et al. 2007 Kocster et al. 



1997 



Kohler| 



and Roos[ |l979l |McGuire et al.\ [2010| |Nicholas et al. 



1980 Pudalov et al. 1985 Takashina et al. 2006) 



As discussed in Section IIII.B.2I a 2-dimensional elec- 
tron gas has two T-valleys, separated by the valley split- 
ting Ey, see Fig. [7j A finite valley splitting influences 
the spin filling as observed in ground-state magnetospec- 



In 2005 

found alternating spin-up and spin-down holes in a many- 
hole quantum dot. The magnetic field evolution of 
the positions of eight consecutive Coulomb peaks in 
Fig. [34^ a) reveals alternating slopes of ±<7/Lts/2, with an 
extracted g-factor of 2±0.2. The few- hole regime dis- 
played similar spin filling of the first four holes in an 



Zhong et al.\ (|2005|) troscopy ( |Hada and Eto| |2003j): the first electron is al- 



ways a spin-down, yielding a slope of the corresponding 
Coulomb peak of — g/^/2, see the experiment by Lim 



empty dot (Zwanenburg et al. 2009b), see Fig. |34|b). 
The even-odd filling suggests that the degeneracy of 
heavy and light holes is lifted due to strain and con- 
finement effects; see, for example, calculations based on 



density functional theory ( Leu et al. 2006 Sorokin et al. 



2008) and tight-binding models (Buin et al. 2008 Ni 



et al. (2011 1 in Fig. |35p. The kink in the second Coulomb 
peak (marked 2a) at ~0.86 T is caused by a sign change 
of the N = 2 ground-state spin: at low magnetic field 
(before the kink), the second electron fills the quantum 
dot with a spin-up. As the magnetic field is increased, 
the sign of the second electron spin changes from up to 
down at B ~ 0.86 T. 

When the valleys and orbits are mixed (section 
III.B.2 1 , there are no pure valleys or pure orbits, and 



quet et al. 2006). SiGe nanowires have been shown to 



exhibit the same spin filling, see (Roddaro et al. 2008) 



and Fig. 34 'c 



the lowest available levels are referred to as valley-orbits. 
The sign change can then be explained with a simple 
model where the two lowest valley-orbit levels are sepa- 



rated by the valley-orbit splitting AEyo , see Fig. 35 'b) 
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FIG. 35 (Color online) Spin filling in valleys in a planar 

MOS Si quantum dot. (a) Magnetospectroscopy of the 
first two electrons entering the quantum dot. The circle 2a 
marks a kink in the second Coulomb peak at ~0.86 T. The 
arrows in the boxes (VOi for valley-orbit 1 and VO2 for valley- 
orbit 2) represent the spin filling of electrons in the quantum 
dot. (b) For B < 0.86 T, the first two electrons fill with 
opposite spins in the same valley-orbit level (left panel). The 
Zeeman energy at the kink is equal to the valley-orbit splitting 
(0.10 meV). Data reproduced from Lim et al. (2011 1 



At zero magnetic field, the first two electrons fill with op- 
posite spins in valley-orbit level 1. When a magnetic field 
is applied, the spin-down and spin-up states are split by 
the Zeeman energy Ez- Above 0.86 T the spin-up state 
of valley-orbit level 1 (VOi) is higher in energy than the 
spin-down state of valley-orbit level 2 (VO2) and it be- 
comes energetically favored for the second electron to oc- 
cupy the latter, i.e. VO2. At the kink the valley-orbit 
splitting equals the Zeeman energy, which is 0.10 meV at 
0.86 T. Comparable kinks were reported simultaneously 
in accumulation mode Si/SiGe quantum dots, yielding 



valley splittings of 0.12 and 0.27 meV (Borselli et al. 



2011b). In 2010, the absence of kinks in the ground- 
state magnetospectroscopy of a planar MOS Si quantum 
dot was explained as a result of a large exchange energy 
and an unusually large valley splitting of 0.77 meV (Xiao 



et al. 2010b |. 



F. Double quantum dots 

Like their counterparts in the Ga-AlGaAs material sys- 
tem, double quantum dots in silicon represent the natu- 
ral extension from a semiconductor 'artificial atom' to an 
'artificial molecule'. As outlined in the previous sections, 
it took until around 2006 for low-disorder silicon-based 
quantum dots to be produced with reasonable repeata- 
bility. Correspondingly, this is also when the first demon- 
strations of double quantum dots in silicon began to be 
reported. 



1. Charge-state control 

One of the earliest reports of silicon double dot oper- 



ation was by Gorman et al. ( 2005 ) , who formed an iso- 



lated double dot by etching a thin (35 nm) layer of bulk- 
phosphorus-doped silicon (Si:P) in a SOI substrate. They 
also integrated a nearby SET, again formed by etching 
the Si:P layer, which they used to monitor charge trans- 
fer in the double dot. By rapidly pulsing a nearby control 
gate they observed oscillations in the charge state of the 
double dot, as a function of pulse length, which they in- 
terpreted as coherent oscillations between the (n, m) and 
(n-1, m+1) charge states of the double dot. Because of 
the high electron numbers in the dots resulting from the 
degenerative doping, and the difficulty of controlling the 
dots size via the etching process, this type of dot struc- 
ture has not progressed significantly since this time, and 
most studies of silicon quantum dots are now based on 
dots induced in undoped silicon layers. 

The starting point for any experimental study of a dou- 
ble quantum dot is the determination of its charge state 
(Ni, AT 2 ) as a function of at least two gate voltages V g \ 
and V g 2 controlling the system. Here N\, (N2) is the elec- 
tron occupancy of dot-1 (dot-2). By directly measuring 
the transport current / through the dot, or by measur- 
ing the local electrostatic potential of the system using a 
nearby QPC or SET charge sensor, one can plot a charge 
stability map as a function of V g \ and V g 2- |Van der Wiel 
et al. ( 2003 ) provide an excellent review of semiconduo 



tor double quantum dots and Fig. 36 a-c) taken from this 
review depicts the charge stability maps expected for dif- 
ferent inter-dot coupling strengths. Figure 36 'b) shows 
a map at intermediate inter-dot coupling, where one ob- 
serves a characteristic 'honeycomb' like structure. The 
points on the map where three different charge states are 
degenerate in energy are known as 'triple points' and it 
is only at these points where a transport current / can 



flow through the system, from source to drain (Fig. 36 i) . 
If we then apply a source-drain bias Vsd across the dou- 
ble dot these triple points expand to take on a triangular 
shape, as shown in Fig. |36^e). It should be noted that 
while in a transport measurement only the triple points 
(or bias triangles) can be observed, when charge sensing 
is employed one can also directly detect the transition 
lines between charge states. Das Sarma et al. (2011) 



have recently applied a Hubbard model approach to de- 
termine the charge stability diagrams for silicon double 



dots, showing excellent agreement with experiments ( Lai 



et al. 2011 Simmons et al. 2009) 



The first clear demonstration of a double quantum dot 



in nominally undoped silicon (Fujiwara et al. 2006) used 
a double-gated silicon nanowire formed from a SOI sub- 
strate, with three independently controlled barrier gates. 
This device structure was previous described in Section 



IV. B. 5 and is depicted in Fig. [26[a,b). By varying the 



voltage Vlgc on the central barrier gate, Fujiwara and 
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FIG. 36 Schematic stability diagrams for a double dot 
system. Maps are shown for (a) small, (b) intermediate, and 
(c) large inter-dot coupling. The equilibrium charge on each 
dot in each domain is denoted by (Nl ,N2). (e) Region within 
the dotted square of (b), corresponding to the unit cell of the 
double dot stability diagram, at finite bias voltage. The solid 
lines separate the charge domains. Classically, the regions of 
the stability diagram where current flows are given by the 
gray triangles. Reproduced from Van der Wiel et al. (2003 k 



co-workers were able to gradually tune the system from 
one large dot to two well-isolated dots, as evidenced 
by the charge stability plots shown in Fig. 37 b-e). In 
Fig. [37^b) the central gate voltage Vlgc is tuned to 
minimize the tunnel barrier between the dots, forming 
one large dot. As the outer barrier gate voltages Vlgs 
and Vlgd are varied, the transport current I plotted in 



Fig. 37 'b) shows Coulomb blockade oscillations as a func- 
tion of the addition voltage (Vlgs + Vlgd), which can 
be compared with Fig. [36]^c) , indicating that one large 
dot is present. However, as the tunnel barrier height is 
increased by tuning Vlgc, two separate dots form, as 
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FIG. 37 (Color online) Evolution from a single dot to 
a double quantum dot in a gated silicon nanowire 
device, (a) Equivalent circuit, (b-e) Contour plots of the 
drain current as a function of the outer barrier gate voltages 
Vlgs and Vlgd- The central barrier gate voltages used were 
Vlgc = -0.75 (b), -1.13 (c), -1.18 (d), and -1.284 V (e). 



Data reproduced from Fujiwara et al. (20061 



revealed in Figs. 37 c-e) with the gradual formation of 
a honeycomb-shaped map of I (Vlgs, Vlgd), consistent 
with that in Fig. 36 'b) . Gate-tuneable double quantum 
dots based on etched silicon nanowires have since also 



been reported by other groups (e.g. Pierre et al. (20101). 

Epitaxially grown nanowires have also been configured 
as double quantum dots. Hu and co-workers from Har- 



vard (Hu et al. 2007) used a Ge/Si core/shell nanowire, 
as described in Section IV.B.l, with a number of surface 
metal gates to demonstrate a highly-tunable double dot 
device, as shown in Fig. 38 Figs. 38 ^b, c) show the evo- 



lution of the charge stability map from a strongly coupled 
single dot to a double dot as the central barrier gate is 
tuned. 

Using a gated etched nanowire device, |Liu et aL\ 



(2008a) were able to explore the excited state energy 
levels within a double quantum dot, which in this case 
was formed using two barrier gates and the presence of 
a third barrier created by local disorder. Figure [39^ a) 
shows the charge stability diagram for this double dot 
in the presence of a source-drain bias, which transforms 
each triple point into a 'bias triangle', as described in 
Fig. 36 'e) . By mapping the bias triangle in more detail, 
Liu and co-workers were able to observe resonant tun- 
nelling through excited states of the double dots. Using 



a planar MOS structure, similar to that in Fig. 24'a,b), 



Lim et al. ( 2009a ) were also able to observe excited state 



transport through a double quantum dot, this time con- 
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FIG. 38 (Color online) Gate tunable double quantum 

dots, (a) SEM image of a Ge/Si nanowire-based hole quan- 
tum dot. The Ge/Si nanowire at top (white in image) is gated 
by metal gates to form a double dot. (b,c) Charge stability 
maps of the conductance as a function of plunger gate volt- 
ages, (d) SEM image of an electron quantum dot defined 
by electrostatic top gates in a Si/SiGe heterostructure. (e) 
Charge sensing measurement showing the difference in the 
charge detection signal from the dot furthest from the QPC 
(4 small steps in /qpc) and the dot closest to the QPC (single 
large step) as a function of gate voltage, (e) Two-dimensional 
plot of the charge sensing current showing the sequential ad- 
dition of electrons to the left and right dots. Data in (a-c) 
repr oduced from|Hu et al.\ ( |2007| ); data in (d-f) reproduced 
from Simmons et al. (2009 1. 



trolled using three independent barrier gates. Figure 
39 'b) shows a pair of bias triangles for two triple points, 
each showing structure in the source-drain current Isd, 
that is further revealed in Fig. 39 c) , which shows a line 



trace of Isd along the detuning axis e. 
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FIG. 39 (Color online) Bias spectroscopy of silicon dou- 
ble quantum dots, (a) Stability map with a source-drain 
bias Vsn = 1 mV for a silicon nanowire double dot, depicted 
in Fig. |26| obtained by plotting source-drain current / as a 
function of two barrier gate voltages. The triple points have 
clearly evolved into bias triangles, (b) Bias triangles for two 
triple points at Vsd — 1 mV, obtained in a Si MOS double 
dot. (c) Line trace of Isd, taken along arrow in (b), show- 
ing resonances corresponding to excited states in the double 
dot. Data reproduced from Liu et al. (2008a I and Lim et al 
( |2009a| ). 



highly tunable double quantum dots in a gated Si/SiGe 
device. The device structure shown in Fig. |38^d) also 
incorporated a QPC charge sensor, which enabled map- 
ping of the double dot charge stability as a function of 
the gate voltages Vj, and Vr controlling the two dots 
(see Fig. |38{e,f)). They were able to tune the inter-dot 
coupling by control of the central gate voltage Vm- The 
same group were able to demonstrate depletion of a dou- 



ble quantum dot to the single electron level ( Thalakulam 



et al. 20101. The data in Fig. 40 b) shows a charge sta- 
bility map of the double dot, measured using the QPC 
sensor. The lack of charge transitions in the lower left 
quadrant of this map demonstrates control of electron 
number down to the (0, 0) charge state. 

Occupation down to the (0,0) charge state in a double 
dot has also been recently demonstrated in an undoped 



Si/SiGe heterostructure device (Borselli et al. 2011b), in 
which a two-layer gate structure analogous to that used 



for Si MOS dots (see Fig. 23 i,b) was used, resulting in a 
very stable system. Very recently, a few-electron double 
dot has been demonstrated in an isotopically-enriched 
Si quantum well (Wild et al.||2012| |. 
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FIG. 40 Single-electron occupancy in a Si/SiGe dou- 
ble quantum dot. (a) SEM of the device, (b) Charge sta- 
bility map of the double dot, obtained by plotting the QPC 
charge sensor output as a function of the control gate volt- 
ages Vl and Vr. The charge configurations (n, m) are marked, 
showing depletion to the (0, 0) state. Data reproduced from 
Thalakulam et aZ. I poTot) . 



2. Spin transport in double quantum dots 

As discussed in Section V.D, in tightly-confined quan- 
tum dots the singlet-triplet exchange energy for an elec- 
tron pair can become appreciable. In a double quantum 
dot system this can lead to an effect known as Pauli spin 
blockade, where transport through the double dot is de- 
pendent upon the spin state of the electron. This phe- 
nomenon was first observed in 2002 in a GaAs/AlGaAs 
double quantum dot (Ono et al. 2002). In 2008 Liu and 



Simmons et al. (2009) demonstrated one of the first 



co-workers observed Pauli blockade in a silicon nanowire 
device similar to that shown in Fig. |26[ a,b), in which a 
double quantum dot was formed using two barrier gates 
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and a third barrier created by local disorder (Liu et al. 



2008b) 
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FIG. 41 (Color online) Pauli spin blockade in a sili- 
con MOS double quantum dot. (a) SEM image and (b) 
cross-sectional schematic of the Si MOS device. Gates LI and 
L2 induce electron reservoirs at the Si/Si02 interface, while 
barrier gates B1-B3 define the double dot potential. Plunger 
gates PI and P2 control the occupancy of each dot. (c,d) 
Current /jc as a function of Vpi and Vp2 for B — T. (c) 
For Vsd = +2.5 mV, the ground state and excited states of 
a full bias triangle are observed. The current flows freely at 
the S(0,2)-S(l,l) transition, as illustrated in the box marked 
by the dot. (b) The same configuration at Vsd = —2.5 mV. 
Here the current between the singlet and triplet states is fully 
suppressed by spin blockade (box marked by star), (e) The 
measured singlet-triplet splitting Ast, plotted as a function 
of magnetic field B. Figure from Lai et al. (2011 1. 



Lai et al. (2011) demonstrated Pauli blockade in a Si 
MOS double quantum dot formed using an AI-AI2O3 
multi-layer gate stack, similar to that discussed in Section 
IV. B. 4 Their device structure, shown in Fig. |4l|a,b), 



incorporated three barrier gates (Bl, B2, B3) and two 
'plunger' gates (PI, P2) which controlled the occupancy 
of the two dots. Pauli blockade occurs in one bias di- 
rection only, as depicted in Fig. [41] The phenomenon 
is revealed experimentally via the difference between the 
bias triangles for positive and negative source-drain bias 



Vsd, observed in the charge stability map as a function 
of the two plunger gate voltages. 



For Vsd > (Fig. 41 ;) a pair of overlapping complete 
bias triangles are observed. Resonant transport through 
the ground and the excited states in the double dot oc- 
curs when the states within the dots are aligned, lead- 
ing to current peaks that appear as lines parallel to the 
triangle base. The non-resonant background current in 
the triangle is due to inelastic tunneling. The non-zero 
current throughout the triangular region indicates that 
electrons from the reservoir can tunnel freely from the 
5(0,2) singlet state to the 5(1,1) singlet state, as de- 
picted in the cartoon (red box in Fig. 41 ). For Vsd < 



the current is suppressed in the region bounded by the 
dashed lines in Fig.|4l|d). The suppression arises because 
the transition from T(l, 1) to S(0, 2) is forbidden by spin 
conservation during electron tunneling. Once the T(l, 1) 
triplet state is occupied, further current flow is blocked 
until the electron spin on one dot reverses its orientation 



via a relaxation process (green box in Fig. 41 ). 

In Fig.|4l|d) it is possible to discern some non-zero cur- 
rent at the bottom of the bias triangle. This 'leakage cur- 
rent" in the spin-blockade regions has been identified as 



resulting from a spin-flip co-tunneling mechanism ( Coish 



and Qassemi |201l||Lai et al. 2011 Qassemi et al. 2009), 



where a spin- up electron from one of the reservoirs swaps 
with a spin-down electron in one of the dots. This effect 
has also been observed in a double dot formed from an 



etched silicon nanowire device (Yamahata et al., 2012). 



Note that in this experiment (Lai et al. 20111, and 



also that performed in the nanowire device (Liu et al. 



2008b), the electron occupancy in each dot was of order 
10 or more, and so the labels (1,1) and (0,2) refer to the 
effective electron occupancy, whereas the true electron 
occupancy is (ra + mo, n + no). Pauli blockade for two- 
electron singlet and triplet states therefore occurs when 
the total electron spin of each dot is zero in the (ra , n ) 
state. More recently, Borselli et al. (2011a) have demon- 



strated a Si/SiGe double dot that exhibits Pauli blockade 
in the true (l,l)-(0,2) limit. The same group have since 
used this structure to demonstrate coherent oscillations 
between singlet and triplet states of the double dot sys- 
tem ( |Maune et al. 2012), as discussed further in Section 
IVI.C.4I 

By applying a magnetic field B, it is possible to modify 
the singlet-triplet splitting Ast, defined as the energy 
difference between the blockaded ground state S(0, 2) 
and the excited state T!_(0, 2). In a magnetic field there 
are four accessible spin states: the singlet S; and three 
triplets TL, T and T+, corresponding to S z = -1, 0, +1. 



Lai et al. (20111 studied the singlet-triplet splitting by 
mapping the bias triangles in the spin blockade regime 
at increasing magnetic fields B = — 8T. They found 
that the splitting Ast decreased linearly with increasing 
B [Fig. |4T|e)], as expected, since the triplet states split 
linearly by the Zeeman energy, Ez — ±S z g(iBB, where 
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Hb is the Bohr magneton and Sz is -1, 0, +1. Here, 
a linear fit through A$t(B) yielded a Lande g- factor of 
2.1 ± 0.2, consistent with electrons in silicon. 



zero current were also observed (see Figs. 42 |d,e)). These 
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FIG. 42 (Color online) Spin blockade and lifetime en- 
hanced transport in a Si/SiGe double quantum dot. 

(a) Measured, and (b) schematic, charge stability map of cur- 
rent / through the double dot, with a source-drain bias of 
Vsd ~ +0.2 mV. The dotted trapezoids in (a) and (b) mark 
the zero current regions due to spin blockade, as depicted 
in the schematics in (c). (d) Measured, and (e) schematic, 
charge stability map of current / with a source-drain bias of 
Vsd = —0.3 mV. In this bias direction there is no blockade 
and current flows throughout the entire bias triangle, how- 
ever, additional tails are observed due to lifetime enhanced 
transport, as depicted schematically in (f) and described in 
the text. Figure from Shaji et al. (2008). 



Spin transport in silicon double quantum dots can also 
be strongly affected by the relative values of the spin life- 
times and the various tunnel rates of the system. For 
example, it is sometimes possible to observe transport 
through a double dot in gate-space regions where cur- 
rent would normally be blockaded, a phenomenon that 
has been termed lifetime enhanced transport (LET), ob- 
served in a Si/SiGe quantum dot structure ( |Shaji et al. 



2008). Shaji and co-workers observed this effect in a de- 
vice similar to that depicted in Fig. [28j which was origi- 
nally configured to operate as a single quantum dot. Un- 
der certain gate bias conditions a double quantum dot 
could be formed that exhibited a charge stability map 
consistent with occupancies down to the single electron 
level, in particular providing access to the (1, 1) and (0, 2) 
states. In the vicinity of the (1, 1) — (0, 2) charge transi- 
tion, a positive source-drain bias (Vsd = +0.2 mV) led 
to bias triangles that exhibited a zero current region due 
to Pauli blockade (see Figs. |42^a,b)), similar to that ob- 
served in Fig. |4Tj 

Under negative bias (Vsd = -0-3 mV) the entire bias 
triangle exhibited a transport current, as expected, but 
outside the bias triangles additional "tail" regions of non- 



can be understood with reference to the green box in the 
schematic diagram in Fig. |42^f) and by considering the 
tunnel rates to the relevant states in the double dot (ISim 



mons et al. 2010). Typically, once the S(2,0) state is 



loaded from the source reservoir (with rate Tls), it would 
remained locked in position due to Coulomb blockade, 
and an electron could only pass to the drain via a co- 
tunneling process (with unloading rate T^). If, however, 
the 5(2,0) loading rate is much slower than the unload- 
ing rate (Tls « ^s) and the triplet-singlet relaxation 
rate r^s from T(2, 0) to S(2, 0) is also much slower than 
Tg, then the most likely transport pathway is via the 



T(2,0) and T(l, 1) triplets, as shown in Fig.ggf), lead- 
ing to a non-zero current. The upper triangle in the data 
and schematic diagram, known as the hole triangle, also 
shows LET behavior. The direction of the "tail" in gate 
voltage space for the hole triangle is different than that 
for the electron triangle, a phenomenon that can be ex- 



plained by a spin-flip co-tunneling process ( |Koh et al. 
2011). This type of spin-lifetime enhanced transport can 



occur in silicon quantum dots in part because of the very 
long spin lifetimes present, as discussed further in Section 
VI.A. 



V. DOPANTS IN SILICON 

A. Dopants in silicon transistors 

1. Early work: mesoscopic silicon transistors 

Low temperature transport experiments in silicon 
transistors have been used since the 1980's to perform 
spectroscopic measurements of dopants and defect states. 
Devices reached a sufficiently small length scale about 25 
years ago to observe mesoscopic transport phenomenon, 
see the review by Fowler et al. (1988). In these devices 
the conductance, G, was found to fluctuate as a func- 
tion of the gate voltage, Vq around the threshold, as 
shown in Fig. [43k. The strongly fluctuating pattern, re- 



producible within the same cooling cycle, did not origi- 
nate from electrical noise but from the presence of a finite 
number of dopant or defect states in the channel. The 
important length scales needed to understand this phe- 
nomenon are the localization length, £ and the device 
dimensions (channel length L and channel width W). As 
device dimensions approached ^100 nm in size, compa- 
rable or less than the localization length, only a limited 
number of defect or dopant states contributed to the cur- 
rent in contrast to the much larger device dimensions 
before the 1980's where these fluctuations were always 
averaged out. 

Three major conduction processes are known to con- 
tribute to the conductance of such small transistors, as 
shown in Fig. [43b. Firstly, there is thermally activated 
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hopping conduction, where the transported electrons hop 
via several dopant states from source to drain. This 
type of transport decreases exponentially with temper- 
ature T, i.e. G oc exp(— AE/k B T), where AE is the 
energy difference between the localized states of the dom- 
inant (most resistant) hop and fee is the Boltzmann con- 
stant. The second process is direct tunneling, which 
scales with the barrier height E-q and L roughly as 
G oc exp — yj2m* E-qL / h 2 . Although typically irrelevant 
in the early 1980's, his conduction mechanism is play- 
ing an increasingly important role in today's nano-scale 
transistor operation. The last and increasingly impor- 
tant conduction mechanism is tunneling through a single 
defect. Due to the large number of conductance fluctu- 
ations, as depicted in Fig. [43^,, the identification of each 
state with a particular defect or dopants in the chan- 
nel region was not possible. Dopants are not the only 
sources of disorder that cause localized states in MOS- 



FETs as discussed in Peters et al. (1998) and Sanqucr 



et al. (2000). However, these early measurements repre- 
sent the first observation of mesoscopic physics in silicon 
MOSFET devices, and show how low temperature trans- 
port data offer a tool to electrically access dopant states 
in the channel region. 
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FIG. 43 Conductance in micron-scale silicon MOS- 
FETs. (a) Typical low-temperature conductance pattern of 
a 1980's generation MOSFET around the threshold regime. 
The strongly oscillating but chaotic pattern that appears at 
low temperature is associated with localized states in the 
channel region, (b) Schematic representation of the three ma- 
jor conduction mechanisms through the channel. From [Fowler| 
\eTal\ [19881 



age were observed, caused by the statistical fluctuation 
in the number of dopants in the channel as a result of 
the random Poisson distribution during doping. These 
results challenged the conventional understanding and 
modeling of micron-silicon devices where continuous ion- 
ized dopant charge with smooth boundaries and inter- 
faces had previously been assumed. Now the granularity 
of the electric charge and the atomicity of matter in- 
troduced substantial variation in individual device char- 



acteristics, as shown in Fig. 44 (Asenov et al. 2003). In 



particular the variation in number and position of the in- 
dividual dopant atoms in the active region of MOSFETs 
were found to make each transistor microscopically dif- 
ferent, introducing significant variations from device to 
device. 




FIG. 44 The importance of discrete dopants in nano- 
scale MOSFETs. (a) The transition from continuously ion- 
ized dopant charge and smooth boundaries and interfaces to 
(b) a 4-nm MOSFET where there are less than 10 Si atoms 
along the channel. From [Asenov et al.\ |2003| 

Recent advances in single ion implantation using a fo- 
cussed ion beam source have shown that device to device 
fluctuations can be suppressed by ordering the dopants 



within the channel. ( |Shinada et al. 2005). The benefit of 
ordering the dopants was to create a homogeneous poten- 
tial distribution in the channel resulting in the formation 
of a uniform current path. The uniformity in the channel 
lowers the voltage required to open the channel from the 
source to the drain allowing for an earlier turn on and 
reduction of the threshold voltage. This contrasted to 
devices with a random distribution of dopants where the 
non-homogeneous potential could lead to a block in the 
current path. 



B. Single dopant transistors 

1. The demand for single dopant architectures 



2. Nano-scale transistors 

Following this early work the purity of silicon MOS- 
FETs steadily continued to improve with a concomitant 
decrease in device size, until the point where discrete 
impurities clearly started to show up in device trans- 



In parallel with the increasingly important role of in- 
dividual dopants in classical silicon CMOS devices has 
been the advent of their importance in quantum com- 
putation. In 1998 Bruce Kane introduced the concept 
of using the nuclear spins of individual donor atoms in 
doped silicon electronic devices as quantum bits or qubits 



port properties as they reached the nano-scale (Mizuno (Kane 1998). Originally Kane envisioned that the quan- 
et al. 1994). Here, fluctuations in the threshold volt- turn information could be stored in the state of the 31 P 
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nuclear spin and accessed by the electron-nuclear hyper- 
fine coupling. Figure 11 a) shows the basic two qubit 



unit cell of the Kane proposal based on an arrays of P 
donors beneath the silicon surface. The addition of a 
group V phosphorus donor to the silicon crystal results 
in electron states close in energy to the conduction band 
but weakly bound to the donor site at low temperatures. 
The electron has spin S = 1/2, whilst the host silicon 
has stable 1 = isotopes. It is important to isolate the 
qubits from any degrees of freedom that may lead to de- 
coherence. Recent results have shown that 28 Si can be 



isotopically refined to a level of 99.98% (Tezuka et al. 



2010 ) making it a good choice for the host material. 
Any proposal for a quantum computer must meet the 



so-called DiVincenzo criteria ( DiVincenzo 2000 ) , includ- 
ing the ability to operate on individual qubits, couple 
qubits into quantum logic gates, read-out the informa- 
tion encoded on the qubit, and to be free of environmen- 
tal effects that destroy qubit entanglement. In Kane's 
architecture control of the qubit states is achieved by a 
combination of gates and globally applied a.c. magnetic 
field. The gates above the donors, labeled A gates control 
the strength of the hyperfine interaction and therefore the 
resonance frequency of the nuclear spins beneath them. 
The J gates are used to turn on and off the electron- 
mediated coupling between the nuclear spins. Finally 
the nuclear spins can be flipped by a resonant globally 
applied a.c. magnetic field. Read-out of the final spin 
state is achieved by spin-to-charge conversion and detec- 
tion through spin-dependent tunneling to a doubly oc- 
cupied donor state. This donor state is a singlet with a 



second electron binding energy of 1.7 meV (Larsen 1981 



Larsen and McCann 1992 1. Consequently by applying 



a differential voltage between the A gates charge motion 
between the donors can only occur if the electrons are 
in a singlet state. This charge motion can then be de- 
tected using a sensitive electrometer nearby, such as a 
single-electron transistor (SET). 

Subsequently a number of proposals for encoding and 
manipulating quantum information based on donor spin 
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ward. Electron spins bound to donor nuclei are particu- 
larly attractive since they have exceptionally long coher- 
ence times and relaxation times relative to the time scales 



for the control of the quantum state (Hill et al. 2005 1 . 
The electron spin coherence time of a P donor is T 2 > 60 
ms at T = 6.9 K in isotopically pure 28 Si (Tyryshkin 



et al. 2003). These times are currently limited by the 
presence of 29 Si which causes spectral diffusion due to 



the dipolar fluctuations of nuclear spins (Witzel et al. 



2005 ) . Dipolar fluctuations in the nuclear spins give rise 
to a temporally random effective magnetic field at the lo- 
calized electron spin leading to irreversible decoherence 



(i.e a T2 process). Isotopic purification in silicon systems 
can in principle overcome this limiting process. 

These developments have lead to the proposal of a com- 
plete 2D donor-based architecture, incorporating trans- 



port, for a truly scalable design in silicon (Hollenberg 



et al. 2006). This design has considered the limita- 



tions and constraints posed by the sensitivity of the ex- 



change interaction due to donor placement (Koiller et al. 



2002a), high gate densities required (Copsey et al. 2003), 



spin read-out based on spin-charge conversion (Kane 



1998 ) and the communication bottleneck for linear near- 
est neighbor qubit arrays. A buried array of ionized 
donors provide pathways for coherent transport of elec- 
tron spins for in-plane horizontal and vertical shuttling 
(dashed border sections) of qubit states into and out of 
the interaction zones. The overall gate density is low 
compared to the original Kane version since coherent spin 
transport is achieved adiabatically, lowering the barriers 
between donors in a well defined sequence to effect co- 
herent transport by adiabatic passage (CTAP) without 



populating the intervening donors (Greentree et al.\ 2004 



Rahman et al. 2010). Logic gates are carried out in in- 
teraction zones with the A and J gates for electron spin 
qubit control and these are distinct from the qubit stor- 
age regions. The design allows space for local B-field 
antennae and SET read-out devices. The introduction 
of coherent spin transport to donor quantum computing 
provides a means to consider scalable, fault tolerant ar- 
chitectures. 

The use of single donor atoms in silicon as qubits 
has demanded tremendous advances in single atom fab- 
rication and engineering. Donor separations of order 
~ 10 — 20 nm are required to ensure significant cou- 
pling between neighboring spins. Currently this is at 
the limit of what is technologically achievable. Indeed 
in the original critique of Kane's paper by |DiVincenzo 
( 1998 ) he recognized that the fundamental and engineer- 



ing obstacles to implementing the scheme were vast stat- 
ing "At the time no existing materials-preparation tech- 
nology will place an array of individual phosphorus atoms 
at desired spots in the interior of a perfect crystal, let 
alone systems free from defects in the semiconductor and 
the overlying oxide layer". Despite these concerns there 
have been concentrated efforts internationally to realize 
a donor-based qubit architecture resulting in a plethora 
of experiments of transport in nano-scale doped quantum 
dots and donors. Significant to these results has been the 
different technologies developed to fabricate donor-based 
devices. 



2. Single dopants in MOS-based architectures 

In this sub-section we discuss single dopant transport 
in ultra-scaled MOSFET structures based on randomly 



and deterministically doped devices. The Kane (1998) 
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proposal sparked the interest in single dopants and small 
MOSFETs were quickly identified as devices that should 
easily allow observation of single dopant transport. Tabe 



et al. (2010) confirmed the impact of dopants on the po- 
tential landscape of a FET and discussed the evolution 
from many dopant to single dopant transport. In the 
newer generation of CMOS devices the issue of random 
device fluctuations is circumvented in a more straightfor- 
ward approach, namely by the use of undoped channel 
FETs. However even in these newest generation of pro- 
totype FETs, fluctuations in device characteristics are 
still evident due to the presence of only a few down to a 



single unintentional dopant(s) in the channel region (Col- 



inge et al. 20071. By now, several groups have reported 
transport through a random, single dopant in a three ter- 
minal configuration ( Calvet et al. 2007a|b Pierre et al. 



2010 Sellier et al. 20061, including microwave assisted 



transport (Prati et al. 2009). Recently this work has 
been extended to double gate structures and 2 8Si de- 



vices (Lo et al. 2009 Roche et al. 2012). Roche et al. 



|2012| even demonstrated controlled sequential tunneling 
through two donors. In all these experiments, the num- 
ber of electrons bound to the dopant atom could be con- 
trolled by the gate electrode. Such spectroscopic trans- 
port experiments reveal vital information on the orbital 
levels, the charging energy and the binding energy of the 
dopant atom and the spin configuration of the bound 
electron(s) (Lansbergen et al. 2011 2008). They thus 



form a powerful characterization tool in the development 
of single-dopant structures in parallel with the develop- 



ment of precision controlled single dopant devices ( Fuech 



sle et al. 2012) 



These experiments all relied on the in-diffusion of 
dopant atoms into nano-scale transistors from the source 
and drain regions. However several groups have concen- 
trated on adapting the industry standard technique of 
ion implantation to implant single dopants into a silicon 
device in a controlled manner ( Batra et al.\ 



et al. 



2008) 



2007 



Biele- 



jec et al. 2010 Jamieson et al.\ 2005 Schenkel et al. 



2003t |Seamons et al.\ |2008[ |Shinada et al.\ |2008[ |Weis| 



There are three main limitations to apply- 
ing the technique to scalable single atom architectures: 
the ability to register individual ion strikes, the over- 
all spatial registration of the ion implant site and subse- 
quent straggling of the ion due to the statistical nature 
of the stopping process. Once the dopants are implanted 
a rapid thermal anneal at high temperature is needed 
to repair the damage and activate the donors, causing 
diffusion and segregation of dopants (Park et al., 2004). 



Despite these concerns silicon nanoscale transistors have 
been fabricated with a low density of local area implants 
where it has been possible to tune individual dopants into 
resonance and observe transport spectroscopy through a 
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FIG. 45 (Color online) Transport through dopants ion- 
implanted in a nanoFET. (a) Schematic of a nanoFET 
where roughly 3 donors have been implanted into the 50 x 30 
nm active area of the device, (b) The stability diagram show- 
ing the differential conductance as a function of the barrier 
gate and dc source-drain bias, highlighting the resonant tun- 



neling peaks a%, b\ and ci of the three donors. From Tan 
[eTallpoTol 



metal-oxide-semiconductor ficld-effect-transistor ( John- 



son et al. 2010 Prati et al. 2012b). In the work by 



Johnson et al. (2010) FinFET devices were fabricated us- 
ing SOI with 20nm of Si on a 145nm thick buried oxide, 
giving nominal channel dimensions of 25 x 70 x 20 nm. 
Using 14 keV there was a 57% chance of a P + ion stop- 
ping within the channel region. The implant resulted 
in an increase in charge in the buried oxide, causing a 
shift in threshold voltage and an increase in series resis- 
tance consistent with the introduction of Frenkel pairs 
in the channel. In the more recent work of lPrati et al.\ 



(2012b) donors were placed in a ID array allowing the 
regime between single electron tunneling and Hubbard 
band formation due to inter-dot coupling to be investi- 
gated. 

There are 3 different regimes for single dopant trans- 
port experiments, as shown in Fig. |47| In the first the 
channel current is influenced by the presence of a neu- 
tral or charged dopant. Ono et al. (2007) identified a 



single acceptor that modified the current through a FET 
depending on it's charge state and refined this technique 



with a dual gate device to realize acceptor mapping (Kha- 



lafalla et al. , 2009 Ono et al. 2008 1 . In the second regime 



there is direct transport through a dopant in the access 
region to a FET channel ( Calvet et al. 2007a|b Hofheinz 



et al.. 2006a). This allowed the study of electric field 



single dopant, such as a P donor in Fig. 45 (Tan et al. and strain effects on the acceptor (Calvet et al. . 2007a|b ) 



2010). More recently single dopant implantation has 



been demonstrated into the channel of a silicon nanoscale 



and the Zeeman splitting ( Hofheinz et al. 2006a ) of the 
donor ground state. The third regime, representing the 
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most direct way to access information about the proper- 
ties of a dopant and its environment, is direct transport 
through the dopant states in the sub-threshold extreme 
of a transistor. Due to the progress in device scaling this 
was recently realized in ultra-scaled MOSFET devices by 



( Pierre et al. 


2010 


Prati et al. 


2009 


Sellier et al. 2006 


Tan et al. 


2010 


. All three transport regimes are illus- 
and will be discussed in more detail in 


trated in Fig 
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the remaining part of the sub-section, 
a) 



upper gate (V UG ) 




lower gate (V La ) Si substrate (V BG ) 



b) 




conduction band edge , 



C) 




FIG. 46 (Color online) Three examples of device layouts 
that illustrate different transport regimes for the de- 
tection of a single dopant, (a) Capacitive coupling to the 
channel which leads to a modification of the channel current 
due to the charge state of a dopant, (b) Tunneling through a 
dopant in the access region in series with transport through 
the channel, (c) Direct tunneling through a dopant in the 
channel in the sub-threshold regime. From (a^ 
[20071 (b) |Hofheinz et aZ.|[2006al (c) |Sellier et~al 



Ono et al. 



2006 



Sellier et al. ( 2006 1 demonstrated transport through a 
donor in a nano-FET and identified excited states as well 
as the doubly occupied D~ state confirmed by magneto 
transport (regime 3). The ionization energy observed in 
this work was consistent with As as expected from indif- 



fusion from the source and drain region. In contrast to 
the ionization energy of the D° state, which was similar 
to bulk, the ionization energy of the D~ state was en- 
hanced. This reduction in the charging energy, i.e. an 
increase in D~ ionization energy at an unchanged D° 
ionization energy, is discussed in comparison to several 



experiments in Sec. V.C.2 Lansbergen et al. (2008) ana- 



lyzed the orbital spectrum based on nonlinear transport 
and studied the impact of the electric field and gate inter- 
face on the donor state. They analyzed the data in com- 
parison to a large scale atomistic model and found good 
agreement even though the spectrum is strongly altered 
from the bulk. Tan et al. (20101 demonstrated transport 



spectroscopy on a transistor that was implanted with 3 
donors based on a timed exposure, thus the first exper- 
iment on an intentionally placed group of donors. The 
intentional placement is confirmed by a sample that was 
not implanted and did not show any resonances. As ex- 
pected the donor peaks show a paramagnetic Zeeman 
shift of the D° state and the opposite for the D~ state. 
Spectroscopic measurements combined with device based 
modeling led to a detailed understanding of the valley- 



orbit coupling of donor or dot states (Rahman et al. 



2011c) and the two electron state of a donor (Lansber- 



gen et al. 2011 Rahman et al. 2011a) which are dis- 
cussed in detail in the next section. Recently, |Pierre| 



et al. (2010) linked low temperature resonant transport 
through donor states to the room temperature perfor- 
mance of a MOSFET. They performed a statistical anal- 
ysis of the threshold voltage of ultra-short channel Fin- 
FETs (less than 20 nm). Transistors with a threshold 
voltage far below the average display resonant transport 
at low temperature due to a donor in the middle of the 
channel, as shown in Fig. |48[ This direct link between 
the room temperature variability and the low tempera- 
ture spectroscopy represents an important contribution 



to device engineering (Wacquez et al. 20101. 



Hofheinz et al. (2006a 



reported transport through a 
donor in the access region of a small MOSFET (regime 
2). The access region is part of the channel that has a 
weaker coupling to the gate which leads to a barrier in 
the band structure, as shown in Fig.[46p. They observed 
sequential transport through the donor and through a lo- 
calized state between the two barriers of the access region 
which has a much lower charging energy. Magnetotrans- 
port revealed a Zeeman shift of the resonance consistent 
with a paramagnetic trap. The complex interaction be- 
tween a single electron transistor and a dopant was only 
recently explained in detail ( |Golovach et al. 2011). Cal- 
vet et al. investigated acceptors in the barrier of a Schot- 
tky FET. They observed a Zeeman shift of the acceptor 
and analyzed the impact of the electric field on the ac- 



(Calvet et al. 2007a) 



ceptor ( Calvet et al. 2007b ) and the effect of local strain 



As mentioned before, single dopants can not only be 
detected by passing current through them but also by 
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FIG. 47 (Color online) Three-dopant transport regimes 
in a transistor geometry, (a) An example of the dopant 
detection regime based on the capacitive coupling the the 
channel for an undoped (left) and doped (right) double gate 
sample. The signature of a single acceptor charging event is 
evident in the doped sample. From Ono et al. (2007 1, (b) 
An example of the second regime where the dopant is in the 
barrier of the access region in series with a quantum dot, re- 
produced from |Hofheinz et al.\ |2006a| with kind permission 
of The European Physical Journal. Top line represents the 
room temperature FET characteristics and the line below the 
low-temperature Coulomb peaks, (c) The third regime with 
direct transport through a dopant in the sub-threshold limit, 
taken from ( jSellier et aH|2006[ ). 



their electrostatic coupling to free electrons in the chan- 



nel (regime 1). Khalafalla et al. (20091 developed earlier 



single gate experiments (Ono et al. 2007) to a multi- 
gate configuration. Two overlapping top-gates allowed 
the modulation of the channel potential profile. This 
made it possible to define the position of a donor along 
the length of the channel which led to the observation 
that an acceptor close to the source has the largest im- 
pact on transport. This device geometry has been further 
developed to allow for charge pumping based on a small 
number of dopants which is attractive due to the large 
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FIG. 48 (Color online) Direct tunneling through a 
dopant in a short-channel FET. (a) Illustration of a 
Monte Carlo simulation of the doping profile in a 20nm chan- 
nel where some dopants diffused into the channel region from 
the source and drain, (b) The dashed curve shows the cur- 
rent averaged over many devices where the black line indicates 
the threshold. Two devices show a drastically lower threshold 
linked to resonant transport at low temperature as indicated 
in (c) for the device with the lowest Vth- From Pierre et al. 
120101 This data shows the clear connection between the low 
threshold of these devices at room temperature and the reso- 
nant transport at low temperature, both mediated by a single 
dopant. 



charging energy of the dopant and its potential for high 



accuracy as current standard (Lansbergen et al. 2012 >. 
A double gate study of n-type channels revealed a clear 
difference between doped and undoped devices. In un- 
doped devices the back gate simply shifts the threshold 
voltage of the device. Devices of the same geometry with 
2 x 10 18 P cm~ 3 show many sub-threshold resonances with 
a large charging energy and a different coupling to the 
top gate. Furthermore, coupling between some of these 
is observed. These observations close the loop between 
the randomly- and deterministically-doped devices. 

Until recently single donor transport was focused on se- 
quential tunneling to probe the spectrum of a donor. Rc- 
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cently, several coherent transport experiments with sin- 
gle donor and double donor systems have been reported. 
Strong coupling to the contacts leads to the formation of 
a Kondo ground state where the parametric donor state is 
screened leading to transport in the traditional Coulomb 



blocked one electron region (Lansbergen et al. 2010). 



The valleys play a key role in Si Kondo physics as theo- 



retically predicted (Shiau and Joynt 2007 Shiau et al 



2007) and experimentally confirmed for donors by the 



presence of spin and orbital Kondo and their combination 
(Lansbergen et al. 2010 Tettamanzi et al. 2012). Fur- 



thermore, quantum interference between two tunneling 
paths has been discussed in the context of a Fano reso- 



nance (Calvet et al. 2011 Verduijn et al, 2010). Dopants 



present an ideal platform to investigate the crossover 
from scale-less mesoscopic effects to atomistic transport. 



3. Single dopants in crystalline silicon 

An alternative technology has been developed to place 
dopants in silicon with atomic precision using a scanning 
probe microscope. Ever since its invention in 1981 by 
Binnig and Rohrer, the scanning tunnelling microscope 
(STM) has gained international recognition by not only 
its capability to image surfaces with unprecedented res- 
olution but also by its potential to modify and pattern 
crystalline surfaces at the atomic scale. Amongst the 
most notable examples are the formation of the letters 
IBM with individual xenon atoms on a nickel surface by 



Eigler's group (Eigler and Schweizer 19901 



The translation of this technology to manipulate atoms 
in silicon was, however, not simple due to the strong, co- 
valent nature of silicon bonds. To position atoms in sili- 
con it was necessary to use a lithographic process, analo- 
gous to conventional optical/e-beam lithography. Here a 
monolayer of hydrogen resist is patterned using the tip of 
the STM to create an atomic-scale template. In the early 



1990s Lyo and Avouris (19901 and Lyding et al (1994) 



proposed the use of such a template to create a pattern of 
highly reactive dangling bonds sites on the silicon surface 
which could subsequently be functionalized with various 
atomic and molecular species. 

The process of STM hydrogen lithography has since 
been adapted to realize a complete fabrication strategy 



for atomic-scale silicon device fabrication (Ruess et al 



2005 2004 Simmons et al 2005). Here dopants are 



placed in the silicon crystal with atomic precision lat- 
erally using scanning probe techniques and atomic preci- 
sion vertically with molecular beam epitaxy (MBE) . The 
important feature of these devices is that, in contrast to 
modern CMOS devices and almost all quantum semicon- 
ductor devices, which use many materials and have het- 
erogeneous interfaces, STM-patterned devices are formed 
in single crystal silicon. Confinement of electrons is thus 
achieved by atomically abrupt changes in the density of 



dopant atoms within the silicon crystal. The doped re- 
gions can have very high planar electron densities in the 
range n2D ~ 2.5 x 10 14 cm" 2 . In three dimensions this 
density corresponds to a value ps 10 21 cm -2 , three orders 
of magnitude above the Mott metal-insulator transition. 
At these high carrier densities one in every four silicon 
atoms is substituted with a phosphorus atom, so that the 
average separation of phosphorus atoms is < 1 nm, much 
smaller than the single dopant Bohr radii (ag ps 2.5 nm). 
As a consequence the doped regions are highly conduct- 
ing and behave like a disordered metal. Electron trans- 
port has been studied in these highly doped single crystal 
donor-based quantum wires ( |Ruess et al] [2008} [2007a|b| 
and recently this technique has been used to realize con- 
ducting wires in silicon with Ohmic behavior down to the 



atomic-scale (Weber et al. 2012a) 



It is also possible to pattern more complex devices, 
such atomically abrupt, epitaxial quantum dots with ID 



source and drain leads (Fuechsle et al. 2010 Fuhrer 



et al. 20091. These studies have allowed the impact 



of vertical and lateral confinement on silicon quantum 
dots to be investigated. In such abruptly confined quan- 
tum dots very small energy level splittings of ps lOOfxeV 
have been observed in electron transport and attributed 
to transport through the valley states of a few-electron 



quantum dot, as shown in Fig. 49 (Fuechsle et al 2010 1. 



STM images of the device in Fig.|49[a,b) show the central 
region of the device into which 6±3 P atoms are laterally 
confined with the STM. These dopants are also strongly 
vertically confined by low temperature silicon molecu- 
lar beam epitaxy. The corresponding stability diagram 
shown in Fig. |49[ c,d) reveals spacing in the energy spec- 
trum of order ps 100/ieV. This very small energy level 
splitting was surprising given the ultra-small size (~ 4 
nm 2 ) of the quantum dot. However it is well known that 
strong lateral and vertical confinement breaks the degen- 
eracy of silicon valley states. 

This is illustrated in Fig.|49{e) where electrons strongly 
confined in a two-dimensional plane result in splitting of 
the 6 bulk valleys into 4 degenerate A-pockets as well as 
two T-pockets at k = 0. The remaining degeneracies can 
be broken in the presence of sharp lateral or vertical con- 
finement. In these highly doped 5-layers strong, abrupt 
quantum confinement in z-direction splits the degener- 
acy of the out-of-plane T-bands to give the lower energy 
Ti and T2 bands. The four A-valleys, two each in the 
x and y directions, are usually degenerate in 2D devices. 
However, since the quantum dot device is also confined 
laterally on the nm-scale valley splitting of these states 
is also observed. It is this valley splitting that gives rise 
to the ss 100/ieV energy level separation observed exper- 
imentally in STM-patterned few-electron quantum dots. 
This contrasts to studies of other few-electron quantum 
dot systems, where it is the large size of the dot itself 
that gives rise to m lOO^ieV (IBeenakkerl 11991b. 



The source and drain leads to these quantum dots are 
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decreasing symmetry 



FIG. 49 (Color online) Few-electron quantum dot. (a) 

An STM image of the central device region of a few-electron 
single-crystal quantum dot acquired during hydrogen lithog- 
raphy, showing a four terminal device with source (S), drain 
(D) and two in-plane gates (G1,G2). The bright regions corre- 
spond to areas where phosphorus donors will be incorporated, 
(b) A close-up showing the central quantum dot containing 
6 ± 3 donors, (c) Stability diagram showing the conductance, 
dl /dVsD through the dot as a function of gate voltage, Vg 
and bias voltage Vsd- (d) A close-up of the transition (white 
square in (c)) reveals a high density of conduction resonances 
with an average energy spacing of « 100/^eV\ (e) The 6-fold 
degeneracy of the conduction band minima of bulk silicon is 
lifted by confining the electrons vertically to two dimensions 
and is then split again by abrupt, lateral confinement. From 
|Fuechsle et~aL\ [20TO| 



not 2D reservoirs but ID leads. This results in the forma- 
tion of resonant tunneling features due to the presence of 
ID subbands. The energy separation of these subbands 
depends on the width of the ID leads and for this de- 
vice was found to be ~10 meV. Such ID states have also 



been observed in other silicon-based quantum dots (Lim 



et al. 2009b Mottonen et al. 2010) and a recent review 



discusses how to distinguish these in electron transport 



(Escott et al. 2010) 



Using such a technology it is also possible to pat- 
tern individual dopants in silicon with atomic preci- 



sion(Schofield et al. 2003 Wilson et al. 2004) to real 
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FIG. 50 (Color online) A single-atom transistor, (a) 

3D perspective STM image of a hydrogenated silicon surface. 
Phosphorus will incorporate in the red shaded regions selec- 
tively desorbed with an STM tip to form electrical leads to 
a single phosphorus atom patterned precisely in the centre, 
(b) the source (S), drain (D) and two gate leads (Gl, G2) 
to the central donor, which is incorporated into the dotted 
square region, (c) The electronic spectrum of the single-atom 
transistor, showing the drain current Isd as a function of 
source-drain bias Vsd and gate voltage Vg applied to both 
gates, (d) The differential conductance dlsD /dVsD as a func- 
tion of Vsd and Vg in the region of the D° diamond shown 
in (c). (e) A comparison of the potential profile between the 
source and drain electrodes in this device (straight line) to an 
isolated bulk phosphorus donor (dashed line), where the D° 
state resides 45.6 meV below E c \>. In contrast, the D° state 
in the single-atom transistor resides closer to the top of the 
potential barrier. From Fue chsle et al.\ |2012| 



are that the dopant, the source-drain leads and the con- 
trol gates are crystalline and all exist within one plane of 

The encapsu- 



ize single dopant atom transistors, as shown in Fig.[50[a) 



(Fuechsle et al. 2012 ). The critical features of this device 



the silicon crystal, as shown in Fig. 50 'b). 
lation of this device in epitaxial silicon removes the con- 
fined dopant states away from the influence of surfaces 
and interfaces. However transport devices by definition 
also contain electrodes, and these electrodes are known 
to have profound effects on the energetics of the single 
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dopant atom. In these epitaxial architectures the elec- 
trostatic potential at the dopant could be tuned using 
two in-plane gates Gl and G2 patterned either side of 
the transport channel defined by the S and D leads. 

'c) presents the measured stability diagram of 



Fig. 50 



the single donor, in which the three charge states of the 
donor can easily be identified: the ionized D + state, the 
neutral D° state and the negatively charged D~ state. 
The diamond below Vg ~ 0A5V does not close, as ex- 
pected for the ionized D + state, because a donor can- 
not lose more than its one valence electrons. The con- 
ductance remains high (on the order of microsiemens) 
down to the lower end of the gating range, making the 
possibility of additional charge transitions unlikely. The 
D + — > D° charge transition occurred reproducibly at 
V G ~ 0.45 ± 0.03F, as shown in FigQd) for multiple 
cool-downs and is attributed to the high stability of the 
device and the inherent influence of the nearby electrodes 
on the position of the donor eigenstates relative to the 
Fermi level of the leads. 

To understand quantitatively how the nearby trans- 
port electrodes affected the electronic properties of the 
donor, the electrostatic potential landscape of the in- 
nermost part of the device was calculated, treating the 
heavily doped gate regions in a self-consistent atomistic 
approach using a Thomas-Fermi approximation. Hav- 
ing established the electrostatic potential of the device, 
the donor electronic states were then calculated using 
a tight-binding approach (Lansbergen et al. 2008). The 



position of the resulting one-electron ground state D° for 
the solitary phosphorus dopant is depicted in Fig.[50{e) 
(blue line). As expected, due to the electrostatic envi- 
ronment, the energy levels of the device are raised sig- 
nificantly from the bulk case (dashed grey line), where 
the unperturbed Coulombic donor potential asymptoti- 
cally approaches the silicon conduction band minimum 
E c i, (red dashed line) and D° has a binding energy of 
Eb ~ —45.6 meV. In contrast, D° in the effective donor 
potential of the single donor transport device resides 
much closer to the top of the barrier (solid line) along 
the S-D transport direction. Despite this, the charging 
energy Ec could be extracted from the transport data 
and was found to be 47 ± 3 meV, remarkably similar to 
the value expected for isolated phosphorus donors in bulk 
silicon (« 44 meV)(Ramdas and Rodriguez 1981). 



These results are in sharp contrast to previous ex- 
periments on single dopant in silicon transport devices, 
which have reported charging energies that significantly 



differ from the bulk case ( 


Lansbergen et al. 2008 


Pierre 


et al. 


2010 


Rahman et al. 


2011a 


). There, the difference 



was attributed either to screening effects resulting from 



strong capacitive coupling to a nearby gate (Lansber- 



gen et al. 2008 ) or strong electric fields ( Rahman et al. 



2011a), or to an enhanced donor ionization energy in the 



proximity of a dielectric interface (Pierre et al. 2010). 



rus dopant, which is symmetrically positioned between 
two gates, encapsulated deep within an epitaxial silicon 
environment. 



C. Discussion 

A dopant in a semiconductor represents the ultimate 
limit of a quantum dot. In contrast to a quantum dot the 
confinement potential is given by the three dimensional 
Coulomb potential of the dopant ion (see section [TT]) and 
not by external gates. In quantum dots the orbital en- 
ergy is small and mainly probed as the energy difference 
between the triplet and single sates in the two electron 
problem. The dopant has hydrogenic level spectrum with 
splittings in excess of lOmeV. The valleys in Si lead to 
interesting corrections due to the restricted momentum 
space of these states. Again in strong contrast to a quan- 
tum dot, the confinement potential of a dopant is strongly 
altered by the amount of charge on the dopant. A shal- 
low impurity can only bind two electrons where the 2nd 
charging energy for the 2nd electron is almost the ion- 
ization energy of the first, i.e. H~ like the two electron 
state is very close in energy to the continuum. Here we 
compare the physical properties of the dopant confined 
states for the different devices. The orbital spectrum of a 
gated donor is illustrated followed by a discussion of the 
charging energy between the one and two electron state. 
Finally, we look at the interaction between the donor and 
the leads and the interaction between two donors. 



1. Orbital structure of a dopant in a nanostructure 



Section [Til . C | discussed the spectrum and orbital struc- 
ture of a bulk dopant. A dopant in a nanostructure pos- 
sesses not necessarily these properties as already calcu- 



lated by Macmillen and Landman (1984). The environ 



ment in a nanodevice has a large impact on the orbital 
spectrum of a donor. An electric field will lift degen- 



eracies (Friesen 2005 Smit et al. 2004) and a triangu- 



lar well, e.g. due to a gate, lowers the excited states of 
the dopant due to the interaction with the interface well 
(iCalderon et all 120071 iMartins et al] 120041). These the- 



oretical predictions have been experimentally confirmed 



in the interface (Lansbergen et al. 2008) as well as the 



bulk regime (Fuechsle et al. 2012 1. Critically important 



Importantly, these effects are small for a single phospho- 



for the ability to model the devices and obtain metrol- 
ogy data with respect to impurity type and depth, was 
the ability to compute the excited state spectra as well 
as the ground state spectra in NEM03D, see section III 
of the supplemental material. The effective mass models 
offer key physical insight into the problem and go hand 
in hand with the tight-binding work which generates ac- 
curate predictions to interpret the experimental data. 
The environment, i.e. the leads, an interface to a gate, 
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or an electrical field, can drastically alter the orbitals of 
a dopant in comparison to the unperturbed bulk con- 
dition. Lansbergen et al. (2008) measured the energy 



spectrum of single donors, located in the channel of Fin- 



FETs by transport spectroscopy, as shown in Fig. 51 
They were not bulk like but agreed well with a multimil- 
lion atom simulations of the complete system. In con- 
junction with the data, the theoretical analysis allowed 
the authors to identify the species of the donors (As) and 
furthermore provided an explicit determination of the de- 
gree of gate-controlled quantum confinement in each de- 
vice. Figure [5l| b) shows the 3 confinement regimes that 
can be distinguished: Coulomb, hybridized, and interfa- 
cial confinement with the charge density and schematic 
potential landscape for these three regimes. At low elec- 
tric fields the electron is located at the donor site and 
its ground state corresponds to a donor in bulk (thus full 
lattice symmetry). At high electric fields the electron is 
pulled inside the triangular potential well at the interface 
reducing the symmetry of the system. The electron is 
still localized near the donor site in the lateral directions 
though, in correspondence with the results of |Calder6n| 
et al. (20071. At the crossover between these regimes, the 
electron is delocalized over the donor- and well-sites. 

Donors in devices fabricated with the STM, as dis- 
cussed in sub-section V.B.3[ exh ibit a more bulk like or- 
bital spectrum (Fuechsle et al. 2012). This is due to 



the fact that in these devices the donor is far away from 
a dielectric interface and the gate does not create large 



fields in comparison with MOSFET structures (Lansber- 



gen et al. 2008). Fuechsle et al. (20121 also showed that 



the charging energy of a dopant in an STM fabricated 
device is comparable to the bulk value. This is to be ex- 
pected for an environment that is close to bulk besides 
the source and drain electrodes with a cross-section of 
only a few square nanometers. 



2. Charging energy of a dopant in a nanostructure 




Coulomb confinement Hybridized regime Interracial confinement 
regime regime 



FIG. 51 (Color online) Excited-state spectroscopy of 
single gated donors, (a) Differential conductance of a 
dopant in a FET. Excited states are indicated by the dots 
and arrows. Inset in (a) shows current 7sd as a function of 
gate voltage at Vb = 40 mV where each plateau indicates ad- 
dition of a quantum channel due to an orbital, (b) Simulations 
of the gated donors eigenstates: wave function density of the 
D°-ground state (| ^gs | 2 ) located 4.3 nm below the interface 
in three different electric field regimes: Coulomb confinement 
regime, OMVm" 1 (left), hybridized regime, 20 MV m 1 (mid- 
dle) and interfacial confinement regime, 40MVm _1 ( right). 
The grey plane indicates the Si/SiCh interface. From Lans- 
|bergen et al] [20081 



In the constant-interaction model (Beenakker 1991) 
the charging energy of a Coulomb island is independent 
of the number of electrons ./V localized on the charge is- 
land. This assumption is valid as long as the confinement 
potential is not affected by N which is not at all the case 
for isolated donors. For shallow donors, only a single 
charge transition (N = 1 — > N — 2) plays a role since it 
is not possible to bind a 3rd electron. The addition of 
a single electron to an ionized donor site will screen the 
positive nucleus and thus strongly alter the confinement 
potential for an additional electron. The Coulomb inter- 
action between an electron on the donor and all other 
electrons in its environment can still be parameterized 
by a single capacitance C, which is specific to N and the 
donor environment. The charging energy, represented by 
e 2 /2C, of donors close to a gate interface is modified due 



to the screening at the interface as well as the applied 
electric field which was experimentally demonstrated by 
Lansbergen et al. (2008). Fuechsle et al. (2012) showed 



that an STM fabricated single dopant device displays 
an unaltered charging energy consistent with the bulk- 
like environment. Recent theoretical work addresses this 
problem and progress has been made using effective mass 



( Calderon et al. 


2010b 


Fang et al. 


2002 


Hao et al. 


2011 


Hollenberg et al. 


2004) and self-consistent field tight- 



binding (Rahman et al. 2011a) treatments 
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3. Interactions between donors 

The interaction between donors plays a central role in 
quantum information science. The goal is to achieve tun- 



able interaction that preserves coherence (Kane 1998). 



This has not been achieved yet but it is within reach. 
Initial experiments have focused on the study of capac- 
itive coupling as well as the coherent coupling between 
dopants. A detailed understanding of tunnel coupling as 
well as capacitive coupling between a donor and a SET is 
a key issue since this is the central read-out mechanism 



for qubits ( Morcllo et al. 2009 1 . This complex coupling 



between a dopant and a quantum dot, i.e. a semiconduc- 
tor SET, has been analyzed in detail in an experimental 
and theoretical effort by Golovach et al. (2011). Coher- 



ent coupling between dopants has been achieved in the 



limit of weakly coupled dopants (Verduijn et al., 2010) 



and strongly coupled dopants (Calvet et al. 2011). Both 



rely on the interference between two coherent transport 



that is sensitive to the phase difference between the two 
transport paths. This phase difference can be modified 
by changing the magnetic flux that is enclosed in the loop 



of the transport paths. Verduijn et al. (20101 studied two 



As atoms in a nano MOSFET and showed that the dis- 
tance between the dopants is about 30 nm based on the 



magnetic field dependence. Calvet et al. (2011) studied 



acceptors in a Schottky FET and also observed a Fano 
resonance which proofs coherent exchange of electrons. 
The lack of magnetic field dependence confirms their ex- 
pectation that the acceptors studied are strongly coupled 
in these devices. 



D. Double dopant quantum dots 
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channels which leads to a specific line shape (Fano 1961 ) Fig- |53 
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FIG. 52 (Color online) Sequential transport through 
a double donor device with independent gate con- 
trol. The left panel shows the two opposing gates similar to 
a conventional FinFET geometry but with a split gate. The 
channel received a background doping of 10 18 P/cm 3 and this 
device demonstrates independent gate control of two dopants. 
The right panel shows a finite bias stability diagram reveal- 
ing bulk-like excited states of the dopant. From Roch e et al.\ 
I20T21 



Few-electron single-crystal quantum dots have recently 
been realized using STM-patterned devices, as shown in 
Here independent electrostatic control of the 



ultra-small dots was achieved by careful modeling and 
optimization using the capacitance modeling tool FAST- 
CAP ( |Nabors and White[ [1991]) and a single-electronics 
modeling tool (SIMON) (|Wasshuber et aL[|1997[). This is 



quite remarkable given the small physical size of the dots 
(~4nm in diameter) and their close spacing (~10 nm). 
At such small dimensions cross-capacitances between the 
quantum dots become considerable and the dots need 
to be positioned at an angle, a ~ 60° with respect to 
one another to achieve independent electrostatic control. 
One of the advantages with donor-based quantum dots 
is that this combination of device modeling and preci- 
sion lithography using scanning probe microscopy allows 
reliable predictive device design, an important tool as 
devices scale to the single donor level. 



The study of transport and interactions in donor-based 
double quantum dots has been motivated by their poten- 



tial for solid state quantum computing applications ( Loss 



and DiVincenzo[|1998[|Taylor et al.\\20Qh\ . Initial stud- 
ies concentrated on ion-implanted devices, where both 



independent gate control (Hudson et al. 2008) on the 



dot occupancies and charge detection using surface alu- 



minum SETs were demonstrated (Mitic et al. 2008) in 



multi-donor devices which contained hundreds of dopants 
in each dot. Characteristic honeycomb structures (see 
section IV. F ) were observed in the charge stability maps. 



However difficulty was encountered going to smaller dot 
sizes due to the inherent straggling in the ion implanta- 
tion process. Non the less, sequential transport through a 
stochastically doped FinFET structure has been demon- 
strated by (Roche et al. 2012). They used a split gate 



geometry to independently control the chemical potential 
of two dopants and probe the excited states by tunneling 
spectroscopy in a similar manner to a double quantum 



dot, as shown in Fig. 52 



E. Charge sensing in few-electron dopants 



Following the principal work by Field et al. ( 1993 ) on 



remote sensing of charge using a quantum point con- 



tact in a AlGaAs/GaAs system, Elzerman et al. (2004) 
adapted this technique to perform single shot detection 
of spin dependent single electron tunneling events in a 
single electron quantum dot. Real time sensing of single 
electron tunneling is fundamental to electrical read-out of 
qubit states in spin quantum computing. In these exper- 



iments ( Elzerman et al. 2004 1 a quantum point contact 
was capacitively coupled to the qubit. However the vis- 
ibility, and therefore fidelity of spin read-out, of these 
charge detectors is greatly enhanced when a single elec- 
tron transistor rather then quantum point contact is used 



as a charge detector ( Morello et al. 2009 1 . Here the SET 



is additionally tunnel-coupled to the qubit and electrons 
can be loaded from the SET-island itself, thus eliminat- 
ing the need for a separate electron reservoir. 
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FIG. 53 (Color online) A donor-based double quantum 

dot in silicon, (a) An overview STM-image of the device 
showing the two quantum dots, tunnel-coupled to the source 
and drain (S/D) leads and capacitively coupled to the gates 
Gi(2)- (b) Close-up of the two quantum dots, ~4 nm in diam- 
eter. The DQD angle a— 60 ± 3° has been optimized for max- 
imum electrostatic control whilst suppressing parallel leakage 
through the dots. The modeled (c) and measured (d) charge 
stability diagrams show excellent agreement, demonstrating 
independent electrostatic control of the individual dots. From 



Weber et al. (2012b) 



High-fidelity spin read-out of a P-donor-bound electron 
in Si has recently been demonstrated within this archi- 



tecture (Morello et al. 2010) by implanting a small num- 



ber of P donors in the vicinity of an electrostatically in- 
duced SET, at the Si-Si02 interface. This has established 
the feasibility of fiducial detection of P-donor-based spins 
qubits, however the uncertainty in the number and posi- 
tion of the donors relative to the SET is of concern for 
scaling up to a practical 2D Si:P quantum computer with 
multiple donor arrays. To reproducibly achieve sufficient 
charge sensitivity and electron tunnel rates, it is impor- 
tant to precisely situate the SET with respect to deter- 
ministically positioned array of Si:P qubits. Towards this 
end STM-lithography has recently developed a charge de- 
tection device lay-out, wherein a SET charge sensor and 
an ultra-small quantum dot are both patterned on the 
same plane of a Si crystal, at an atomically-precise sepa- 
ration, as shown in Fig. [54] 

In this way two of the key design parameters, i.e., the 
sensitivity of charge detection and electron tunnel rates 
can be made sufficiently high to enable projective spin 
readout of individual P donors in Si. Since the fabrication 
technique is essentially identical to the established ap- 
proach for deterministic placement of P donors in Si, this 
charge sensing layout can be readily integrated in scalable 
Si:P spin quantum computing architectures. Recent re- 
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FIG. 54 (Color online) Charge sensing using a donor- 
based single electron transistor coupled to a small 
donor dot. (a) Filled-state STM image of the overall device 
pattern, showing (in lighter contrast) the regions where the 
hydrogen resist monolayer has been desorbed to create the 
source (S) and drain (D) contacts of the single-electron tran- 
sistor, and the two gates (Gi, G2). (b) High-resolution image 
of the device pattern within the white box in (a), showing the 
SET-island (Di) and the quantum dot (D2) (c) Charge stabil- 
ity plot showing the dependence of Isd on the gate voltages 
(Vgi, Vg2), for a constant Vsd = -50 fiV. The high current 
lines correspond to the Coulomb peaks of the SET. Inset: 
High resolution map of a small section of (c) showing discon- 
tinuity of a current line, due to a particular charge transition 
of D2. The triangles in the main map indicate a total of 7 
such transitions of D2. From Mahapa tra et al.\ |2014| 



suits demonstrate that single electron tunneling between 
the quantum dot and the SET island occurred on a time- 
scale (r ~ ms) two-orders-of-magnitude faster than the 
spin-lattice relaxation time of a P donor in Si and suit- 



able for projective readout of Si:P spin qubits (Mahapa- 



tra et al. 2011 1. Another work by Mazzeo et al. (2012) 



reported on the charge dynamics of a single donor cou- 
pled to a few-electron silicon quantum dot. Single-spin 
sensitivity is discussed in more detail in section VI. A. 



VI. RELAXATION, COHERENCE AND 
MEASUREMENTS 

In the previous sections we have reviewed the quantum 
electronic properties of silicon, and how such properties 
can be harnessed in nanoscale structures. That field of 
research is fairly mature, partly thanks to its technolog- 
ical links with classical silicon nanoelectronics. 

Once the confinement of a single electron in silicon has 



45 



been achieved, it is possible to start exploring the proper- 
ties of the electron spin. An exciting possible application 



relies on well-established techniques that belong to the 



vast field of magnetic resonance (Slichter 1990). Even in 



is quantum computing (Ladd et al. 20101, where the 



electron spin is used as a quantum bit (Morton et al. 



2011). Another vast field of research is spin transport 
^ Appelbaum et al. 2007[ Dash et al. 2009) and spin- 
tronics ( Awschalom and Flattel 2007 Jansen 2012 Zutic 



et al. 2004), which we will not discuss in this review. 



In this section we present a theoretical introduction 
to the behavior of electron spins confined in silicon, and 
a snapshot of the current state of experimental progress. 
The discussion below is meant to provide the reader with 
an outlook on the future direction of the field. A thor- 
ough review of spin control in silicon will only be possible 
several years from now. 



A. Spin relaxation and decoherence 

The suitability of a physical system to encode and pre- 
serve quantum information is quantified by parameters 
such as the relaxation and coherence times. The first, 
called T\, describes the timescale over which an energy- 
excited state decays to the ground state. The second, 
called T2, describes the timescale over which the phase 
coherence between different branches of a quantum su- 
perposition can be preserved. 

A spin S — 1/2, such as an electron spin confined in 
a Si structure by a natural or artificial potential, can be 
described by a 2 x 2 density matrix p = \^){tjj\. The state 
vector belongs to a 2-dimensional Hilbert space with 
basis vectors |!), normally chosen as the eigenstates 
of the spin operator S z if a static magnetic field B \\ z is 
applied. At thermal equilibrium, the diagonal elements 
of the density matrix are related by a Boltzmann fac- 
tor pfl = exp^Ez/ks^plz, where E z = gp B B is the 
Zeeman splitting. Therefore, the diagonal elements of p 
are related to the degree of spin polarization. After a 
perturbation involving exchange of energy with an exter- 
nal field or reservoir, the spin returns to equilibrium in 
a typical time scale T\, e.g. pu(t) — p*j cx exp(— t/T\). 
T^ 1 represents the spin-lattice relaxation rate. T% mea- 
surements in bulk samples are performed by observing 
the timescale over which the thermal equilibrium is re- 
covered after either an 'inversion pulse', which swaps the 
populations of the ground and excited spin states, or a 
'saturation comb', which equalizes the populations. Sec- 
tion |VI.C.3| discusses how to obtain the T\ of a single 
spin in a nanostructure from a measurement of the prob- 
ability of detecting the spin excited state as a function of 
the waiting time after the excited state preparation. 

A coherent superposition of the |t)>|-l) basis states 
results in nonzero off-diagonal elements ('coherences') 
P12 = 9%x- The preparation and manipulation of such 
coherent superpositions is at the heart of quantum in- 



the absence of energy exchange with the environment, the 
coherence may decay in time like p\-2.(t) oc exp[— (t/T 2 ) a ], 
where T^ 1 is the decoherence rate, and a is an exponent 
that depends on the details and the dynamics of the envi- 
ronment coupled to the spin. For electron spins in solid 
state, a major contribution to decoherence is given by 
the hyperfine coupling between the electron and the sur- 
rounding nuclear spins. The nuclear spins exhibit com- 
plex dynamics, driven by the interplay of their mutual 
interactions and the coupling with the electron. The time 
fluctuations of the hyperfine field randomize the electron 
spin precession frequency and destroy its coherence. 

Even in the presence of a perfectly static nuclear spin 
bath, a macroscopic ensemble of spins would exhibit a 
spread of precession frequencies as a consequence of inho- 
mogeneity in the local magnetic field, caused e.g. by the 
difference in the instantaneous value of the local hyper- 
fine field at every electron site. The resulting dephasing 
time T 2 * represents the timescale over which a free induc- 
tion decay occurs, i.e. the vector sum of all the spins in 
the ensemble averages to zero. For a single spin, the free 
precession cannot be observed in a single experiment, and 
must be obtained through repetition averaging. There- 
fore, a T 2 * -process arises when the quasi-static value of 
the local magnetic field changes from one repetition to 
the next. The 'true' decoherence time T2 is obtained in 
the experiments by 'refocusing' the quasi-static inhomo- 
geneity (in space, for a spin ensemble, or in time, for a 
single spin) of the magnetic field through a Hahn echo 



technique (Slichter 19901. 



The definitions and discussion above can be readapted 
to the case where a two-level system is obtained from 
the truncation of the Hilbert space of two exchange- 
coupled spins ( |Levy[ |2002| |Petta et al\ |2005| . The 
basis states then become the singlet and triplet states, 
\S) = (It!) - \lt))/V2 and \T ) = (\U) + llt))/A and 
the energy splitting caused by the exchange interaction J 
replaces Ez in the expressions above. This way to define 
spin-based two-level systems has been proposed to allow 
the control of the qubit purely by electrical means, i.e. 
without resorting to magnetic resonance techniques. 

Relaxation and decoherence of spins in semiconductors 
has been the subject of intense research, and an accessible 
review is given in (|Hanson et al. 2007 I . Here we highlight 



formation technology (Nielsen and Chuang 2000), and 



the specific phenomena that arise in silicon, in particular 
due to the valley degeneracy of the conduction band. 

The spin-lattice relaxation, i.e. the return of the diag- 
onal elements of the spin density matrix to their equilib- 
rium value, requires the coupling of the spin to a phonon 
reservoir. Silicon lacks piezoelectric effect, which is of- 
ten the dominant source of spin-phonon coupling in III- 
V materials. The only type of phonons present in Si is 
the "deformation potential" , i.e. a local change in lattice 
spacing which propagates with wave vector q. A defor- 
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mation potential phonon alters the band gap in an inho- 
mogeneous and time-dependent way, with repercussions 
on the exact mixture of spin, valley and orbital nature of 
the electronic wave functions. 

The relaxation rate T^ 1 is obtained in a 'Fermi golden 
rule' approach as: 

Tr 1 « \Ke- P h;So\l)\ 2 N(E z ) (12) 

where N(Ez) oc E\ is the density of phonon states at the 
energy splitting Ez, and T-L e -ph-so is a Hamiltonian term 
that includes the electron-phonon interaction and the 
spin-orbit coupling. It should be noted that the electron- 
phonon interaction does not directly couple Zeeman-split 
pure spin states. However, a nonzero coupling is obtained 
if the true eigenstates contain admixtures of other orbital 
or valley states, mixed in by the spin-orbit interaction. 

Another way to look at the problem - more familiar 
to the spin resonance community - is to think of the 
spin as being subject to an effective magnetic field, whose 
magnitude and direction can be modulated by a lattice 
phonon. Then T^ 1 is proportional to the spectral den- 
sity, at frequency w e = Ez / fi, of the component of the 
phonon-induced fluctuating local field perpendicular to 
the spin quantization axis. 

Let us recall a simple expression for the electron g- 



factor in a semiconductor (Kittel 1963 Roth 1960): 
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2Aso,vb 
3E g + 2Aso,vb 



(13) 



where Ago.VB is the spin-orbit splitting of the valence 
band, E g is the band gap, m and m* are the free elec- 
tron and the effective mass, respectively. In Si, Aso.vb ~ 
40 meV is relatively small, due to the small atomic num- 
ber. The large band gap E g — 1.12 eV results in electron 
g-factors very close to 2. Accordingly, spin relaxation in 
Si is relatively slow, since the modulation of the g-factor 
due to phonon scattering is very small. 

There are four main differences in the spin relaxation 
behavior between Si and III-V semiconductors such as 
GaAs ( |Blakemo"re 1982): (i) Si has no piezoelectric ef- 
fect, therefore only deformation potential phonons are 
present; (ii) Si has no bulk inversion asymmetry, there- 
fore is immune from Dresselhaus spin-orbit coupling ef- 



fects (Dresselhaus 1955 Hanson et al. 2007); (iii) The 



small atomic number and large band gap of Si produce 
a weak spin-orbit coupling; (iv) The physical mechanism 
and the magnetic field dependence of T^ 1 depends on the 
nature (valley or orbital) and the details of the excited 
states above the valley-orbit ground state. 



1. Electron spin relaxation in donors 

Let us consider first the case of a shallow donor such 
as P, As or Sb, where doublet and triplet valley-orbit ex- 
cited states (see section III.B.2 ) lie ~ 10—15 meV above 



the singlet (spin-degenerate) ground state. All of these 6 
states share the same hydrogenic Is orbital nature, and 
lie well below the 2p orbital states. In this situation, 
the dominant contribution to spin relaxation arises from 
valley effects, which can take two forms: (i) "valley re- 
population" or (ii) "one-valley" mechanisms. 

(i) "Valley repopulation" : within one valley, the g- 
factor of an electron is slightly different (anisotropic) for 
magnetic field parallel (gy) or perpendicular (<?j_) to the 
valley axis. In the unperturbed Is singlet ground state of 
a donor, all 6 valleys contribute equally, and the overall 
g-factor of the donor-bound electron is isotropic. How- 
ever, the local strain produced by a phonon has the ef- 
fect of disrupting the symmetry of the 6 valleys, lowering 
certain valleys with respect to others. Now the (/-factor 
may assume an anisotropic character, which can be inter- 
preted as a phonon-induced modulation of the effective 
local field, with a component perpendicular to the spin 
quantization axis. The resulting relaxation rate becomes 



(Hasegawa 1960): 



T 1 -\B,T) = f Sl ( 
1 
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= KaB^T. 



(14) 
(15) 



Here 2 is the deformation potential parameter ( Bardeen 



and Shockley 1950 Herring and Vogt 1956), represent 



ing the energy shift of the valleys due to a deformation 
of the crystal lattice, g' — (gi — <?t)/3 describes the 
anisotropy of the (/-factor along the principal axes of the 
effective mass tensor for each valley, E vo is the energy 
difference between the first excited valley-orbit state and 
the ground state, p = 2330 kg/m 3 is the density of Si, and 
v t = 5860 m/s and v\ — 8480 m/s are the transverse and 
longitudinal sound velocities, respectively. /si(^i ( t ) ) is an 
angular factor that goes to zero for 9 = (B || [001]), 
and is maximum for B II 111 11. 



Eq. 14 was derived in the high-T limit, appropriate 
for typical X-band (~ 10 GHz) ESR experiments at 
T > 1 K ( |Feher and "Gere] p59| . The full expres- 



sion contains the term (1 + n p h) ~ hsT/g^sS, where 
n>ph = (exp((7/xsB/fcsT) — 1) 1 is the Bose occupation 
factor of the phonon mode at the Zeeman energy. For this 
reason, T^ 1 oc T in the high-T limit. Conversely, single- 



spin experiments in nanostructures (Morello et al.\ 2010) 



are conducted in the low-T limit, where (1 -I- n p h) ~ 1- 
This indicates that only spontaneous emission of phonons 



can take place. Eq. (15) becomes 

9HbB 



T^(B)\ low _ T = K^-f^B 4 T = K 5 B 5 



(16) 



The Tj -1 oc B 5 dependence arises from the following fac- 
tors: (i) The density of phonon states N(Ez) is propor- 
tional to B 2 ; (ii) In the matrix element (f \He-ph;So\V) ■> 
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a factor proportional to B accounts for the need to break 
time-reversal symmetry, while another factor propor- 
tional to yB arises from the Jq dependence of the strain 
caused by a deformation potential phonon, where q oc 
gf^BB is the wave number. Therefore, |(f \H e -ph-so\l) | 2 
is proportional to B 3 . 

(ii) u One-valley" mechanism: A phonon-induced strain 
introduces a coupling between the T band and the nearest 



Donors 



Quantum dots 



A band (see section III.B.2 1 . This yields an anisotropic 
modulation of the <?-factor even for an electron confined 
to a single valley. This one-valley mechanism also yields 
T^ 1 oc B 5 in the low-T limit, but has a different angu- 
lar dependence, with fastest relaxation for B || [001] and 
slowest along [111]. Unlike the valley repopulation, the 
one-valley mechanism always produces a nonzero relax- 
ation rate. 

A detailed discussion and experimental study of these 
relaxation channels for Si:P was given by | Wilson and Fe^ 



her (1961), in the high-T limit (see Fig. 55 'a)). Since 
both the valley repopulation and the one-valley mecha- 
nism are generally active at the same time and have a 



comparable strength, Wilson and Feher ( 1961 ) included 
the analysis of the angular dependence of T{~ to unravel 
the different contributions. The low-T limit has been in- 
vestigated in the sin gle-shot spin readout experiments of 



Morello et aL| ( 2010 |, where the Tf 1 oc B 5 law was veri- 
fied, and the experimental values of T^ 1 were found to be 
in quantitative agreement with the prediction of Eq. [16] 
to within factors <~ 2 (see Fig. 55 'b)). The longest ob- 
served relaxation time for a single spin was T\ rj 6 s 
at B = 1.5 T. Because of the very strong field depen- 



dence of Ti, Feher and Gere (19591 were able to observe 
Ti « 5000 s at T = 1.25 K and B = 0.32 T in a bulk 
sample. 



2. Electron spin relaxation in quantum dots 

The spin relaxation mechanisms for an electron con- 
fined to a quantum dot obtained differ slightly from those 
in a donor, due to the different valley and orbital nature 
of the electron states (see section III.B and Fig. [7b. The 
ground and first excited electron wave functions are sym- 
metric or antisymmetric combination of the ±z valleys, 
because of the strong vertical confinement in the quan- 
tum well from which the dot is formed. It can be shown 
that, under this circumstance, the "valley repopulation" 



mechanism does not contribute to spin relaxation (Glavin 



and Kim 2003 Tahan 2007). 



The "one-valley" mechanism, on the contrary, is active 



and yields a relaxation rate T 1 cx B 5 ( Glavin and Kim 



2003). Notice that, unlike in the donor case, the one- 
valley mechanism in dots obtained from a [001] quantum 
well gives vanishing relaxation for B | [001] and [110] 



(Tahan 2007) 
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TABLE I Summary of the magnetic field dependence of the 
spin relaxation rates Tj~ (B), for different mechanisms appli- 
cable to donors and quantum dots in Si. It is assumed that 
the dots are formed from [001] quantum wells. 



turn dots arises from the structural inversion asymme- 
try of the quantum well in which the dot in confined, 



known as Rashba spin-orbit coupling (SOC) ( 


Hanson 


et al.\ 2007 Khaetskii and Nazarov[ |2000 ; Rashba 1960 


Tahan and Joynt 


2005 1 . Spin relaxation due to Rashba 



SOC can become dominant in Si quantum dots if the 
dot geometry gives rise to low-lying excited states of 
different orbital symmetry as compared to the ground 
state, or when the "one-valley" mechanism vanishes due 
to B || [001] or [110]. The magnetic field dependence of 
the Rashba-SOC spin relaxation channel is T-f cx B 7 



for deformation potential phonons (Hanson et al. 2007 



Raith et al. 2011 Tahan 2007), the only ones present in 
Si. The additional factor B 2 as compared to the valley- 
related mechanisms, arises from the linear dependence 
on q of the matrix element for deformation potential 
phonons to couple states of different orbital nature. Two 



experiments (Hayes et al. 2009 Xiao et al. 2010a) have 



indeed found a behavior consistent with T x oc B 7 in 



gate-defined Si dots (see Fig. |55[c,d)). Relaxation times 
as long as T% w 2.8 s at B = 1.85 T have been measured 
in a Si/SiGe dot (Simmons et al. |2011 ). The relaxation 
rate depends on the size and shape of the dot, and is 
inversely proportional to the square of the orbital level 
spacing, producing longer Ti for smaller dots. Impor- 
tantly, the Rashba SOC relaxation channel gives nonzero 
contribution for any magnetic field direction. A summary 
of the magnetic field dependencies of the spin relaxation 
mechanism for donors and dots in Si is given in Table [I] 



3. Singlet-triplet relaxation 

The Zeeman-split states of a single electron bound to 
a donor or a dot constitute a natural qubit. An alterna- 
tive scheme has been proposed, where the logical qubit is 
obtained from the two-electron singlet and triplet (S/T) 
spin states of a double quantum dot (Levy 2002). Co- 



An additional mechanism for spin relaxation in quan- herent manipulation of S/T qubits has been pioneered 
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FIG. 55 (Color online) (a) Spin-lattice relaxation rate T 1 1 
of P donors in bulk Si, at B « 0.3 T and T = 1.2 K, as a 
function of the field orientation. The angular dependence 
allows the separation of "valley repopulation" and "single 
valley" contributions. From Wil son and Feher| ( |1961[ ). (b) 
Tf for single P donors in two different devices. Both 
show a Tj -1 oc B 5 contribution, but Device A also exhibits a 
B-independent plateau, attributed to dipolar flip-flops with 
nearby donors . Also shown is T 1 _1 (3.3T) in bulk Si:P. From 



Morello et al. 



dot 



(2010) 



(c) T 1 _1 (B) in a gate-defined Si/SiGe 
• J, compared to data for a InGaAs dot (■,♦). From 
Hayes et aZ.|([2009|. (d) T X (B) in a gate-defined Si MP S dot 

From 



for the 1-electron (■) and 2-electron (o) states. 
et al] (12010a |. 



Xiao 



in GaAs double dots (Petta et al. 2005), and has re- contains 7.7% 



cently been demonstrated in Si/SiGe double dots as well 
(|Maune et al] |2012|). A 'digital' scheme applicable to 



donors in Si has also been proposed ( |Skinner et al. 2003 1. 



The first measurement of S/T relaxation in a Si double 
quantum dot ( Prance et al.\ 2012), discussed in Section 
VI. C. 4 gave a value of 2i s/t ~ 10 ms &t B = 0, two 



et al. 2005b). In an applied magnetic field, the spin life- 
time of the T_ state grows up to values of order 3 s at 
B = 1 T. The theory of S/T spin relaxation in the (1,1) 
charge configuration of double dots in Si/SiGe was dis- 
cussed in ( |Raith et aLl|2012||Wang and Wu| ~ 



two groups ( |Prada et al] |2008| |Wang et al 



2011), while 



2010) have 



analyzed the relaxation mechanisms in a two-electron 
single dot, which is relevant for the understanding of 
lifetime-enhanced transport (Shaji et al. 2008[ ) or the 



direct measurement of spin relaxation in a 2-electron dot 



(Xiao et al. 2010a). 



A theory for exchange-coupled donor pairs exists 



(Borhani and Hu 2010), and predicts a complicated de- 
pendence of the triplets-singlet relaxation rate on the 
exchange interaction J, ranging from oc J to oc J 3 . Val- 



ley interference plays a crucial role, since J can vary over 
orders of magnitude by simply changing the direction of 
the axis joining two donors. 



4. Spin decoherence 

The attractiveness of silicon for quantum computing 
applications arises in large part because of the predicted 



weakness of the decoherence mechanisms (Dc Sousa and 



Das Sarma[ [20031 |Tahan et al] [20021 |Tahan and Joynt 
2005 ) . The main source of decoherence for electron spins 



in solid state is the coupling to the bath of nuclear spins 
in the host material. In the spin resonance literature 
this goes under the name of "spectral diffusion" , to in- 
dicate that the time evolution of the state of a bath of 
nuclear spins coupled to an electron spin causes the elec- 
tron spin resonance frequency (the "spectrum") to "dif- 
fuse" over a certain range (Klauder and Anderson 1962 1. 



Natural silicon has only a 4.7% concentration of spin- 
carrying (I = 1/2) 29 Si isotope, greatly reducing the ef- 



fects of nuclear fields compared to GaAs devices ( Witzel 
and Das Sarma[ 20061. The isotopic purification to sil- 
icon consisting of only spinless 28 Si has been demon- 
strated for some time (Ager et al. 2005), and further 



pursued to extreme levels in the context of the Avogadro 



project (Andreas et al. 



2011), where a 28 Si sphere with 
less than 5 x 10 -5 2y Si concentration has been produced 
with the goal of redefining the kilogram. Purified mate- 
rial originating from the Avogadro project has been used 
to demonstrate exceptional electron spin coherence times 
T 2 



> 10 



2011). Natural germanium 



(Tyryshkin et al. 

Ya Ge (/ = 9/2), all other isotopes being 
spinless. An isotopically purified Si/SiGe heterostruc- 
ture has been demonstrated ( Sailer et al] 2009). Even 
with unenriched Ge, the effects of the Ge nuclear spins 
are greatly suppressed because in the typical heterostruc- 
tures used for Si/SiGe qubits, only about 0.6% of the 



electron density resides in the SiGe barrier (Shi et al. 



orders of magnitude longer than in GaAs dots (Johnson 2012) 



The theory of electron spin decoherence arising from 
nuclear spins in Si is well understood. The correct or- 
der of magnitude of the decoherence time T 2 was al- 
ready predicted by early studies, where the dynamics 
of the nuclear spin bath was approximated as Marko- 
vian, i.e. neglecting time correlations (De Sousa and 
Full agreement with experimental 
2003) requires a 



Das Sarma 2003) 



data in bulk Si:P (Tyryshkin et al. 



more sophisticated treatment of the spin bath dynam- 
ics, where non-Markovian time correlations are taken 
into account. Quantum mechanical solutions of the dy- 
namics of electron spins in a nuclear bath included nu- 
clear spin pair correlations (Yao et al. 2006), all the way 



to higher-order cluster-expansion techniques (Witzel and 



Das Sarma 2006 Witzel et al. 2005). The echo decay 
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takes the form: 

V(2t) ex exp(-2r/T R x cxp[-(2r/r SD ) Q ] (17) 
where r is the time interval before and after the refocus- 



ing 7r-pulse in the Hahn-echo sequence (Slichter, 19901, 
Tr is a relaxation time that accounts for both instan- 
taneous diffusion and T\ processes, T$n is the spectral 
diffusion time, which depends on the internal dynamics 
of the nuclear bath, and the exponent a takes the value 
2.3 (I Witzel et al] 120071). Full matching between clus- 



ter expansion theory and experimental data (Tyryshkin 
et al. |2003 ) (see Fig. |56|a)) was obtained by includ 



ing the Electron Spin Echo Envelope Modulation (ES- 
EEM) effect ( |Rowan et aTj|1965| ), which arises from the 



anisotropic component of the hyperfine coupling to the 



29 Si nuclei) 


Ivey and Mieher 


1975b 


Park et al. 


2009 


Saikin and Fedichkin 


2003). 



The cluster expansion technique has been extended 
to the study of decoherence upon isotopic purification 
fWitzel et al] |2010[) (Fig. [56lb)). An interesting phe- 



nomenon that must be taken into account to match Si:P 
bulk data is the interplay between 29 Si spectral diffusion 
and dipole-dipole coupling among electron spins. The 
latter gives an additional contribution to decoherence, 
which is always present because any realistic Si crystal 
contains some level of background doping. The dipolar 
coupling between donor electron spins contains terms of 
the form S^S^ + S'jfSj, which allow the excitation of 
spin 1 and de-excitation of spin 2 ( "flip-flop" process) or 
vice-versa, while conserving total energy to within the 
strength of the coupling term. This process produces 
additional magnetic noise on a spin qubit. However the 
presence of some 29 Si nuclei may actually be beneficial in 
this context. When the inhomogeneity of the local hyper- 
fine fields is stronger than the electron dipole-dipole cou- 
pling, flip-flop processes are suppressed by energy con- 
servation. The concept can be extended to any source 
of local field inhomogeneity, and the field inhomogene- 
ity does not degrade the potential implementation of ex- 



change gates (De Sousa et 
On this basis, 



al. 2001 Hu et al. 2001). 



Tyryshkin et al. (2011) have measured 



T 2 > 10 s by using a highly purified 2B Si:P crystal, and 
deliberately introducing a magnetic field gradient across 
the sample to prevent neighboring spins from undergo- 
ing energy-conserving flip-flop processes. Another way to 
suppress flip-flop is lowering the temperature such that 
g^B fceT, thereby polarizing the electron spins. This 
would lead to an exponential suppression of the dipolar 



decoherence channel (Morello et al. 2006 Witzel et al. 



2010 ) because of the scarcity of spins in an excited state. 

Once the nuclear spin and dipole-dipole decoherence 
mechanisms have been thoroughly suppressed, one may 
expect the remaining dominant decoherence channel to 
be charge noise, particularly in the case where exchange 



coupling is used to implement quantum logic gates ( Cul 



cer et al. 2009b Gamble et al. 2012). Dephasing from 
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FIG. 56 (Color online) (a) Experimental echo decay (black) 
and cluster expansion theory (orange) for nat Si:P at differ- 
ent angles of the magnetic field with respect to the crystallo- 
graphic [001] axis. Notice the echo envelope modulation aris- 
ing from anisotropic hyperfi ne coupling between donor elec- 
tron and 29 Si nuclei. From \ Witzel et al] |2007[ ). (b) Deco- 
herence time T2 for Si:P as a function of ^ a Si concentration, 
Cn, for different dopant concentrations, Cb. Symbols are ex- 
perimental data points. From (Witzel et al. 20101. 



charge noise is expected to be more pronounced in quan- 
tum dot qubits than in donor qubits, but less pronounced 
than in superconducting qubits, because the character- 
istic size of quantum dot qubits is intermediate between 
the sizes of impurity qubits and superconducting qubits. 



B. Orbital and valley relaxation 

So far we have discussed the relaxation processes for 
the electron spin confined to a donor or a gate-defined 
dot, with the intention of describing the lifetime of an 
excited qubit state encoded in the spin Hilbert space. 
Excited orbital or valley states then act as intermedi- 
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ate states for perturbations involving lattice phonons and 
spin-orbit coupling to cause spin relaxation. 

However, the orbital and valley excited states can also 
be used actively, for instance to mediate strong interac- 
tion between nearby donors. An early proposal suggested 
the use of the excited 2p orbital states of a deep donor 
to induce a superexchange interaction between pairs of 
shallow donors placed on either sides of the central one 
(Stoneham et al. 2003). This involves the coherent ma- 



relaxation times at low temperature, with a longest mea- 
sured Tie « 1.4 hours at 1.25 K and 0.3 T. These ex- 
periments were crucial in the development of the general 
theory of spin relaxation in semiconductors, as well as for 
the understanding of the electronic structure of donors 



(Fcher 19591. 



Bulk spin resonance has also been used to study elec- 
tron gases in modulation-doped Si/SiGc quantum wells 
Oantsch et all 119981 iTyryshkin et all 120051), and the 



nipulation of hydrogenic Rydberg states, a well estab- 
lished practice in atomic physics. The 2p Rydberg state 
lifetime for P donors in Si was found to be T\ ps 200 ps, 



paramagnetic defects that occur at Si/Si02 interfaces 
(Brower 1989 Poindexter and Caplan 1983) and in 



attributed to the spontaneous emission of phonons ( Vinh Biegelsen 19831. The temperature dependence of T\ 



et al. 2008). Coherent control of the Rydberg states 



amorphous silicon (Askew et al 1984 |Stutzmann and 



has also been achieved, with an orbital coherence time 
T 2 28 ps ([Greenland et aL||2010). 



in paramagnetic dangling bonds points to a relaxation 
mechanism where the electron spin is coupled to the 
charge fluctuations of the defect, which acts as a tunnel- 



Valley states are expected to have much longer life- 
times and coherence, due to the unlikelihood of processes 
that cause inter-valley transitions. A recent proposal de- 
scribes the use of singlet and triplet valley states of a 
double quantum dot to encode and manipulate quail- 



ing 2- level system (Askew et al 1984 De Sousa 2007). 



turn information with reduced sensitivity to noise ( Culcer 



et al 2012). Recent experiments have shed light on the 



valley physics and its effect on electronic states. Through 
transport spectroscopy measurements of donor states in 



Thus, ESR can be used as a non-invasive diagnostic tool 
to characterize the distribution of tunneling energies of 
defects in or near amorphous interfaces. 

The possibility of enhancing the electron spin coher- 
ence of donors in Si by reducing the concentration of the 
spin-1/2 29 Si isotope (Abe et al 2010 ) was demonstrated 
as early as 1958 (Gordon and Bowers 1958). More re- 



cently, the quality of isotopic purification has been fur- 



FinFETs, Lansbergen et al (2011) showed that under ther improved (Ager et al 2005 ) and reached a pinnacle 



certain conditions relaxation of excited states into lower 
manifolds is suppressed due to a combination of both 
spin and valley blockade. This enhanced lifetime results 
in an additional transport path through the excited state, 
and appears as a current step in the stability diagram. 
The phenomena dubbed as 'lifetime enhanced transport' 
(LET) was first observed in a silicon double quantum dot 
( |Shaji et al] |2008[ ) due to a blocked relaxation of a spin 
triplet into a ground state spin singlet, arising from the 
long spin relaxation times in silicon (see section |IV.F.2| . 
In the experiment, LET enabled [Lansbergen et aZ.| ( |201 1[ ) 
to identify a blocked transition between states that have 
different valley symmetries. They confirm this observa- 
tion (i) by extracting the tunnel rates in and out of the 
donor states through a temperature dependent measure- 
ment and analysis, and (ii) by computing the low-energy 
two-electron spectrum of the system from a multimillion 
atom tight-binding method to compare and identify the 
measured excited manifolds. 



C. Control and readout of spins in silicon 

1. Bulk spin resonance 

The dynamics of spins in bulk materials has been tra- 
ditionally studied by electron spin resonance (ESR) and 
nuclear magnetic resonance (NMR) techniques. Pioneer- 
ing experiments on ESR of donors in Si by |Feher and Gere] 
( 1959 ) measured exceptionally long electron spin-lattice 



with the Avogadro project (Becker et al. 2010), to re- 
define the kilogram as a sphere of pure ^ 8 Si. Tyryshkin 



et al (2003) showed that a 28 Si:P sample with P doping 
concentration n ~ 10 15 cm -3 exhibits a coherence time 
T'2e = 60 ms, by using a conventional Hahn-echo tech- 



nique (Slichter 1990), but accounting for the effect of 



instantaneous diffusion. In bulk experiments, the dipole- 
dipole coupling between the spins introduces an arte- 
fact whereby the refocusing pulse has the effect of flip- 
ping the coupled spins, therefore instantaneously chang- 
ing the local magnetic field and artificially suppressing 
the echo. The "true" T 2 must be obtained by extrapo- 
lating the echo decay time constant in the limit 62 — > 0. 
This extrapolation method, however, does not eliminate 



the dynamical effect of dipole-dipole coupling (Witzcl 



et al. 2010) during the wait time r. The decoherence 



due to dipolar interaction can be suppressed by intro- 
ducing a magnetic field gradient across the sample, of 
magnitude larger than the spin-spin coupling strength. 
With this method, and using a bulk sample with ex- 
treme isotopic puri ty (< 50 ppm 29 Si) and low doping 
(n ~ 10 14 cm" 3 ), |Tyryshkin et al\ fl201l| ) obtained a 
record value of T 2e ps 10 s. The combination of nar- 
row ESR absorption lines, very long spin coherence and 
the presence of a nuclear spin with / = 1/2, make the 
28 Si:P system an ideal candidate to explore sophisticated 
techniques to encode, retrieve and manipulate non-trivial 



quantum states. Morton et al (2008) demonstrated the 



ability to store and retrieve an arbitrary quantum state 
of the P electron onto the 31 P nucleus, obtaining a quan- 
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turn memory with coherence time T2„ > 1 s. 



2. Electrically-detected magnetic resonance 

In a bulk spin resonance experiment, the precession 
of a spin ensemble is detected through the electromo- 
tive force induced in a cavity or pick-up coil. In this 
way, one can only detect a macroscopic number of spins, 
typically > 10 15 . In semiconductors, however, it is possi- 
ble to make localized spins and mobile electrons coexist. 
This allows the detection of spin resonance by electri- 
cal means, and yields a significant improvement in detec- 
tion sensitivity. Electrically-detected magnetic resonance 
(EDMR) exploits spin-dependent scattering between free 
carriers and localized spins. A change in the current (or 
the conductance) of a suitably designed nanostructure 
is observed when a resonant oscillating magnetic field al- 
ters the equilibrium magnetization of localized spins onto 



which the free carriers are made to scatter (De Sousa 



et al. 2009). The free carriers can be generated cither 



by illumination (Boehme and Lips 2003) or by electro- 



statically inducing an electron layer in a MOSFET struc- 



ture ( 


van Beveren et al. 


2008 


Ghosh and Silsbee 




1992 


Lo et al. 


2007 


) . EDMR has been successfully applied to 



the detection of spin resonance and coherent control of 



31 P dopant spins in Si ( 


Huebl et al. 


2008 


Lu et al. 


2011 


Stegner et al. 


2006) (Fig. 


57 


) , and to demonstrate a very 



long-lived classical spin memory (McCamey et al. 2010). 



The detection sensitivity has been pushed to the level of 



100 donors in ion- implanted nanostructures (McCamey 



et al. 2006), and it has been proposed that reaching the 
single-spin limit is possible, and would yield a quantum 
nondemolition measurement of the donor nuclear spin 



(Sarovar et al. 2008) 



(a) 4 




(c) 




FIG. 57 (Color online) (a) Sketch of the spin-dependent tran- 
sition between a donor-bound electron and an interface trap, 
following the creation of free carriers through illumination, 
(b) Schematics of an EDMR device. P donors close to charge 
traps at the Si/Si02 interface contribute a spin-dependent 
scattering mechanism for the electrons traveling between the 
Au contacts. A resonant microwave excitation alters the po- 
larization of the donor-bound electrons, causing a measurable 
change of the overall device resistance (c) Electrically detected 
Rabi oscillations of P-donor electrons at different values of the 



3. Single-shot readout of a single electron spin 



driving power. From Stegner et al. (20061 



To reach single-spin sensitivity, it is necessary to inte- 
grate single-charge detection with a spin-dependent dis- 
placement of the charge. This idea was already incor- 
porated in the Kane proposal for a Si:P quantum com- 
puter (Kane 1998), where the readout of the electron 



spin state would take place by detecting the transfer of 
an electron from the D° state of a donor to the D~ 
state on its neighbor. This transfer is only allowed if the 
two electrons form a spin singlet state. The detection 
of the spin-dependent charge transfer would occur via 
a single-electron transistor (SET) on the surface of the 
device. Subsequent proposals pointed out that it is possi- 
ble to detect the displacement of a single charge through 



the change in conductance of a small transistor (Vrijen 



et al. . 2000 ) . In addition to the spin-dependent addition 
of a second electron to an already occupied donor, |Mar-| 
tin et al. (2003) pointed out that, in the presence of a 



large magnetic field, a charge center (not necessarily a 
donor) can change its occupancy state when the excited 
spin state lies above the Fermi level of a nearby elec- 
tron reservoir (e.g. the channel of a transistor), while the 
ground spin state lies below. This process corresponds 
to an energy-dependent spin-to-charge conversion. If a 
resonant magnetic field is applied to drive transitions be- 
tween the spin states, one expects to observe a switching 
behavior in the current through the transistor, as the sys- 
tem goes through the cycle: excite spin-up state — » ionize 
— > load spin-down electron. This method was employed 



by Xiao et al. ( 2004 ) to detect the spin resonance of a sin- 



gle charge trap coupled to a small Si transistor (Fig. 58). 
The same type of spin-to-charge conversion lies at the 
heart of the single-shot readout of a single electron con- 



fined to a GaAs quantum dot (Elzerman et al. 2004). 
In that case, the spin state of the single electron was 
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FIG. 58 (Color online) (a) Schematics of a single charge trap 
coupled to the channel of a Si transistor, (b) Single-electron 
spin resonance measurement, obtained by monitoring the av- 
erage current through the transistor as a function of mag- 
netic field, while applying a microwave excitation at 45 GHz. 
The excess current at the resonance frequency arises from the 
change in charge occupancy of the trap, made possible by the 
driven flipping of its electron spin. From Xiao et al. (2004 1. 



FIG. 59 (Color online) (a) Spin-to-charge conversion scheme 
for a single donor tunnel-coupled to the island of an SET. 
The presence of quantized states inside the SET island can 
be ignored if the single-particle ener gy level spacing is sm aller 

(b) 



Morello et al. 



(2009). 



than the thermal broadening. From '. 

Single-shot readout of a donor electron spin. The individual 
traces show the evolution of the readout signal as a function of 
the donor electrochemical potential with respect to the Fermi 
level. FromlMorello et al.\\201Q\. 



detected in a single-shot manner, i.e., with no need for 
repetition averaging, thanks to the large electrical signal 
obtained by monitoring the change in conductance of a 
quantum point contact with strong electrostatic coupling 
to the quantum dot. 

The spin-to-charge conversion, and therefore the 
single-shot spin readout, is considerably more challeng- 
ing in Si than e.g in GaAs quantum dots. This is because 
the large effective mass requires tighter electron confine- 
ment, and decreases the transparency of tunnel barriers. 
Averaged spin readout experiments were performed in 



Si/SiGe (Hayes et al. 2009) and Si MOS (Xiao et al. 



The donor and the SET island effectively form a hybrid 



double quantum dot (Huebl et al. 2010) connected "in 
parallel" (Hofmann et al. 19951, where one dot is cou- 



pled to source and drain leads, and the other (the donor 
in this case) is only coupled to the main dot ( |Morello 



et al. 2009). This results in a very compact structure and 



charge transfer signals large enough to completely switch 
the SET from Coulomb blockade (/set = 0) to the top of 
a Coulomb peak (/set ~ 2 nA), resulting in single-shot 
readout of the donor spin with > 90% visibility (Fig. 59 1 . 



2010a) quantum dots, yielding the spin relaxation time Single-shot spin readout has also been achieved in a gate- 
T\. The first successful single-shot electron spin read- defined Si quantum dot (Simmons et al. 20111, using a 



out in Si was obtained by Morello et al. (Morello et al. 
2010), where the electron was bound to a al P donor and 



tunnel-coupled to the island of an induced Si-SET. The 
readout scheme is a modification of the energy-dependent 



spin-to-charge conversion used by Elzerman et al. ( 2004 ) . 



QPC as charge sensor and the 2DEG in a Si/SiGe het- 
erostructure as the charge reservoir. There, the weaker 
(purely capacitive) coupling between sensor and dot lead 
to a current signal ~ 20 pA upon spin-dependent dis- 
placement of a single electron charge. 
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4. Readout and control of singlet-triplet states in double 
quantum dots 

Some of the most successful implementations of spin- 
based qubits in semiconductors have made use of two- 



electron systems (Levy 2002), where quantum informa- 



tion can be encoded into the singlet and triplet (S/T) 
states of exchange-coupled electrons, instead of the Zee- 
man split spin states of a single electron. Coherent con- 
trol ( Petta et al.\ |2005 ), single-shot readout (Barthel 



et al. 



et al. 



2009 ) and dynamical decoupling methods ( Bluhm 

2010) for S/T qubits have been demonstrated in 
GaAs double quantum dots. In the quest to implement 
S/T qubits in Si, the large effective mass plays again a 
role in requiring very tight electron confinement and re- 
ducing the tunnel couplings, which in this case also have 
the essential role of determining the spin exchange cou- 
pling J. In addition, most S/T qubit implementations 
in GaAs have made use of a gradient of hyperfine field, 
AB Z , between the two dots to be able to control the qubit 



along two orthogonal axes in the S/T basis (Foletti et al. 



2009 ). Because of the much smaller hyperfine interaction 



in Si (Assali et al. 2011[ ) as compared to GaAs, the two- 
axis control of a S/T qubit through J and AB Z becomes 
more challenging. On the other hand, the weak coupling 
to the nuclear spin bath allows for substantially longer 
coherence times. 

As long as the valley degeneracy of the Si conduc- 
tion band (Section III. B. 2 1 is completely lifted, the sin- 
glet /triplet spin states in Si double quantum dots can be 
detected and manipulated in the same way as in GaAs 
dots. The readout mechanism involves Pauli spin block- 



ade (Section IV.F.2). Figure 


60 shows the single-shot 


readout of the singlet and triplet states of a Si/SiGe dou- 


ble quantum dot ([Prance et al. 


20121. The state of the 


two-electron system is detected 


by pulsing the detuning e 



from negative - where the (1,1) charge state is stable - to 
positive - where the electrons can occupy the (0,2) state, 
provided their spin state is a singlet. Switching between 
the (1,1) and the (0,2) state produces a signal on the 
QPC current which can be measured in single-shot. This 
experiment also yields the triplet spin relaxation time T± , 
which is found to be ~ 10 ms for all triplets at B = 0, 
but extends up to ~ 3 s at B = 1 T for the T_ state, 
whose splitting from the singlet state is reduced by the 
applied field. 

The coherent control of singlet/triplet states has been 
achieved in an accumulation- mode Si/SiGe (Section 
IV.B.3| ) double quantum dot flMaune et~al\ |2012[ ). In 



the S/T qubit basis, one can represent the singlet and 
triplet states as the poles of a Bloch sphere, with the 



|t4), ||t) states on the equator ( |Petta et aTj pOOS). The 
exchange interaction J acts equivalently to an effective 
field along z, while a gradient of hyperfine field AB Z be- 
tween the two dots acts as an effective field along x. Fig- 
ure 161] shows the measurement of Rabi oscillations in the 
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FIG. 60 (Color online) Single-shot readout of singlet- 
triplet states in a Si/SiGe double quantum dot. (a,b) 
QPC current traces, /qpc, while pulsing the detuning with a 
square wave. Singlet states are identified when Jqpc returns 
to a high value, as in (b). (c) Charge stability diagram and 
pulsing levels, (d - f) Time traces of -Zqpc at different mag- 
netic fields, as indicated. Increasing B extends the lifetime 
of the Tn (constant current) state, (g) Control sequence, 
pulsing outside the spin blockade region. From Pran ce - et al.\ 

poT2| . 



S/T basis of a Si/SiGe DQD. The system is prepared 
in the (0,2) singlet state by exchange with the leads. A 
subsequent adiabatic pulse (i.e., slow with respect to the 
interdot tunnel rate, fast with respect to the S/T mix- 
ing time arising from hyperfine interactions) prepares a 
state close to the equator of the Bloch sphere. The ex- 
change oscillations are then initiated by pulsing closer 
to the zero detuning line, causing J to increase and the 
spin state to undergo a rotation around an axis that de- 
pends on the instantaneous value of J and AB Z . A final 
adiabatic pulse brings the state back to e > 0, where 
the electrons occupy the same dot if they returned to a 
singlet state. The oscillations of the singlet return prob- 
ability constitute a demonstration of coherent control of 
the two-electron spin states. The dephasing introduced 
by the randomness of the hyperfine field can be measured 
with a modified pulse sequence, where the electrons pre- 
pared in the (0,2) singlet state are rapidly separated and 
left to dephase at J ss before being brought back to the 
(0,2) region for readout. The measured dephasing time 



To ~ 360 ns (Maune et al. 2012 ) represents an improve- 



ment by nearly two orders of magnitude over the value 
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FIG. 61 (Color online) Coherent manipulation of 
singlet-triplet states in a Si/SiGe double quantum 
dot. (a) Charge stability diagram of the double dot system. 
Arrows describe the trajectory in gate space during the puls- 
ing sequence shown in panel (b). The (0,2) singlet state is 
prepared at point F. Adiabatically moving to point S, where 
the exchange coupling is very weak, brings the system to the 
(1,1) singlet. Pulsing to point E turns on the exchange and 
causes the two-spin state to oscillate between the (1,1) singlet 
and triplet. M is the measurement point where the electrons 
recombine in the (0,2) state if in singlet state, (c) Rabi oscil- 
lations of the singlet probability, as a function of the exchange 
pulse duration (time spent at point E) and (0,2) - (1,1) detun- 
ing e. (d) Bloch sphere representation of the trajectories of 
the two-spin states for different initial values of the hyperfine 
fields. FromlMaune et al.\ (|2012h. 



observed in GaAs dots (Petta et al. 20051, as expected 



from the much weaker hyperfine coupling in Si as com- 
pared to GaAs (Assali et al. 2011). It should be noted, 



however, that the smallness of the hyperfine field poses 
a challenge when attempting to reach the regime where 
J < AB Z (see Fig. [MJd)). 



5. Single-atom spin qubit 

The single-shot spin readout techniques discussed in 
Section |VI.C.3| have been recently combined with co- 
herent spin control via microwave pulses, yielding the 
demonstration of a spin qubit based on the electron 
bound to a single 31 P atom implanted in a silicon MOS 



device (Pla et al. 2012) (Fig. 62»). The use of energy 



selective spin-dependent tunneling as the readout pro- 



cess (Morello et al. 2010) has the consequence that the 
qubit must be operated in a magnetic field B such that 
the Zeeman splitting gfJ-sB is larger than the thermal 
broadening of the charge reservoir at electron tempera- 
ture T e i, ~ 5fc_BT e i. For instance, T c \ ~ 200 mK requires 
B > 1 T, and consequently a qubit operation frequency 
v > gfisB/h w 20 GHz. Careful microwave design is 
necessary to achieve a strong oscillating magnetic field 
B\ at the qubit, while minimizing the electric field dis- 



turbance to the charge detector (Dchollain et al. 2013) 



(Fig. [62JV)). 

Coherent control of a donor-bound electron spin was 
achieved by electrically initializing the qubit in the |4_) 
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FIG. 62 (Color online) Single-atom electron spin qubit 
based on an implanted 31 P donor, (a) Optimized design 
of an on-chip planar transmission line capable of delivering co- 
herent microwave pulses at frequencies up to 50 GHz. From 
|Dehollain et al.\ ( |2013[ ) (b). Scanning electron micrograph of 
the spin qubit device, (c) Rabi oscillations of the electron spin 
state, with 10 dBm driving power at 30 GHz. (d) Measure- 
ment of spin coherence with Hahn echo and XYXY dynamical 
decoupling, (b-d) from Pla et al. (20121. 



state, applying short microwave pulses to rotate the spin 
to and reading out the final state in single-shot. The 
probability to measure a |f) electron reveals the charac- 
teristic Rabi oscillation as function of the pulse duration 
(Fig.[62jc)). A microwave power of 10 dBm (« —20 dBm 
at the chip) resulted in a Rabi frequency of 3.3 MHz, i.e. 
an oscillating field B\ sa 0.12 mT. The coherence time of 
the electron spin qubit was measured with a Hahn echo 
technique, yielding T 2 ~ 210 /is, close to the value ob- 



tained in bulk ESR experiments in natural Si (Gordon 



and Bowers 1958). Coherent qubit rotations around two 
orthogonal axes were also demonstrate d, with an XYXY 
sequence yielding T* YXY 
(Fig.ptd)). 



410 us (Pla et al. 2012) 



VII. OUTLOOK 

We have presented a comprehensive review of the 
progress of quantum electronic devices in silicon. This 
field of research has reached maturity in the theoret- 
ical understanding and experimental demonstration of 
the confinement, manipulation and measurement of sin- 
gle electrons in nanometer-scale structures. These devel- 
opments have the potential to inform and support the 
progress of classical silicon-based nanoelectronic devices, 
as necessary to keep pursuing Moore's law in computer 
chips. Most importantly, they underpin the manipulation 
of quantum degrees of freedom such as those of electron 
and nuclear spins, and open the perspective of process- 
ing quantum information using the same technological 
platform that has allowed the extraordinary progress of 
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microelectronics over the last half-century. 

While the manipulation of individual spin qubits in 
silicon is still in its infancy, the results obtained so far 



(Maune et al. 2012 Pla et al. 20121 have confirmed 



that silicon is an excellent host for spin qubits, as ex- 
pected on the basis of the weak spin-orbit coupling and 
the abundance of isotopes with zero nuclear spin. To 
build a large-scale quantum computer in silicon, the next 
milestones will involve harnessing the exchange interac- 
tion (Section III. C. 2 1 and demonstrating 2-qubit logic 
gates, as well as coherent transport of quantum spin 
states. This is arguably the most challenging goal in the 
whole program, and several theoretical ideas have been 



put forward to achieve it (Friesen et al. 2007a Green- 



tree et al. 2004 1 . Among solid-state qubits, great success 



in multi-qubit coupling has been obtained by combin- 



ing superconducting qubits and microwave cavities (You 



and Nori 20111, and similar methods are being applied 



to semiconductor quantum dots (Frey et al. 2012 Pe 



tersson et al. 2012). Once again, silicon appears to be 



the ideal host to attain strong coupling between spins 



qubits and cavities (Hu et al. 2012), thanks to the ex- 



tremely narrow resonance lines achievable with isotopic 
purification. The ideas and methods for silicon quantum 
electronics described here lay the foundations for future 
efforts to bring the potential of silicon to full fruition in 
the quantum information era. 
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